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APRESENTAÇÃO 


A Engenharia Mecânica pode ser definida como o ramo da engenharia que aplica 
os princípios de física e ciência dos materiais para a concepção, análise, fabricação e 
manutenção de sistemas mecânicos. O aumento no interesse por essa área se dá 
principalmente pela escassez de matérias primas, a necessidade de novos materiais que 
possuam melhores características físicas e químicas e a necessidade de reaproveitamento 
dos resíduos em geral. Além disso a busca pela otimização no desenvolvimento de projetos, 
leva cada vez mais a simulação de processos, buscando uma redução de custos e de 
tempo. 

Neste livro são apresentados trabalho teóricos e práticos, relacionados a área 
de mecânica, materiais e automação, dando um panorama dos assuntos em pesquisa 
atualmente. A caracterização dos materiais é de extrema importância, visto que afeta 
diretamente aos projetos e sua execução dentro de premissas técnicas e econômicas. Pode- 
se ainda estabelecer que estas características levam a alterações quase que imediatas no 
projeto, sendo uma modificação constante na busca por melhores respostas e resultados. 

De abordagem objetiva, a obra se mostra de grande relevância para graduandos, 
alunos de pós-graduação, docentes e profissionais, apresentando temáticas e metodologias 
diversificadas, em situações reais. Sendo hoje que utilizar dos conhecimentos científicos 
de uma maneira eficaz e eficiente é um dos desafios dos novos engenheiros. 


Boa leitura. 


Henrique Ajuz Holzmann 


João Dallamuta 
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RESUMO: Discussões relacionadas ao 
desenvolvimento sustentável estão cada vez 
mais em voga buscando aliar desenvolvimento 
tecnológico com consumo consciente de bens 
e serviços. Nessa perspectiva, combustíveis 
utilizados pelos consumidores e eficiência 
energética dos veículos com esses combustíveis 
contribuem significativamente no planejamento 
de medidas mais sustentáveis tanto no mercado 
de combustíveis como no mercado automotivo. O 
exemplo dos biocombustíveis, metas de redução 
de poluentes e incentivos públicos à fabricação 
de automóveis mais eficientes concorrem para 
estimular ecossistemas sustentáveis em toda 
cadeia produtiva. Assim, a ciência de dados foi 
utilizada como ferramenta de mapeamento e 
análise do progresso de iniciativas sustentáveis, 
utilizando fases do método estatístico, incluindo 
coleta e processamento de dados e análises 
dos resultados através de medidas de tendência 


Coleção desafios das engenharias: Engenharia mecânica 2 


BRASILEIRO 


central e variabilidades, bem como análises 
gráficas. Notou-se assimetrias em consumo de 
combustível e eficiência energética, podendo 
estimular a consciência dos consumidores sobre 
o papel da tomada de decisão orientada por 
dados. 

PALAVRAS-CHAVE: Desenvolvimento 
Sustentável; Setor Automotivo; Consumo de 
Combustível; Ciência de Dados. 


EFFICIENCY AND ENERGY 
CONSUMPTION: A COMPARATIVE 
ANALYSIS IN THE BRAZILIAN 
AUTOMOTIVE SECTOR 


ABSTRACT: Discussions related to sustainable 
development are increasing, seeking to 
combine technological development with 
conscious consumption of goods and services. 
In this perspective, fuels used by consumers 
and energy efficiency of vehicles with these 
fuels significantly contribute to planning more 
sustainable measures in both fuel and automotive 
market. The example of biofuels, goals to 
reduce pollutants and public incentives for the 
manufacture of more efficient cars contribute to 
stimulate sustainable ecosystems throughout the 
entire production chain. Thus, data science was 
used as a tool for mapping and analyzing the 
progress of sustainable initiatives, using phases 
of the statistical method, including data collection, 
processing and analysis of results through 
measures of central tendency and variability, as 
well as graphical analysis. Asymmetries in energy 
efficiency and fuel consumption were noted, 
which could stimulate consumer awareness 
about the role of data-driven decision-making. 
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KEYWORDS: Sustainable Development; Automotive Sector; Fuel Consumption; Data 
Science. 


11 INTRODUÇÃO 


Para ajudar consumidores a terem um consumo mais consciente e estimular 
fabricantes de veículos a produzirem carros mais eficientes e mais econômicos, diversos 
acordos internacionais, programas governamentais e de incentivo à produção sustentável 
foram desenvolvidos no intuito de estimular esse desenvolvimento mais consciente, tanto 
no âmbito internacional como no nacional. No nacional, em especial, destacam-se o 
RenovaBio, o Inovar-Auto e o Rota 2030. 

O RenovaBio é um programa de incentivo do governo que visa estimular o cenário 
de biocombustíveis em território nacional, por meio de metas em que se busca uma maior 
introdução desses combustíveis na matriz energética brasileira, reduzir as emissões de 
poluentes provenientes dos combustíveis e atender as demandas de descarbonização 
do setor de transportes, contribuindo também para o desenvolvimento sustentável. Esse 
programa busca desenvolver a produção de biogás, biodiesel e principalmente etanol, 
como alternativa para utilização de combustíveis fósseis, como por exemplo, a gasolina. 

O Programa Inovar Auto foi desenvolvido para atuar em todo o território nacional 
durante os anos de 2013 e 2017. Visou criar condições para aumentar a competitividade no 
setor automotivo, produção de veículos mais econômicos e mais seguros e investimentos 
em empregos, pesquisa e desenvolvimento no setor. Somado a isso, as empresas que 
conseguissem cumprir as metas de eficiência energética estabelecidas, reduzindo em no 
mínimo 15,46% o consumo energético de seus veículos produzidos, tiveram reduções de 
até 2% no imposto sobre produtos industrializados (IPI). 

O programa Rota 2030, em vigor desde 2018, também tem como objetivo 
aumentar a eficiência energética dos veículos, integrá-los com as novas tecnologias na 
direção e acessibilidade, promover a capacitação de trabalhadores, investir em Pesquisa 
e Desenvolvimento, incentivar o uso de biocombustíveis e valorizar a matriz energética 
brasileira. Porém, no programa Rota 2030, ao invés de promover descontos no IPI para 
as montadoras envolvidas no programa, ele se estendeu para todo o setor automotivo, 
incluindo, por exemplo, as fábricas de autopeças e serviços estratégicos na produção de 
veículos. 

Analisando a frota de veículos, percebe-se diferenças entre as várias categorias 
de carros, assim como entre os próprios modelos desses veículos. Dessa forma, buscou- 
se analisar se esses dados disponibilizados de consumo e eficiência dos veículos estão 
cumprindo com as metas e acordos estabelecidos. Assim, por meio dos dados e das 
diferenças percebidas, propõe-se uma análise do cenário da frota de veículos nacional, 


ajudando os consumidores no abastecimento, exaltando o consumo consciente, mas 
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além disso, ajudando-os a perceberem que os carros consomem e possuem eficiências 
energéticas diferentes, o que leva também a emissões de poluentes diferentes. 


21 METODOLOGIA 


Analisou-se o desempenho dos veículos após o programa Inovar Auto com foco no 
consumo de combustíveis e eficiência energética veicular, assim como foram utilizados 
tabelas, resumos e anuários informativos, disponibilizadas por órgãos governamentais 
e agências ligadas ao setor automotivo brasileiro. As tabelas, que constituem a maior 
parte dos dados analisados, foram disponibilizadas pelo Instituto Nacional de Metrologia, 
Qualidade e Tecnologia (INMETRO). Essas tabelas contêm informações sobre os modelos 
dos veículos, seus fabricantes, suas motorizações, transmissões, combustíveis utilizados, 
suas emissões de poluentes por combustível, seus consumos de combustível na cidade e 
na estrada e suas eficiências energéticas ao longo dos anos de 2010 a 20109. 

De posse desses dados, primeiramente focou-se na análise da eficiência energética 
de cada veículo, medida em Mega Joules por Quilômetro (), analisando as variações de 
eficiência tanto entre as diferentes categorias de veículos como nos veículos de mesmo 
modelo e configurações distintas, como por exemplo, carro A com motor 1.0 e carro A com 
motor 1.6. Essa eficiência energética se relaciona também com o Programa Brasileiro de 
Etiquetagem (PBE), que desde 2009 afixa nos diversos veículos uma etiqueta, a Etiqueta 
Nacional de Conservação de Energia (ENCE), informando ao consumidor o consumo 
energético de cada modelo, vide figura 1. 


Energia (contusiver 
mo 0000 





Figura 1: Etiqueta Nacional de Conservação de Energia (ENCE) veicular. 


Fonte: INMETRO. 
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Outra medida das tabelas disponibilizadas pelo INMETRO utilizada para comparar 
os veículos foi o consumo de combustíveis tanto na cidade como na estrada, de forma que 
posteriormente poderia ser criado uma razão entre esses dois consumos. 

Os consumos de etanol e gasolina na cidade, em Quilômetros por Litro (), também 
foram analisados, servindo de comparativo entre os rendimentos dos dois combustíveis. 
Após essa comparação, seria feita uma razão entre esses consumos de forma a obter 
através de uma única medida uma relação entre os dois, esquematizado na figura 1. 


o ANMETRO 


E castosnercénco 











INFORMAÇÕES SOBRE VEÍCULOS mm —s ANÁLISE ESTATÍTICA mma RESULTADOS 


INOVAR AUTO 
A qd —— 


Figura 1: Esquema de Processamento e Análise de Dados em Estudo 


RENOVABIO 


Fonte: Própria. 


Dessa forma, com esses dados, busca-se comparar os consumos e as eficiências 
dos veículos com os mais diferentes tipos de combustíveis, ajudando o consumidor a 
perceber essas diferenças e ajudando-o escolher a melhor opção para seu veículo. 


31 RESULTADOS E DISCUSSÃO 


Após a data limite para atingimento das metas do Inovar Auto, nem todas as montadoras 
conseguiram reduzir os consumos e aumentar as eficiências de seus veículos. Porém, os 
consumos que estavam seguindo uma tendência decrescente durante todo o período do 
programa, apresentaram aumento logo após seu término, podendo ser um indicador de que 
as montadoras pretendiam somente atingir as metas estabelecidas de reduzir no mínimo 
15,46% do consumo energético para obter os descontos. 

A partir dos dados obtidos nas tabelas, pode-se perceber que houve uma grande 
variação no consumo ao longo dos anos. Essas variações se propagaram não somente entre 
as categorias de veículos, mas também entre os modelos analisados, exemplificados pelas 
medidas estatísticas da tabela 1. 
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Sub Extra Utilitário Fora de 


Categoria Compacto Médio Grande Minivan Comercial Picape 


Compacto Grande Esportivo Estrada 
Média 1,56 1,69 1,70 1,95 2,14 2,04 2,97 2,08 2,05 2,42 
Mediana 1,52 1,69 1,70 1,92 2,14 1,99 2,97 2,00 2,04 2,53 
Desvio 
Padrão 0,12 0,12 0,13 0,15 0,12 0,28 0,21 0,25 0,17 0,42 


Tabela 1: Consumo Energético por Categoria de Veículo em 2019 (em) 


Fonte: Própria. 


Percebeu-se um comportamento semelhante com a eficiência energética, a qual 
também apresentava variações tanto entre as categorias como nos modelos, exemplificado 
pelas variações no Desvio Padrão na tabela 1, com 0,12 para compactos e 0,42 para picapes, 
por exemplo. Os dados revelaram que essas eficiências são proporcionais as categorias 
veiculares, de forma que as categorias de veículos maiores — maiores áreas - apresentaram 
maiores gastos energéticos. Além disso, as variações nos modelos com características de 
fabricação diferentes mostraram que, mesmo sendo de uma categoria de veículos específica, 
certos modelos apresentaram variações significativas em relação a seus modelos similares, a 
ponto de apresentarem consumos energéticos como se fossem de outra categoria, mostrado 
na figura 2. 


SUB COMPACTO -— E 


COMPACTO ——1TL LI ' 
MÉDIO -—— LE , 
GRANDE -— TE , 
E! 
x EXTRA GRANDE = mm 
mr 
E UTILITÁRIO ESPORTIVO ' [— [1 [++ + : 
FORA DE ESTRADA -—BHB - 
MINIVAN -— TE ———: 
PICAPE ——— TT — ———————— 


1.50 1.75 2.00 2.25 2.50 2.75 3.00 
CONSUMO ENERGÉTICO (Mj/km) 


Figura 2: Eficiência Energética por Categoria de Veículo em 2019. 


Fonte: Própria. 


Exemplificando melhor esse conceito de veículos similares com consumos distintos, 
temos a tabela 2, a qual mostra as diferentes configurações que um mesmo modelo de 
veículo possui. 
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EFICIÊNCIA RAZÃO 
MODELO VERSÃO MOTOR CÂMBIO ADE CAE ENERGÉTICA ETANOL- 
( GASOLINA 
DRIVE (COM 
ARGO START & 1.0 M-5 9,9 14,2 1,45 0,69 
STOP) 
ARGO 1.0 1.0 M-5 9,1 13,0 1,55 0,70 
ARGO DRIVE 1.0 M-5 9,1 13,0 1,55 0,70 
DRIVE (COM 
ARGO START & 1.3 M-5 9,2 12,9 1,57 0,71 
STOP) 
ARGO DRIVE 1.3 M-5 87 12,6 1,60 0,69 
ARGO DRIVE GSR 1.3 MTA-5 8,9 12,7 1,60 0,70 
ARGO HGT 1.8 A-6 6,9 10,0 1,96 0,69 
ARGO HGT 1.8 M-5 78 14 1,79 0,68 
ARGO PRECISION 1.8 A-6 71 10,1 1,91 0,70 
ARGO PRECISION 1.8 M-5 80 1,5 1,74 0,69 
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Tabela 2: Veículo Fiat Argo com diferentes características no ano de 2019 


Fonte: Própria 


Nessa tabela 2, além das diferentes versões de Fiat Argo do ano de 2019, apresenta-se 
o câmbio que esses veículos apresentam, sendo eles manual de cinco marchas, automático 
de 6 marchas e automatizado de 5 marchas; o rendimento desses veículos tanto com etanol 
quanto com gasolina na cidade; a eficiência energética que esses veículos possuem e, por 
fim, uma razão entre os rendimentos, que por meio de porcentagem mostra o quanto o 
rendimento de etanol representa do rendimento da gasolina. 

Ainda analisando a tabela 2, percebe-se a variação que veículos de mesmo modelo e 
versões diferentes podem possuir devido a características de fábrica distintas. Como exemplo, 
temos as versões Argo Drive com Start & Stop 1.0 e Argo 1.0, que apresentaram eficiências 
de 1,45 e 1,55, respectivamente, além de apresentarem também diferentes rendimentos 
com etanol e com gasolina. Assim como as versões HGT 1.8 Automático de 6 marchas e HGT 
1.8 Manual de 5 marchas, que apresentando esses diferentes tipos de câmbio implicaram em 
diferentes eficiências energéticas, 1,96 e 1,79, respectivamente. 

De maneira geral, percebe-se que os veículos possuem diferentes características e 
atributos de fábricas atribuídas a suas versões, atributos esses que podem ser motorização 
e câmbio, que implicam em maiores e menores consumos de combustível e eficiência 
energética, dependendo da versão. Mas além disso, até naqueles veículos em que esses 
atributos de fábrica são iguais, sistemas desenvolvidos e tecnologias embarcadas, como o 
caso do Start & Stop, pode levar a variações nesses parâmetros estudados. 
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41 CONSIDERAÇÕES FINAIS 


Por meio do estudo comparativo pode-se inferir que os carros apresentam 
variabilidades importantes, consumindo combustíveis de maneiras diferentes, com eficiências 
energéticas assimétricas. Assim, não é possível generalizar consumos e eficiências dentro 
das categorias visto essas variações e amplitudes analisadas. Uma vez que as categorias, 
incluem diversos modelos, podemos inferir sobre a variabilidade entre os modelos para os 
parâmetros avaliados. 

Dessa forma, esse estudo visa ajudar os consumidores a perceberem essas 
diferenças e compreenderem esse panorama complexo dentro do setor automotivo com 
parâmetros variados. Além de ajuda-los esquematizando todas essas ramificações possíveis 
dentro desse setor, fazendo-os perceberem os dados e os comportamentos de seus veículos 
específicos. 

Pretende-se dar continuidade ao estudo buscando ainda mais parâmetros sobre os 
veículos e seus diferentes comportamentos, como por exemplo poluentes emitidos, em prol 


de uma análise cada vez mais completa sobre os veículos e sobre a frota brasileira. 
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RESUMO: Alternativas para produção de energia 
limpa e redução no consumo de combustível 
são constantemente estudadas, com o intuito 
de aprimorar a eficiência energética e a 
diminuição na emissão de poluentes. Portanto, 
foi desenvolvido um estudo de rendimento na 
geração de hidrogênio para um automóvel, 
a qual será obtida através da eletrólise. O 
objetivo deste trabalho é verificar quais tipos 
de soluções apresentam melhor eficiência na 
geração do hidrogênio, mediante testes com 
conjunto comercial automotivo e análise de 
dados experimentais. Com os resultados obtidos, 
pôde-se constatar que tanto a água residencial 
quanto a mineral, adicionadas de hidróxido de 
potássio, apresentaram bom fluxo de gás através 
da medição do fluxômetro. A vazão mensurada 
no sistema corresponde aos requisitos mínimos 
para redução no consumo de combustível e na 
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liberação de gases poluentes. 
PALAVRAS-CHAVE: Geração de hidrogênio, 
energia limpa, eletrólise, eletrólitos, tipos de 
água. 


COMPARATIVE ANALYSIS OF WATER 
AND ELECTROLYTE TYPES OF A 
HYDROGEN GENERATOR FOR A 

CONVENTIONAL VEHICLE 


ABSTRACT: Alternatives are being studied for 
the production of clean energy and reduction in 
fuel consumption, with the objective of improving 
energy efficiency and decrease the emission 
of pollutants. Therefore, a study of the yield in 
the generation of hydrogen for an automobile 
was developed, which will be obtained through 
electrolysis. The objective of this work is to 
verify which solutions present better efficiency 
in the generation of hydrogen, through tests 
with a commercial automotive kit and analysis of 
experimental data. With the results obtained, it 
was possible to verify that residential and mineral 
water, added to potassium hydroxide, presented 
good gas flow in the flowmeter measurement. 
The system's flow rate meets the minimum 
requirements for reducing fuel consumption and 
decreasing polluting gases. 

KEYWORDS: Hydrogen generation, clean 
energy, electrolysis, electrolytes, water types. 


11 INTRODUÇÃO 
É de suma importância a preocupação 
com questões relacionadas à sustentabilidade e 


que não poluam o meio ambiente. Os automóveis, 
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através da combustão, aumentam a concentração de gases que provocam o efeito estufa. 
Com esse fator, surge a necessidade de se desenvolver formas de obter energia limpa e 
que não prejudicam o meio ambiente. Uma das opções é através da geração de hidrogênio 
para veículos convencionais, que diminuirá a injeção de combustível feita pelo motor e 
consequentemente reduzirá os gases poluentes emitidos. 

A geração de hidrogênio é atualmente uma forma alternativa de energia em grande 
desenvolvimento (Glenk e Reichelstein, 2019) e, devido a esse fato, têm-se a importância 
de estudos relacionados à sua geração e dos meios para sua utilização (Ilbrahin, 2014). 
Mundialmente, a produção anual de hidrogênio é de cerca de 5x10!! Nmº, a qual corresponde 
a 2% da demanda de energia primária (Veras, 2015). Somado a esses fatores, de acordo 
com o Centro de Gestão e Estudos Estratégicos — CGEE (2010), o Brasil, por exemplo, 
investiu em tecnologias do hidrogênio, entre 1997 e 2007, cerca de R$ 134 milhões. 

Ademais, Wollmann (2013) ressalta que o hidrogênio possui a maior taxa de energia 
por unidade de massa se comparado com demais combustíveis, o que significa que sua 
queima possui um grande aproveitamento do seu poder calorífico. 

A eletrólise, juntamente com a reforma a vapor do etanol, reforma do gás natural e 
a gaseificação da biomassa, é considerada como uma das principais formas de obtenção 
do hidrogênio (Veras, 2015; Anantharaj et al., 2018). Aplicada neste trabalho, a eletrólise 
da água consiste na quebra da molécula de hidrogênio em dois elementos: H* e OH (Fig. 
1). O processo ocorre por meio da passagem de uma corrente elétrica através da solução 
aquosa, ionizando as partículas. Esta solução, por si só, não permite a passagem elétrica 
com facilidade, sendo necessário adicionar um eletrólito para a realização da quebra da 


molécula de forma mais eficiente. 


Fonte de energia e solução iônica 


Bateria do automóvel (energia elétrica) H,0 (água) e eletrólito 


TS 


Processo de eletrólise 





Separação dos elementos: H* e OH 


Yw 


Queima de hidrogênio em conjunto com o combustivel 





Figura 1 — Geração de hidrogênio para um veículo convencional através da eletrólise da água. 


Asolução com substâncias diluídas na água, sejam sais, ácidos ou bases, é denominada 
solução iônica. Contudo, pode-se encontrar eletrólitos sólidos e também gasosos (Marques, 
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2015). A água com um sal solúvel, por exemplo, tendo a passagem da corrente elétrica, 
permite a quebra da ligação química entre os átomos de hidrogênio e oxigênio, separando- 
os. À separação atômica produz íons, que são partículas carregadas eletricamente. 

A geração de hidrogênio através da eletrólise, para ser utilizada nos veículos, surgiu 
como um tema a ser estudado e analisada a sua viabilidade para redução da emissão de 
poluentes e, por conseguinte, diminuir a poluição do ar. 

Com os fatores apresentados, tem-se a possibilidade de redução no consumo de 
combustível, utilizando-se da queima e utilização imediata do gás, não sendo necessário um 
tanque de gás como reservatório e contribuindo para a segurança do sistema. 

Através das diferentes soluções para realização da eletrólise, neste trabalho será 
analisado qual o melhor composto para geração de hidrogênio, juntamente com os tipos de 


água. 
21 METODOLOGIA 


As medições de alta condutividade são realizadas quando as amostras são 
soluções aquosas, sendo amplamente concentradas em sais, ácidos ou bases (Radiometer 
Analytical SAS, 2004). Além disso, a condutividade depende parcialmente do pH (medida da 
concentração de íons H+ na água), da temperatura e da quantidade de dióxido de carbono 
atmosférico, que foi dissolvido na água para formar íons. 

Entre os cinco tipos analisados, a que apresenta o mais alto índice de condutividade 
é a água salgada, por conter em sua composição química 3,5% de sais. O valor de 
condutividade deste tipo de água é de 42.896 mS/cm (micro-Siemens por centímetros), em 
uma temperatura de 15 ºC (Radiometer Analytical SAS, 2004). Os demais tipos de água que 
foram analisados são: deionizada, residencial, destilada e mineral. 

O conjunto de geração de hidrogênio automotivo foi preparado em laboratório para 
análise de seu rendimento. Como apresenta a Figura 2, o conjunto possui uma célula de 
geração de hidrogênio, um reservatório e tubos para o transporte do líquido e do gás. 
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Entrada de água 







Reservatório | | Saída do ar 


Mangueira 


Célula de geração de 
hidrogênio 


Figura 2 — Conjunto comercial automotivo de geração de hidrogênio. 


Para a realização do experimento, utilizou-se 500 ml de água no reservatório 
e 500 ml para dissolver a substância em outro utensílio. Inicialmente, foram colocadas 7 
g da substância na água. Após a dissolução, a mistura foi colocada no reservatório para 
funcionamento da célula de geração. Posteriormente, acrescentou-se mais 5 g da substância, 
a fins de comparação. 

Como a eletrólise depende da passagem de corrente elétrica, utilizou-se uma fonte de 
energia de 12 V de corrente contínua para simular a bateria de um veículo convencional. A 
quantidade de energia elétrica que passa pelo sistema foi aferida através de um multímetro. 

A quebra de molécula da água (H,0), obtida com a passagem da corrente elétrica, 
separa os átomos de H* e OH. Após a quebra, o gás passa por um filtro para retirada de 
possíveis impurezas e em seguida pelo fluxômetro, indicando o nível de vazão do gás. 

O fluxômetro (PROTEC) utilizado para medir o fluxo de ar é aprovado pelo INMETRO 
(Instituto Nacional de Metrologia, Qualidade e Tecnologia). Esse dispositivo possui válvula 
para regulagem, esfera flutuadora de aço inoxidável e escala graduada de O a 15 I/min. 


31 RESULTADOS E DISCUSSÃO 


3.1 Análise de diferentes tipos de água 


Para a verificação do rendimento na produção de hidrogênio em relação aos 


Coleção desafios das engenharias: Engenharia mecânica 2 Capítulo 2 o 


diferentes tipos de água, os testes em laboratório apresentaram os resultados dados a seguir. 
Em primeiro instante, no teste com água residencial, foram utilizadas 7 g de hidróxido de 
potássio. A temperatura inicial da água era de 32,0 ºC. Com o funcionamento do conjunto 
experimental, a corrente elétrica apresentou 4,9 A, enquanto o nível de vazão variou entre 
2 e 3 I/min. A temperatura final no reservatório estava em 33,0 ºC. Logo em sequência, na 
segunda parte do teste com água residencial, acrescentou-se 5 g de hidróxido de potássio, 
resultando assim em 12 g utilizadas. Foram feitas as seguintes medições: 5,5 A de corrente 
elétrica e nível de vazão entre 3,5 a 4 I/min. A temperatura final do sistema foi de 33,5 ºC. O 


resumo dos dados é apresentado na Tabela 1. 


Temperatura Substância e sua Corrente elétrica Temperatura no Nível de vazão 
inicial (ºC) quantidade (A) reservatório (ºC) (m3x10-5/s) 





potássio 


de potássio 
Tabela 1 — Resultados dos testes realizados com água residencial. 
A água mineral foi o segundo tipo utilizado no experimento, com temperatura de 
27,1 ºC. Para 7 g iniciais de hidróxido de potássio, a corrente elétrica do sistema foi de 4,3 
A, enquanto a temperatura final indicou 28,9 ºC. O nível de vazão, através do fluxômetro, 
apresentou valor de 4,5 I/min. 
Para 12 g de hidróxido de potássio, a água mineral obteve aumento nos valores de 


corrente elétrica, temperatura e vazão. A Tabela 2 demonstra os resultados obtidos: 


Temperatura inicial | Substância e sua | Corrente elétrica | Temperatura no Nível de vazão 
(ºC) quantidade (A) reservatório (ºC) (m3x10-5/s) 


27,1ºC Fase aos ASA 28,9 ºC 7,50 
potássio 
12 g de hidróxido 5 


Tabela 2 — Resultados dos testes realizados com água mineral. 





O terceiro teste foi realizado com água salgada em três instantes: No primeiro 
momento, não foi acrescido hidróxido de potássio na água, devido aos sais presentes na 
sua composição. Dessa forma, apresentou corrente elétrica de 1,6 A, temperatura final em 
23,1 ºC e nível de vazão de 4 I/min. 


Após essa etapa, no segundo instante acrescentou-se uma quantidade de 7 g de 
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hidróxido de potássio. Os resultados obtidos foram os seguintes: a corrente elétrica que 
passava no sistema, medida por um multímetro, foi de 3,5 A, enquanto que a temperatura 
aumentou para 24,3 ºC. O nível de vazão resultou entre 4 e 4,5 I/min. 

A Tabela 3 apresenta os dados do experimento: 


Temperatura 


Ed Nível de vazão 
no reservatório 


(mºx10*/s) 


Temperatura Substância e sua Corrente 
inicial (ºC) quantidade elétrica (A) 





12 g de hidróxido 
os de potássio E e dia 


Tabela 3 — Resultados dos testes realizados com água salgada. 


Para água destilada com 19,2 ºC, os dados do experimento são apresentados da 
seguinte forma na Tabela 4: 


Temperatura 


E Nível de vazão 
no 8 (m3x10/s) 


potássio 
12 g de hidróxido 
de potássio 


Tabela 4 — Resultados dos testes realizados com água destilada. 


Temperatura Substância e sua Corrente 
inicial (ºC) quantidade elétrica (A) 





O último teste, utilizando-se água deionizada à 18,5º C, indicou corrente de 4,9 A 
no primeiro instante do funcionamento do conjunto. A temperatura final no reservatório foi 
medida em 27,2 ºC e nível de vazão através do fluxômetro apresentou cerca de 2,5 a 3 |/ 
min. A Tabela 5 demonstra os resultados obtidos do experimento com 7 e 12 g de substância. 


Temperatura 


su Nível de vazão 
no reservatório 


(m3x10*/s) 


Temperatura | Substância e sua Corrente 
inicial (ºC) quantidade elétrica (A) 


7 g de hidróxido 
de potássio 





pe 
27,2 Pegando 29,1 5,00 a 5,83 
de potássio 


Tabela 5 — Resultados dos testes realizados com água deionizada. 
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3.2 Análise de substâncias e temperaturas 


A próxima etapa refere-se ao estudo das substâncias químicas, os eletrólitos, para a 
realização da eletrólise na água. Conforme Chang e Goldsby (2013), a solubilidade pode ser 
definida como a quantidade máxima de soluto a ser dissolvida em um solvente, a determinada 
temperatura. A Tabela 6 apresenta a solubilidade (em 9/1) de substâncias consideradas para 
a realização dos experimentos. 


Substância | Solubiidado — 


Cloreto de sódio 358 gl a 20 ºC 


Tabela 6 — Substâncias químicas e respectiva solubilidade. 





Fonte: Dados de IFA (Institute for Occupational Safety and Health of the German Social Accident 
Insurance), 2021. 


Por fim, foram selecionados conforme facilidade de obtenção, custo e solubilidade: 
sulfato de cobre, hidróxido de sódio, carbonato de potássio e hidróxido de potássio. 

Verificou-se os resultados para sulfato de cobre dissolvido em água residencial, o qual 
resultou em um fluxo nulo de gás, seja para 7 ou 12 g da substância (Tabela 7). 


ir Temperatura ; = 
Temperatur nci rren e Nível de vazão 
emperatura Substância e sua Corrente no reservatório 


inicial (ºC) quantidade elétrica (A) (ºC) (mºx10-5/s) 


7 g de sulfato de 
cobre 

12 g de sulfato de 
cobre 


Tabela 7 — Resultados dos testes realizados com sulfato de cobre e água residencial. 





No segundo momento foi utilizado hidróxido de sódio para o experimento com água 
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residencial. A Tabela 8 apresenta os resultados obtidos. 


Temperatura Substância e sua Corrente Temperatura no Nível de vazão 
inicial (ºC) quantidade elétrica (A) reservatório (ºC) (mºx10-5/s) 
7 g de hidróxido de 
e 


Tabela 8 — Resultados dos testes realizados com hidróxido de sódio e água residencial. 





Para o teste com carbonato de potássio, 7 g da substância apontou fluxo insignificante 
de gás, enquanto que com 12 g o fluxômetro apresentou entre 0,5 a 0,8 I/min. Pode-se 
observar pela Tabela 9 que o máximo de corrente elétrica utilizada no sistema foi de 3,3 A e 


o ganho de temperatura também foi baixo. 


Temperatura | Substância e sua Corrente Temperatura no Nível de vazão 
inicial (ºC) quantidade elétrica (A) reservatório (ºC) (mºx10%/s) 
7 g de carbonato de 
12 g de carbonato 


Tabela 9 — Resultados dos testes realizados com carbonato de potássio e água residencial. 





Após a realização do experimento com água residencial, foi feita a análise das 
substâncias quando diluídas em água mineral. Foram utilizados carbonato de potássio e 
hidróxido de sódio. 

O primeiro teste com água mineral e 7 g de carbonato de potássio apresentou vazão 
quase nula, enquanto que cerca de 40 g da substância foram necessárias para gerar um fluxo 
de 1,5 I/min (Tabela 10). 


Temperatura Substância e sua | Corrente elétrica Temperatura no Nível de vazão 
inicial (ºC) quantidade (A) reservatório (ºC) (m3x10-5/s) 


7 g de carbonato de 

26,4 potássio 1,8 26,5 Quase nula 
40 g de carbonato de 

26,5 potássio 3,5 30,1 2,50 


Tabela 10 — Resultados dos testes realizados com carbonato de potássio e água mineral. 
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A Tabela 11 apresenta os resultados obtidos para a diluição de hidróxido de sódio em 
àgua mineral. Verificou-se que o nível de vazão de gás foi semelhante à sua utilização com 


àgua residencial, de no máximo 2 I/min. 


Temperatura Substância e sua Corrente Temperatura no Nível de vazão 
inicial (ºC) quantidade elétrica (A) reservatório (ºC) (mºx10-5/s) 


25.8 7 g de hidróxido de 39 30,6 2509 3,33 
sódio 
12 g de hidróxido 
26,5 SR 48 321 250 


Tabela 11 — Resultados dos testes realizados com hidróxido de sódio e água mineral. 





Por fim, foi realizado um teste com água residencial e 12 g de hidróxido de potássio. 
Durante 50 minutos o conjunto esteve em funcionamento para serem analisados os 
parâmetros de temperatura, corrente elétrica e fluxo de gás. 

A temperatura inicial da água era de 32,1 “C. Com 10 minutos de operação, subiu 
para 34,1 ºC. Após o início, a temperatura teve pouca variação até o final do teste. Com 50 
minutos, o reservatório do sistema estava com 35,4 ºC, tendo uma variação de apenas 3,3 
ºC entre temperatura inicial e final, de acordo com a Fig. 3. 


Temperatura - em “C 





35.9 
35.0 
» 34,5 
du 
ai 
[= 
- 34.0 
- 
33.5 
33.0 


10 minutos 20 minutos 30 minutos 40 minutos 50 minutos 


Figura 3 — Resultado do teste realizado com hidróxido de potássio e água residencial. 
Com 10 minutos de funcionamento do conjunto, o sistema apresentou corrente 


elétrica de 4,6 A. Após 40 minutos, a corrente baixou para 4,2 A e se manteve constante 
até o final do teste (Fig. 4). 
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Corrente - em A (amperes) 


4.8 


4.6 


4.4 


Corrente elétrica 


4.2 





2.0 
10 minutos 20 minutos 30 minutos 40 minutos 50 minutos 


Figura 4 — Resultado do teste realizado com hidróxido de potássio e água residencial. 


Correspondendo aos testes realizados anteriormente com água residencial e hidróxido 
de potássio, o fluxômetro apresentou cerca de 5 I/min com 10 minutos de funcionamento 
do sistema. Com 20 minutos, a vazão baixou para 4,5 I/min e continuou constante até 40 
minutos do experimento. A última medição mostrou que o fluxo de gás estava com 4 I/min. 


A Fig. 5 apresenta os resultados obtidos: 


Fluxo de gás - em m*(10"")/s 





Fluxo de gás 
| 


6,66 


10 minutos 20 minutos 30 minutos 40 minutos S0 minutos 


Figura 5 — Resultado do teste realizado com hidróxido de potássio e água residencial. 


A quantidade ideal de eletrólito quanto ao hidróxido de potássio foi determinada 
considerando uma média aceitável de vazão em m/s, visando não exceder o limite de 


corrente elétrica suportada pela fonte de energia. 


41 CONCLUSÃO 


A realização dos testes utilizando conjunto convencional automotivo permitiu 


a identificação dos melhores tipos de água e eletrólitos selecionados para a geração de 
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hidrogênio, conforme objetivo da pesquisa. Pode-se observar com os experimentos que a 
água mineral obteve um bom fluxo de gás produzido, juntamente com a água residencial. É 
possível concluir, por hipótese, que a condutividade elétrica e o fluxo de gás produzido têm 
relação com os minerais presentes na água no processo de eletrólise. A água destilada, por 
exemplo, que não possui sais em sua composição, apresentou baixo nível de vazão. 

A água residencial destaca-se por ser de fácil obtenção. A água salgada, apesar de 
apresentar bom fluxo de gás nos testes, é de difícil aquisição, deixa resíduos nos reservatórios 
e pode ocasionar uma possível corrosão no sistema conforme seu uso constante. 

Em relação aos eletrólitos, verificou-se que o hidróxido de potássio possui boa 
solubilidade e auxilia na passagem de corrente elétrica. A substância apresentou maior fluxo 
de gás se comparado aos demais eletrólitos. 
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ABSTRACT: The purpose is to review iterative 
and series methods applied to the solution of 
Kepler's equation, which is solved over the entire 
range of elliptic motion. The method whose results 
will work as a reference is the Newton-Raphson's 
numerical method. The results will be discussed 
around the number of iterations required until 
the convergence criterion is satisfied, that is, 
residual error in eccentric anomaly lower than 
rad (for iterative methods) or rad (for series- 
based methods) and the processing time. The 
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FOR KEPLER'S EQUATION 


advantages and drawbacks of each method will 
be presented. 

KEYWORDS: Kepler's Equation; Numerical 
Methods; Iterative Solutions; Series Solutions. 


11 INTRODUCTION 


Artificial satellites are employed in many 
activities, such as space exploration, land 
mapping, microgravity experiments and 
telecommunication. Regardless of the mission 
for which the satellite is designed, the knowledge 
of its most accurate possible orbital position is 
critical for the mission success. And here lies 
the importance of Kepler's equation: it gives a 
relation between the position of the satellite and 
time (Battin, 1999). 

The elliptical form of Kepler's Equation is given 
by: 


M=E-—esín (E) (1) 


where the three quantities are related to the 
orbit Keplerian elements: M is the mean 
anomaly; E, the eccentric anomaly; and , the 
eccentricity. Kepler's equation is transcendental 
in E; therefore, the solution for this quantity, 
when M is given cannot be expressed by a finite 
number of terms (Battin, 1999). The solutions 
for Kepler's equation can only be approximated, 
generally using computational methods. 

Since Kepler's equation is one of the most 
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famous transcendental equations, it has inspired many developments in mathematics during 
the last decades. 

An algorithm based on simple initial cubic approximations and a slight generalization 
of the Newton-Raphson method, was presented for the solution of Kepler's equation (Ng, 
1979). In Danby and Burkardt (1983), methods of iteration are discussed in relation to Kepler's 
equation, considering various initial “guesses”, with possible strategies for their choices. 
Several of these iterative methods are compared; the one used in the comparisons has 
local convergence of fourth order. If in the first study, Danby and Burkardt considered the 
solution of the conventional form of Kepler's equation for elliptic orbits, in the second, they first 
considered hyperbolic orbits equation, then generalizations for elliptic and hyperbolic orbits 
(Burkardt and Danby, 1983). 

Serafin (1986) analytically examined techniques for selecting the interval within which 
the root of the Kepler's equation of satellite motion is to be sought. In 1986, Odell and Gooding 
reviewed starting formulas and iteration processes for the solution of Kepler's equation, giving 
details of two complete procedures that operates with an iterative process. Mikkola (1987) 
derived a method to obtain an approximate solution for Kepler's equation that could be used 
for all orbit types, including hyperbolic. 

Markley (1996) solves Kepler's Equation over the entire range of elliptic motion by 
a fifth-order refinement of the solution of a cubic equation. This method requires a square 
root, a cube root, and two trigonometric functions (four transcendental function evaluations). 
In Fukushima (1997), two approximations of the Newton-Raphson method were developed. 
The firstis a sort of discretization, namely to search an approximate solution on pre-specified 
grid points. The second is a Taylor series expansion. A combination of these was applied to 
solving Kepler's equation for the elliptic case. Later, he developed a procedure to solve a 
modification of the standard form of the universal Kepler's equation, which is expressed as a 
nondimensional equation with respect to a nondimensional variable (Fukushima, 1999). 

Condurache and Martinusi (2007) present an exact vectorial solution to the 
Kepler problem. A vectorial regularization linearizes Kepler's equation, using a Sundman 
transformation. A unified approach to the classic Kepler problem is offered, by studying both 
rectilinear and non-rectilinhear Keplerian motions with the same instrument. In Davis et al. 
(2010) seven sequential starter values for solving Kepler's equation for fast orbit propagation 
are proposed. These methods have constant complexity (not iterative), do not require pre- 
computed data, and can be implemented in a few lines of code. 

More recently, Reza and Ghadiri (2014) focused on Newton-Raphson's method for 
solving Kepler's equation. In order to increase the stability of Newton's method, various 
guesses were studied. Based on time of implementation, an appropriate choice is presented: 
first guesses that increase the isotropy and decrease the solution time of implementation. 
Starting algorithms for the iterative solution of elliptic Kepler's equation are also considered 
in Calvo et al. (2013), where new global efficiency measures are introduced and several well- 
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known starters with minimum computational cost are analyzed on the light of these efficiency 
measures. And Avendano et al. (2015) used Smale's a-theory to prove that Newton's method 
starting at the defined approximate zero produces a sequence that converges to the actual 


solution at quadratic speed. 


21 ITERATIVE AND SERIES-BASED METHODS 


In this paper, the approximations of six different methods are compared: four are 
iterative methods, and two, methods based on series approach. The iterative methods are: 
Newton-Raphson (the solutions it produced were used as a reference), Halley, Regula-Falsi 
and Successive Approximations, and were computationally implemented. The methods 
based on series solutions comprise Lagrange Expansion Theorem and Fourier-Bessel Series 
Expansion (Battin, 1999). 


2.1 Newton-Raphson's method 


Newton-Raphson's method is an iterative method usually applied to numerical solution 
of equations of the form , where is differentiable. The iteration function of this method is 
(Franco, 2006): 


o f(Xn) 
Xn+1 = Xn — FO) (2) 


Kepler's Equation can be conveniently written as: 


E-esin(E)—M=f(E) (3) 


and solved by Newton-Raphson's Method via: 


En— esin(E,)— M 


E. =E—--———— 
dire 1 -ecos(E,) (4) 


The iteration process in Eq. (4) stops when the root accuracy reaches a specific value, 


determined by each problem. 


2.2 Halley's method 


Halley's Method is a generalization of Newton's method that aims at finding the root 
of a nonlinear equation and requires analytical and numerical computation of higher-order 
derivatives of the function. The algorithm adapted to Kepler's equation for any fixed value of 
to iterate for is given by (Gander, 1985): 


Coleção desafios das engenharias: Engenharia mecânica 2 Capítulo 3 HER 


(rã) 


1 > 


J (E) 


For this application, it has been considered , then: 


flEJF CE) 


a = Eb TESP — FED E) (6) 


2.3 Regula-falsi method 


The Regula-Falsi method is a simple iterative technique, which consists in considering 
two initial approximations X1 and X2, such that f(X1) and (X2) have opposite signs, i.e.: 


f(X1)- (X2) < 0 (7) 
so that can be determined by considering the equation of secant line as the function in the 


interval, as follow: 


pe x f(x2) — xo f(x) 
2 f(x) — f(x) 


Aa — Xz 


(8) 




















<E (9) 


Ka ka 
for a given €, then X3 is the root searched. Otherwise, f(:) is calculated from the choice of a 
value X1, between X1 and X2, so that f (X) has the opposite sign of f (X3). From this point, 
X3 is calculated, then X4, and so on. The process should be repeated until a root with a 
specific accuracy is obtained. The formulation for the iterative process can be summarized 


as. 


— Es (Ei — esin(E;)) — Ei(Ei-, — esin(Ei-1)) 


E = 
a E; — Ei, + e(sin(E;- 4) — sin(E;)) (10) 


2.4 Successive approximations method 


In this iterative method, a function , continuous in the range where the root must be 
found, is chosen. Let be rewritten as (Franco, 2006): 


Co Ro EO o an Um) 


Considering: 
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q0)=x+A (9) *F (o) (12) 


and that when x is the root of f (x), i.e., f (x) = 0, follows that x = q (x), for all A (x) z 0. 
Considering q (x) as defined in Eq. (12), if x is the root of f (x), then: 


q (x) = x (13) 


lt means that, on the point where X is the root of f (x), replacing the value of X in the 
function Q (X) will return the very X value. 
Therefore, this method consists in finding the numerical value that, when placed in Q (x), 


returns the X value. The iterative function is: 


Rag O (x) (14) 


in which n is the actual nº” iteration. 


2.5 Lagrange inversion theorem 


The Lagrange Inversion Theorem is an exact analytical method that not relies on 
numerical manipulations. This method provides an analytical solution for non-linear equations 
in terms of an infinite series (Rathie et al., 2013). Consider the functional equation, for which 


the Kepler's equation is a special case (Battin, 1999): 


v=x+a*o (y) (15) 


in which Ot is considered a small parameter (identified as the orbit eccentricity in the Kepler's 
equation). It follows that was a function of x, can be expanded as a Taylor series with OL =(): 


(60) =yw0)+a D+ CD 
E a “da 2100? 





(16) 


which can be turned into the following power series, given by the Lagrange Inversion 


Theorem: 
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on dra 
y=1+) a PO" (17) 


n=1 


2.6 Fourier-bessel series expansion 


Another approach to determine an approximate solution to Eg. (1) is a representation in 
a power series (no iterative approach). This expansion, called Lagrange Expansion, presents 
the eccentric anomaly in terms of a power series of the eccentricity . However, this series 
does not converge for all values of the eccentricity. The Fourier-Bessel Series Expansion 
solves this problem, because it is convergent for all eccentricity values. The expansion is 
defined below and details can be found in Battin (1999) and Colwell (1992): 


E = n+25 1 1 f cos(kE — ke sin(E) de sin(kM) 18 
+ AR E s(kE : : | ( ) 
in which X is the number of terms required for the expansion to provide the value of E with 


the desired accuracy. 


31 RESULTS 


In this paper, four iterative methods were investigated regarding the calculation of the 
eccentric anomaly E. The reference solution was obtained via Newton-Raphson's Method, 
and its results were compared with three methods: Halley, Regula-Falsi, and Successive 
Approximations. This study also aims at evaluating two series-based methods: Lagrange 
Expansion Theorem and Fourier-Bessel Series Expansion. Since their approaches are 
different from the iterative methods, these results were compared directly to each other. 

There is no uniform way in which the various authors evaluate the efficiency of 
methods to solve Kepler's equation (Nijenhuis, 1991). Herein, the evaluation was carried 
out using a convergence criterion: (i) the number of iterations required until the residual 
error (in 'E) is lower than 102 rad for iterative methods and lower than 104 rad for series- 
based methods; and (ii) the processing time to achieve such tolerance. For the iterative 
methods, the choice of a residual error lower than 10º rad is due to the guarantee that 
the error in the orbit will be lower than the order of centimeters; while 10 rad is due to the 
series-based methods limitations. 

In order to assess whether the behavior of these variables differs according to 
eccentricity and mean anomaly values, for each method, the calculation of the eccentric 
anomaly was implemented as follows: 


* | 21 values were taken into account for the eccentricity, over the entire elliptical 
interval ; 
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* | 11 values were considered for the mean anomaly, over the interval [O,r1. 


Therefore, Kepler's equation was solved for a grid of points in the (e, M)-plane. The 
increments for the grid size were A e = — and AM = — rad, which leads to 231 pairs of points 
(e, 'M). The calculations covered the value € = 1, which is orbitally, but not mathematically 
pointless. 

It is important to explain that the results could have been discussed around the error 
(the absolute difference between each method's solution and Newton-Raphson's solution). 
However, as the admitted errors are low, especially for the iterative methods, such analysis 
was not significant in this application. 

Regarding the processing time, Tab. (1) carries the information from the computer 


used to perform the simulations. 


Intel É Core TM |7-7500U CPU O 2.70 GHz 2.90 GHz 
8.00 GB 

4.00 GB 

64-bit operating system, x64-based processor 
MATLAB & 


Processor 
Installed memory (RAM) 
Cash memory 


System type 





Programming language 
Table 1. Specifications of the computer used for the simulations. 


3.1 Number of Iterations Assessment 


Iterative Methods 


The first analysis will be around the number of iterations that each method needed 
to reach the specified accuracy (according to the method's nature: iterative or in series), for 
each pair of points used as the initial guess. In all graphs, the grid of 231 pairs of previously 
established points and its corresponding performance are mapped. 

For better analysis and understanding of the graphs, two elements were used to rank 
the resulis: 


* | Size of the circles that represent each pair results: directly proportional to the 
number of iterations, that is, the smaller its size, the smaller the number of ite- 
rations, and vice versa. 


* | Color of the circles that represent each pair results: the color gradient goes from 
blue (fewer number of iterations) to yellow (greater number of iterations), and 
the green color portrays the intermediate values. 

Figure (1) shows the analysis related to the Halley's method. For all initial guesses, 
the Halley's method requires none, 3 or 4 iterations to achieve accuracy lower than 10! rad. 
The mean anomaly M represents the conversion to angle of the time elapsed since the body 
passed through the orbit perigee (M = 0). f M =0O or M =r,itisknownthat E =0or E =7, 
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respectively, and if e = O, automatically M = E; in these cases, the numerical solution becomes 
unnecessary (the method does not require iteration), as shown in Fig. (1). For all other regions, 
3 or 4 iterations were sufficient for Halley's method to achieve accuracy, including e = 1, which 
does not represent elliptical motion (parabolic orbit). The results for this method were very 
similar to those of the reference (Newton-Raphson method) and, for this reason, the reference 
was withdrawn from the graphs and its results will be presented in Tab. (2). 


Halley's Method: Number of lterations 
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a 

3 E] a = a a al had a a a * . a a s . 1] a » O) 
= 

qa n = z a] e a o] . . » z a z O] A] E] . n 

EI n da a O] o) + LL] ] = & . O) dd õ & E] . a a 

| a = od n m = » . a . » o a = z a » LL] = » 


Eccentricity (ECC) 
«paga 


Figure 1. Number of iterations required for Halley's Method to achieve accuracy lower than rad. 


Figure (2) exhibits the results regarding Regula-Falsi method. For the entire grid of 
initial guesses, 1 to 21 iterations were needed in order to reach accuracy lower than 102 rad. 
As with Halley's method, if M = 0, M =m or e = 0, the numerical solution becomes unnecessary 
and, according to the nature of the method, iteration was necessary, according to Fig. (2). The 
Regula-Falsi method presents greater sensitivity in relation to the eccentricity variation, that 
is, as e — 1, the number of iterations necessary to score accuracy increases, for any value of 
M in the interval [0, 71). 


Regula-Falsi Method: Number of Iterations 
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Figure 2. Number of iterations required for Regula-Falsi Method to achieve accuracy lower than 'º'2 rad. 
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Figure (3) presents the method of Successive Approximations performance. For each 
pair (e, 'M), it took from 1 to 1790 iterations so that the method reaches the required accuracy. 
EM =0, M=more=o0, the behavior is exactly the same as the Regula-Falsi method, 
as shown in Fig. (3). The number of iterations in the Successive Approximations method is 
even more sensitive to the variation of the eccentricity, and when 0.8 < 1 and — <M<rT, 
simultaneously, a problematic region starts to stand out, where convergence still occurs, but 
it is very costly from a computational point of view, requiring more than 100 interactions to 
achieve accuracy smaller than 10!2 rad. 


Regula-Falsi Method: Number of Iterations 


z 
5 LE ” à é * * = a é a & a é a Õ à à a a * + 
E | 
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Figure 3. Number of iterations required for Successive Approximations Method to achieve accuracy 
lower than 10º rad. 


From the results presented, regarding the number of iterations necessary to achieve 
an accuracy lower than 1072 rad, it is clear that, among the iterative methods, the Halley's 
Method is much more efficient than the other two methods, and the Successive Approximations 
methods, on average, is the one which requires the greater number of iterations to reach 
the established accuracy. While a maximum of 21 interactions were required in the Regula- 
Falsi method, the Successive Approximations method performed up to 1790, in order to 
achieve convergence. From the graphs, there is a difference in the critical regions in which 
the methods required a greater number of iterations: the Successive Approximations method 
needs a greater number of iterations when eccentricity is greater than 0.90 and mean anomaly 
is greater than 1.88, while the Regula-Falsi method reaches its maximum for eccentricity 
between 0.5 and, and mean anomaly between 0.63 and 1.57. 

In 34 combinations of eccentricity and mean anomaly values (15% of the possible 
combinations), both methods obtained the same number of iterations, in 78 combinations 
(34%), the Regula-Falsi method required a greater number of iterations compared to the 
Successive Approximations method. However, the largest difference in these cases was 15 
iterations. Finally, in 119 combinations (52%), the Successive Approximation method had a 
greater number of iterations than the Regula-Falsi, with a difference of up to 1783 iterations 
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between them (of these, in 7 cases, the difference was greater than 50 iterations). 


Series-Based Methods 


Figures (4) and (5) discuss the methods based in series performance: Fourier-Bessel 
Series Expansion and Lagrange Inversion Theorem, respectively. Due to the distinct nature 
of these methods, they were compared with each other, separately from iterative methods. 
The difference of performance between the two is huge: for each pair (e, M), while the 
Fourier-Bessel Series Expansion needs a maximum of 13 iterations to achieve an accuracy 
less than 107, the Lagrange Inversion Theorem needs up to 4900 iterations (- 376 times 
more than Fourier-Bessel Series Expansion). In all 231 combinations of eccentricity and 
mean anomaly, the Lagrange Inversion Theorem required a greater number of iterations 
than the Fourier-Bessel Series Expansion. For entries, M=0, M=— and M=m(Ve €E[0,1]) 
ore=0(Y M€ [0, ri), both methods converge quickly. The Lagrange Inversion Theorem 
high number of iterations elsewhere in the grid indicates that a more accurate analysis 
needs to be made on this expansion. 


Fourier-Bessel Series Expansion: Number of Ieratlons 
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Figure 4. Number of iterations required for Fourier-Bessel Series Expansion to achieve accuracy lower 
than 10% rad. 
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Lagrange Inversion Theorem: Number of Iterations 
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Figure 5. Number of iterations required for Lagrange Inversion Theorem to achieve accuracy lower than 
10%rad. 


In Tab. (2) some descriptive measures of the number of iterations are presented, both 
for iterative and series-based methods, even though the series-based methods have been 
compared separately. As the average is a measure of position sensitive to outliers, it was 
decided to analyze its results together with the median; the standard deviation is the measure 
of data dispersion around its average. In the first line of Tab. (2) are the statistics obtained for 
the Newton-Raphson's method (reference) and the similarity of its results with those obtained 
via Halley's method stands out, as already discussed. 

So far, ffom the results presented and Tab. (2) analysis, it is possible to list the 
methods, identifying which ones are the most efficient in relation to the number of iterations. 
The iterative Halley's method showed the best results, followed by the Fourier-Bessel series 
expansion and, later, the iterative Regula-Falsi method. These three methods presented close 
mean and median values and a low standard deviation value, which indicates that there are 
no outliers, that is, for all the values of eccentricity and mean anomaly used in this work, there 
was not a large number of iterations necessary for convergence. The methods that presented 
less satisfactory results were the iterative Successive Approximations methods and, finally, 
the Lagrange Inversion Theorem. Both presented mean values much greater than their 
respective medians, in addition to a high standard deviation value. These results indicate the 
presence of outliers in the data, that is, for certain values of eccentricity and average anomaly, 
these methods required an anomalous number of iterations, quite different ffom the usual one. 
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Standard 


Number of Iterations Mean deviation Median 
Newton-Raphson 3.446 1.921 

Halley 2.931 1.608 

Regula-Falsi 9.883 6.025 

Successive Approx. 33.468 125.631 13 
Fourier-Bessel 5:52 2.926 

Lagrange Theorem 617.420 1276.7083 81 


Table 2. Statistics of number of iterations analysis. 


In order to improve the analysis and to measure the impact of the findings related to 
the number of iterations in the processing time, additional results are presented in the next 
subsection 


3.2 Processing Time Assessment 


Iterative Methods 


The processing time was measured posteriori and, therefore, involves the actual 
execution time of the algorithams. It depends on factors related to the machine, the 
programming language used and, sometimes, it is a function of additional aspecis of a 
particular input (Linder, 2021). 

Figure (6) shows the measured processing time for the reference method, Newton- 
Raphson, over the entire grid of 231 possible inputs. It is possible to perceive that one of 
the limit regions, in which the orbits are circular or quasi circular (e —> 0), has the longest 
processing time, for any value of 'M. Thus, as the algorithm is the same for all initial guesses, 
the method is sensitive to this entry of eccentricity. For the other points, the processing time 
starts to increase if — <M<7m, Ve E (0,1]. All the other iterative methods studied here 


showed shorter processing time in this specific situation. 


Newton-Raphson's Method: Processing Time(s) 
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Figure 6. Processing time for Newton-Raphson's Method (iterative reference). 
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Figure (7) presents the processing time required by the method of Successive 
Approximations, taken into account the 231 pairs of. It is noticeable that the processing time 
in the limit region where is higher, a behavior that is also observed in Fig. (6). As for the other 
points, the longest processing times occurred , when: 


*  1.88<M<3.14 for the Successive Approximations Method. 
* 251<M<3.14 for Regula-Falsi Method. 
* 1.88<M<3.14 for Halley's Method. 


Based on these results, it is concluded that particular aspects of the input values need 
to be better analyzed, which is beyond the objective of this study. 

As a behavior very similar to that shown in Fig. (7) occurred for the methods of Halley 
and Regula-Falsi, the respective graphs will be omitted and more results related to these 
methods will be discussed in the statistical analysis of Tab. (3). 


Successivo Approximations Method: Processing Time(s) 
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Figure 7. Processing time for Successive Approximations Method. 


Series-Based Methods 


Here, it was used the same condition applied in the number of iterations assessment: 
the methods whose solution is based on the series approach were analyzed separately. 
Figures (8) and (9) show the measured processing time of the Fourier-Bessel Series 
Expansion and the Lagrange Inversion Theorem, respectively. These methods called for a 
considerably longer processing time, which was expected, due to the nature of the approach 
and the number of operations that needed to be carried out to obtain the methods resulis. 
The minimum processing times are very close to the equivalent times required by the iterative 
methods, and the differences emerged in the maximum processing times. When analyzing 
Fig. (8), it is noted that the largest order of Fourier-Bessel Series Expansion processing time 
is similar to the reference ones, although its occurrence is greater, that is, most of the poinis 
in the grid demanded the maximum time to reach convergence. In addition, Fig. (9) shows 
that the maximum processing time obtained for the Lagrange Inversion Theorem is up to 
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times greater than the Newton-Raphson's solution equivalent time, which designates this as 


the least efficient method among all. Furthermore, the same problematic region, previously 


observed, starts to stand out ( and, simultaneously), in which convergence is detected, at the 


cost of a higher computational burden. 


Nemwton-Raphson's Method: Processing Time(s) 
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Figure 8. Processing time for Fourier-Bessel Series Expansion. 


Lagrange Inversion Theorem: Processing Time(s) 
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Figure 9. Processing time for Lagrange Inversion Theorem. 


Table (3) contains data from the statistical analysis of the six methods, which again 


include mean, standard deviation and median. The first line of Tab. (3) presents the statistics 


obtained for the Newton-Raphson method (reference) anditis clear that, in terms of processing 


time, the iterative methods were competitive with each other, due to the behaviors shown in 


the graphs, confirmed through Tab. (3) statistics. Series-based methods were computationally 


more costly, which is an expected conclusion, given their approaches. Even so, Fourier- 


Bessel Series Expansion sustained a competitive performance, compared to the iterative 


methods, while Lagrange Inversion Theorem results did not, which statistically corroborated 


the discrepancies detected in the graphs. In theory, if more terms are added to the Lagrange 


Inversion Theorem series, an increase in processing time is expected (although this analysis 
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needs to be done with more criteria). 


Processing Time (s) Mean Standard deviation Median 

Newton-Raphson 1.652107 5.77710* 3.78010% 
Halley 1.3481074 4.261104 3.990105 
Regula-Falsi 1.5511074 5.095101 4.30010% 
Successive Approx. 1.4751074 4.2731074 5.12010º% 
Fourier-Bessel 4.33910% 1.0481072 9.620101 
Lagrange Theorem 1.327 3.481 3.940107? 





Table 3. Statistics of processing time analysis. 


41 CONCLUSIONS 


The goal of this paper was to review an application of six methods, iterative and 
in series, for the solution of the Kepler's equation, over a grid of points in the elliptical 
motion interval. The analysis is based on the number of iterations that each method requires 
to achieve a stipulated accuracy and on the post-processing time (execution time). The 
comparison was separated according to the nature of the methods. The solution obtained 
using the Newton-Raphson's method worked as reference for the iterative methods. 

Regarding the number of iterations, Halley's method performed as competitively as 
the reference solution. Regula-Falsi and Successive Approximations methods, on the other 
hand, were more sensitive to the variation of the eccentricity input value, in the interval (0,1]. 
For the latter method, even, a problematic region began to be laid out: «e E (0.81) (quasi- 
parabolic orbits) and M € — sTT |, simultaneously, where convergence still occurs, but it 
is very computationally costly, in terms of processing time and/or number of iterations. It 
implies that the method of Successive Approximations is not stable for very eccentric orbits 
(e — 1), considering that e = 1 does not configure orbital elliptic motion. 

The processing time analysis indicated competitiveness among the four iterative 
methods and the Fourier-Bessel Series Expansion, observed both in the graphical and in 
the descriptive measures” analysis. The only method that presented results that were very 
distant from the others, with maximum processing time of the order of 100 times greater than 
the others, was the Lagrange Inversion Theorem. Consistently, it had already exhibited signs 
of convergence trouble in the number of iterations analysis for the same region, because the 
same region stood out in the analysis of the iterative Successive Approximations methods. 

The difference presented by the Lagrange Inversion Theorem results for both 
number of iterations to reach the required accuracy and processing time indicates that a 
more accurate analysis needs to be made on this expansion and the impact of the number 
of terms in the series and the truncation order should be studied further. 
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RESUMO: Nos países desenvolvidos a demanda 
de energia per capita tende a ser muito superior 
daqueles que possuem um modelo econômico 
ineficiente. Também é fato que a energia vem 
sendo consumida em maior quantidade com o 
passar do tempo. Encaixando-se sob a forma 
de armazenamento de energia potencial, o 
armazenamento de energia por compressão 
de ar de forma adiabática — A-CAES (do inglês, 
Adiabatic Compressor Air Energy Storage) pode 
funcionar como regulador de sistemas de energia 
elétrica gerada por sistemas renováveis ou 
instáveis e torná-la disponível no local e na hora 
adequada de utilização, sem precisar enviá-la para 
a rede de distribuição, que pode estar indisponível 
ou inapta para recebê-la. Existem algumas plantas 
industriais feitas no mundo, sendo a de Huntorf 
na Alemanha (início em 1978) uma das principais 
elencadas pelo seu pioneirismo. Muitas utilizam 
o conceito de compressão estagiada por turbo- 
máquinas de grande porte e geração com o 
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mesmo equipamento. As utilizações são um pouco 
diferentes, porém visam garantir o fornecimento 
contínuo e estável de energia. Este trabalho visa 
a apresentar um protótipo laboratorial concebido 
e construído na Universidade Tecnológica 
Federal do Paraná (UTFPR), Campus Ponta 
Grossa, para estudo da eficiência do A-CAES em 
microssistemas, com extensão de resultados para 
concepções de diversos portes. 
PALAVRAS-CHAVE: A-CAES, 
Armazenamento,  Huntorf, Ar 
Cavernas, Renovável. 


Energia, 
Comprimido, 


A-CAES (ADIABATIC COMPRESSOR AIR 
ENERGY STORAGE): EXPERIMENTAL 
EQUIPMENT ON LABORATORY SCALE 


ABSTRACT: In developed countries, the demand 
for energy per capita tends to be much higher 
those with an inefficient economic model. It is also 
a fact that energy has been consumed in greater 
quantity over time. Fitting in the form of potential 
energy storage, adiabatic air compression energy 
storage (A-CAES) can function as a regulator for 
electrical power systems generated by renewable 
or unfit systems and make it available on-site and 
at the appropriate time of use, without having to 
send it to the distribution network, which may be 
unavailable or unfit to receive it. There are some 
industrial plants made in the world, the Huntorf plant 
in Germany (beginning in 1978) one of the main 
ones listed for its pioneering. Many use the concept 
of compression aged by large turbo-machines and 
generation with the same equipment. The uses 
are somewhat different, but aim to ensure the 
continuous and stable supply of energy. This work 
aims to present a laboratory prototype designed 
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and built in UTFPR/Ponta Grossa to study the efficiency of A-CAES in micro-systems, with 
extension of results for conceptions of various sizes. 
KEYWORDS: A-CAES, Energy, Storage, Huntorf, Compressed Air, Caves, Renewable. 


11 INTRODUÇÃO 


Há um crescente aumento de consumo e geração de energia a nível mundial, 
especialmente após a passagem do milênio. Além disso, é um desafio permanente ter 
energia para utilização no tempo exato da necessidade. Conforme Castelloes (2016), as 
formas mais comuns de armazenagem não são tão amistosas ecologicamente e todas 
possuem delimitações de uso no seu tamanho físico e aplicação — Figura 1. O sistema de 
armazenagem através da energia potencial hidráulica de um rio promove a inundação de 
àreas de mata, deslocamentos de pessoas, alterações de biomas, não obstante seja bastante 
utilizado, especialmente no Brasil onde as hidrelétricas usam esta funcionalidade para gerar 
e ter o sistema sempre disponível. Os sistemas de armazenagem através da energia química 
são muito densos em termo de energia, porém limitados em quantidade, tamanho e geração 
de resíduos sólidos potencialmente perigosos ao meio ambiente. E também, há direções de 
utilização de energias renováveis e de geração distribuída de energia, de forma mundial, 
assim, em pouco tempo, haverá diversos sistemas médios que necessitarão de reguladores 
de energia, mantendo-se o sistema estável. Dessa forma, será possível retornar energia 
ao sistema produtivo quando houver uma diminuição no rendimento, ou mesmo, reinserir 
toda a energia durante um intervalo onde não haja produção nenhuma. Exemplos disso 
podem ocorrer em falta de ventos ou à noite nos sistemas fotovoltaicos de energia solar 
(Ramadan et al., 2016). Com o aumento de produção de energia solar, uma cidade de médio 
porte poderá incentivar a produção através de sistemas renováveis fotovoltaicos ou eólicos, 
independentemente da rede de distribuição, que possam ser usados durante o dia ou horário 
e estabelecimento de mecanismos estabilizadores de produção, como armazenagem em 
grande reservatório de ar comprimido. Com isso, cumpre-se a evolução técnica eco ambiental 
de produção de energia de uma forma mais prática para a geração distribuída (Ibrahim et al., 
2015). 
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Espectro de Aplicações e Tecnologias de Armazenamento de Energia Elétrica 
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Figura 1. Espectro de aplicações e tecnologias de armazenamento de energia elétrica (Castelloes, 
2016) 


A energia elétrica armazenada sob a forma potencial de pressão de ar comprimido 
poderá ser transformada em energia elétrica em turbinas eólicas de alto rendimento, em 
uma planta térmica de compressão e de transmissão (Milewski et al., 2016). A formação 
do Regulador Hidráulico Operacional (RHO) pode ser dividida em diversos vasos de 
pressão ou em um grande vaso único de pressão subterrâneo, chamado de reservatório 
apenas. Com isso, poder-se-á estudar novos formatos mais compactos de turbinas 
aerodinâmicas impulsionadas por ar comprimido, com a vantagem de não emissão ainda 
de gases poluentes ou residuais de combustão (Hammann et al., 2017). Também o sistema 
poderá ser estendido a dois tipos de porte: médio e micro. Os sistemas de médio porte 
como mencionados poderão abastecer boa parte de uma cidade em regime noturno, por 
exemplo. Enquanto os sistemas de micro porte poderão atender locais onde não possa 
haver perdas ou interrupções de fornecimento de energia, como hospitais. Também 
poderão ser utilizados em horários de pico, diminuindo os custos com energia elétrica e 
também com combustíveis utilizados em geradores térmicos que tem a finalidade de cobrir 
estes horários de alto valor financeiro da energia elétrica (Diyoke et al., 2018). Há de se 
reconhecer que existem fatores limitantes, como por exemplo, que a produção deverá ser 
superior ao consumo, de tal forma que haja estocagem para uso no período adequado. 

Nesse contexto, estudar o Armazenamento de Energia em Ar Comprimido de Forma 
Adiabática — A-CAES (do inglês, Adiabatic Compressor Air Energy Storage) como forma 
para armazenamento de energia renovável proveniente de fontes solares para disponibilizá- 
la nos horários noturnos ou de baixa captação de radiação solar, de forma a estabelecer um 


regulador, é de extrema importância. 
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21 METODOLOGIA 


Este artigo está aderente à temática de Engenharia Mecânica e de Energia como 
linhas de pesquisa, utilizando às Ciências Térmicas (Mecânica dos Fluidos, Termodinâmica e 
Transferência de Calor), a Eletricidade, os motores e geradores elétricos concomitantemente 
as áreas de Energias Renováveis (Crawley, 2017). Tratam-se de grandes áreas aglutinadoras 
que resultam em um projeto de pesquisa com objetivos específicos para atuação na área de 
Armazenamento de Energia por Ar Comprimido de Forma Adiabática (A-CAES). Por meio de 
uma pesquisa bibliográfica de artigos científicos em importantes bases de dados científicas, 
foi realizada uma revisão bibliográfica para estabelecimento dos fundamentos científicos 
visando ao projeto do aparato experimental em escala laboratorial. Para formulação dos 
cálculos específicos do sistema A-CAES, a metodologia empregada consistiu em informações 
técnicas de livros próprios do tema (Rollins, 2004; Rosa, 2015; Crawley, 2017; Moreira, 
2017). Sistemas experimentais, como por exemplo, os utilizados em Budt et al. (2016), Wang 
et al. (2017) e Castellani et al. (2018) analisados com objetivo de estabelecer o nível mínimo 
de atendimento das expectativas com relação ao projeto executivo e instalação do aparato 
experimental em escala reduzida no Laboratório de Meios Porosos e Eficiência Energética 
(LabMPEE) vinculado ao Programa de Pós Graduação (Mestrado) em Engenharia Mecânica 
(PPGEM) da Universidade Tecnológica Federal do Paraná (UTFPR), Câmpus Ponta Grossa. 


31 DESCRIÇÃO GERAL DOS SISTEMAS CAES E A-CAES 


De acordo com Jafarizadeh et al. (2020) é improvável fazer um estudo relacionado ao 
CAES sem elencar a Planta Térmica de Huntorf, localizada na Alemanha. Concebida para 
regularizar a rede de energia elétrica proveniente de usinas eólicas localizadas no Norte 
daquele país, foi a primeira a entrar em operação comercial, no ano de 1978. Trata-se de um 
sistema onde se armazena energia pneumática potencial em cavernas salinas subterrâneas, 
com volume aproximado de 310.000mº. Essa planta tem 290MW de potência e opera segundo 
um ciclo diário (Kaldemeyer et al., 2016). As principais características operacionais dessa 
planta estão apresentadas na Tabela 1. Esse sistema tem sido aperfeiçoado desde a sua 
inauguração e os resultados simulados relatados na literatura dão conta de um rendimento 
global de 42% (Luo et al., 2015). 

Os novos sistemas propostos encontrados na literatura, fazem a consideração da 
possibilidade de recuperação de calor gerado durante o ciclo de compressão e utilização no 
processo de expansão dasturbinas. Essa energiatérmica gerada pelo processo de compressão 
é captada por um sistema simples composto por trocadores de calor (Heat Exchangers ou 
HEXs), armazenada em uma planta térmica (TES — Thermal Energy Storage) composta por 
um reservatório quente e outro frio, basicamente (Chen et al., 2017). O ar é armazenado a 
baixa temperatura e expande a alta temperatura, recebendo calor dos trocadores de calor. 
Normalmente o fluido de trabalho selecionado do TES é a água. A composição do sistema 
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CAES com o TES é chamado de A-CAES (do inglês, Adiabatic Compressor Air Energy 
Storage ou, em português, Armazenamento de Energia por Compressão de Ar de Forma 
Adiabática). Isto quer dizer que somente o fluido de trabalho do TES e o Ar estão envolvidos 
no sistema (Luo et al., 2015). Não há queima de combustível (Gás Natural ou GLP ou outro) 


envolvido no sistema mostrado topologicamente na Figura 2. 
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Item Denominação Valor Unidade 
Compressores: linhas de baixa e alta pressão (conexão em série) 
01 Vazão mássica de ar dos Compressores 108 kg/s 
02 Potência de todos compressores 60 MW 
03 Temperatura de saída do ar do after cooler 50 ºC 
04 Máxima de Ar de saída do after cooler 72 kgf/em2 
Cavernas: dois reservatórios subterrâneos 
05 Volume total (físico) das cavernas (aproximado) 310.000 m3 
06 Faixa de temperatura das paredes das cavernas 40 a 50 ºC 
07 Pressão operacional dos reservatórios 43 a 66 kgf/em2 
08 Pressão máxima admissível 70 kgf/cm2 
09 Queda de pressão operacional 10a 15 kgf/em2 
Turbinas: linhas de alta e baixa pressão (conexão em série) 

10 Vazão mássica de gás nas turbinas 417 kg/s 
11 Potência de todas as turbinas 290 MW 
12 Pressão de entrada na turbina de alta pressão 43 kgf/em2 
13 Pressão de entrada na turbina de baixa pressão 8 a 11 kgf/em2 
14 Razão temporal de carga/descarga de todo sistema 4 - 


Tabela 1. Especificações básicas da Planta Huntorf/Alemanha (Luo et al., 2015). 


O funcionamento básico desse sistema considerando as especificações da Planta 
HuntforlALE é: a rede elétrica instável é conectada ao sistema CAES e gerida pelo Operador 
Elétrico, que controla a entrada e saída de energia para o chamado grid on, ou seja, o operador 
é que verifica e gera as demandas necessárias ao CAES. No Sistema de Compressão o Ar 
é admitido a pressão atmosférica pelo Compressor 01 (BP — Baixa Pressão) e é impelido 
para outros compressores em série denominados boosters (AP — Alta Pressão) que elevam 
a pressão do ar a aproximadamente 72kgf/cm?, conduzindo o ar até o Reservatório de Ar 
Comprimido. Durante os estágios de compressão a entalpia, a pressão e a temperatura do ar 
vão se elevando. Em cada estágio há um trocador de calor que refrigera o ar, recebendo calor 
do mesmo e transmitindo ao fluido de trabalho (geralmente água). A temperatura máxima 
de saída no último estágio de compressão, do ar, é de aproximadamente 163ºC e a mínima 
após o último estágio de expansão (turbinas) é de aproximadamente 101ºC. A temperatura 
da água nos reservatórios quente e frio são 109ºC e 69ºC, respectivamente. 
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Figura 2. Desenho esquemático de um sistema A-CAES (adaptado de Luo et al. (2015). 


A água armazenada no reservatório quente provém dos trocadores de calor e after 
coolers instalados nos estágios de compressão. Esse reservatório fica isolado visando evitar 
perdas térmicas. Esse fluido quente é bombeado para as turbinas onde há trocadores de calor 
em cada estágio de expansão do ar comprimido, com a finalidade de aumentar a temperatura 
do ar, entregando calor ao mesmo e retornando ao reservatório frio. A variação média de 
temperatura entre frio e quente é de 40ºC. O reservatório frio não é isolado, bem como 
o Reservatório de Ar Comprimido. A finalidade do TES (sistema composto por trocadores, 
reservatórios frio e quente e sistema de bombeamento) é evitar a perda de calor através do 
Reservatório de Ar Comprimido. De forma lógica, o Reservatório de Ar Comprimido poderia 
ficar a alta temperatura se não houvesse o sistema de trocas de calor, porém, devido às 
dimensões e irregularidades nas paredes do mesmo é praticamente inviável acontecer o 
isolamento térmico do mesmo (cavernas de grandes dimensões). Dessa forma, esse calor 
que seria perdido é “recuperado” e aproveitado no aumento da entalpia na saída do sistema 
de expansão (turbinas). A temperatura das cavernas é em torno de 50ºC. A Planta Huntorf 
tem funcionado segundo um ciclo de 8 horas de compressão e 2 horas de descarregamento 
(produção elétrica), pois uma avaliação mostrou que esse ciclo diário apresentaria os 


melhores resultados. 


41 BASE DE CÁLCULO DO SISTEMA A-CAES 


O local de instalação do aparato experimental é a cidade Ponta Grossa no estado do 
Paraná no Brasil com Latitude Sul de 25º4'8”, Longitude Oeste de 50º8'7” e com altitude de 
954m do nível do mar. O projeto de armazenamento de energia através de ar comprimido 
de forma adiabática foi baseado na seleção de equipamentos de uso comercial, procurando- 
se o máximo rendimento global para uma bancada experimental de escala reduzida a ser 
instalada na edificação do Laboratório de Meios Porosos e Eficiência Energética (LabMPEE) 
da UTFPR/Ponta Grossa. A Figura 3 mostra uma fotografia de parte da instalação do sistema 
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construído — Casa de Compressão e Geração. Os equipamentos básicos disponíveis para a 


montagem experimental foram: 


* | 01 compressor Chiaperini'” CJ15 duplo estágio, 176 psi, 424 I/min, 2,2 kW / 
1100 RPM; 


* | 01 motor de palhetas (turbina) Kazi'” — KZA750 7,5 HP, 80 I/s, 17 N.m; 


* | 01 alternador elétrico Toyama!” 3,5 kVA / 6 pontas / trifásico / cos q = 0,80, m = 
63%), 1800 rpm; 


* | 01 reservatório de ar comprimido 12 bar, 200 litros, horizontal acoplado ao com- 
pressor (APV200); 


* | 01 Conjunto de painéis fotovoltaicos 5,6 kW-p / com inversor 220 V/ trifásico. 


5 | PROJETO AS BUILT 


Executados os cálculos teóricos fundamentais para início do projeto de um sistema 
experimental A-CAES em escala laboratorial, os esforços foram direcionados para cálculos 
finais seletivos de equipamentos comerciais para instalação de um sistema compatível. 


o. c 
ves = 


“INDUSTRI 


CMIADERr 1 





Figura 3. Casa de Compressão e Geração: compressor, reservatório, turbina, alternador 
elétrico sistema pneumático de controle, válvula elétrica, válvula, registro, regulador e lubrificador. 


b.1 Seleção dos Painéis Fotovoltaicos 


Para viabilizar de forma cabal a captação através de uma fonte de energia renovável 
foi executada a seleção dos painéis fotovoltaicos. Para a captação fotovoltaica foram 
montados dois grids, sendo 10 painéis fotovoltaicos Yingli Energy'” de 150W e 24 painéis 
fotovoltaicos Upsolar'" de 170W, perfazendo 5,58 kWp, podendo atender a carga solicitada 
pelo compressor. Esses 34 painéis fotovoltaicos foram instalados no telhado da edificação do 
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LabMPEE/UTFPR (Figura 4). A seleção prévia foi realizada pelo software shareware PVSYST 
V6.84 (Pereira & Oliveira, 2015). Esses cálculos estão nas Tabelas 2, 3,4 e 5. 


Stand alone system: Parâmetros de Simulação 
Sistema Mini A-CAES Laboratorial - LabMPEE/UTFPR/Ponta Grossa 
Modulo FV Si-poly Poly 110 Wp 72 células 


Utilização: 5 dias/semana | Quantidade: 01 | Potência: 2200 W | Utilização: Gh/dia | 13200 Wh/dia 
Compressor pneumático 220 V Standby consumer 24 Wh/dia 


Tabela 2. Simulação fotovoltaica para seleção dos painéis fotovoltaicos (exigência do consumidor). 





Sistema Mini A-CAES Laboratorial UTFPR/Ponta Grossa 
Latitude -25.09º S Longitude -50.16 W 


Tempo definido como Tempo legal Fuso Horário Altitude 954 m 
TU-03 
Variante da Simulação SIMO4 Data da 15/10/2019 21h34min 
Simulação 


Transposição Perez Difuso Perez, Meteonorm 


Consumo doméstico Constantes ao longo do ano 


Características do Grupo Fotovoltaico 


Número de módulos FV Em paralelo 23 strings 


Número total de módulos Número de Potência 110 Wp 
módulos Unitária 
4501 Wp (5070) 


Tabela 3. Simulação fotovoltaica para seleção dos painéis fotovoltaicos (características). 





Sistema Mini A-CAES Laboratorial - LabMPEE/UTFPR/Ponta Grossa 
Energia disponível | 6181 kW-h/ano | Produção Específica 


ea ia utilizada | 3451 kW-hlano | — Excedente | 
sistema Energia utilizada 3451 kW-h/ano Excedente 
42,87% Fração Solar 


Tabela 4. Simulação fotovoltaica para seleção dos painéis fotovoltaicos (resultados principais da 
simulação). 
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Técnica (China) (EUA) 


| âeaonemadomóaiom) [100 | om 


Tabela 5. Características dos painéis fotovoltaicos selecionados 








Figura 4. Conjunto de painéis fotovoltaicos do sistema A-CAES 


5.2 Sistema Elétrico 


O projeto de painel elétrico com gabinete e componentes visa o controle, medição 
e armazenamento de dados de sistema de cogeração elétrica e do sistema fotovoltaico, 
procurando-se à aquisição de dados da bancada experimental instalada. Como mencionado 
anteriormente, a potência geral instalada foi de 5,58 kW elétricos. Esse sistema elétrico é 
composto basicamente por: 


a. Um gabinete em aço, uso interno IP21, com porta com fecho triangular, com 
assentamento em base de alvenaria, placa de fundo, barramento em cobre 
eletrolítico e sistema de identificação de circuitos e componentes; 


b. (Componentes assentados em trilhos DIN, fiações compatíveis, anilhas identi- 
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ficadoras em cada cabo, cabos coloridos segundo o padrão ABNT, cabos de 
força, disjuntor para compressor (2,2 kW / 220 V/ trifásico) e para sistema de 
cogeração elétrica (3,5 kVA / 220/127V) e disjuntor geral de seccionamento (6,0 
kW) (Figura 5); 


c. CLP para controle lógico do sistema e Inversor Solar 6,0 kWPp — trifásico (220V) 
integrados (Figura 5); 


d. Botoeiras liga/desliga, lâmpadas de sinalização; 
e. Contactoras e sistema de controle da válvula elétrica proporcional 1”; 


f. Sistema QDCC com String e Combiner Box (forma de montagem) (Figura 6). 


Essa contabilização de energia produzida pela carga demandada pelo motor do 
compressor ocorre no inversor solar. Parando-se o motor, um sinal é enviado pelas portas 
de comunicação com o CLP, computando e armazenando a energia produzida neste período 
de funcionamento. A contabilização da energia pneumática (potencial) ocorre através da 
curva de pressão, cujo sinal é gerenciado por um sensor de pressão IFM'M e pelo tempo 
de abertura da válvula elétrica. A válvula abre com 9,0bar de pressão na linha e fecha com 
4,0bar de pressão na linha. As pressões foram estabelecidas conforme o projeto pneumático 
da turbina, de tal forma que se tenha a rotação de 1800 RPM na mesma, vencendo-se o 
torque necessário para girar o alternador elétrico. Com a geração elétrica pelo alternador, um 
conjunto de resistências elétricas (3x 700W) dissipam a produção elétrica — Figura 8. 
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Figura 5. Quadro de comando com gabinete 





Figura 7. Inversor WEG!” e CLP WEG!M 


5.3 Sistema de Compressão e Cogeração Elétrica 
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Figura 6. Croqui elétrico de geração 


fotovoltaica 





Figura 8. Banco de resistências elétricas 


O sistema de tratamento/fluxo/compressão de ar atmosférico e o sistema de 


cogeração elétrica foram instalados em um abrigo de alvenaria anexo às instalações internas 
do LabMPEE/UTFPR/Ponta Grossa conforme apresentado anteriormente na Figura 3. A 
interligação entre os sistemas é feita por tubulação executada em tubos de aço conforme 
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DIN2448. Como mencionado anteriormente, o sistema é composto por compressor recíproco, 
com refrigeração dos cabeçotes de admissão e descarga, pistões e sistemas mecânicos 
completos com válvulas. Esse é acoplado diretamente ao motor elétrico trifásico, 220V, de 
SHP de potência máxima. Reservatório 200 litros, pressão máxima admissível 12,0bar. O 
sistema de cogeração é composto por sistema abaixador de pressão por fluxo pneumático, 
que consumirá 80/s de a 4,0bar, produzindo 17N.m de torque e 7,5CV de potência no eixo, 
de tal forma que o mesmo gire a 1800RPM sincronamente para acoplamento a sistema 
alternador monofásico 220/127V, seis pontas, com acoplamento direto. Ambos equipamentos 
foram assentados em base metálica construída para nivelamento e alinhamento dos mesmos 
conjuntos girantes e interligados pneumaticamente e eletricamente com os demais sistemas 
(Figura 3). O acoplamento do sistema de turbina e gerador ocorre através de um elemento 
elástico, acoplamento tipo flexível. A turbina recebe ar proveniente do reservatório do 
compressor, que está a 12bar de pressão, este ar passa pelo lubrificador e regulador de 
pressão e é liberado a 9,0bar sobre a válvula elétrica. Com a pressão e o fluxo, esta energia 
potencial atinge as palhetas da turbina, fazendo-a girar acima do torque inercial do alternador. 
A ação da turbina sobre o gerador excita-o, produzindo energia elétrica. Quando a pressão, 
lida pelo sensor de pressão, atinge o mínimo de 4,0bar, a válvula elétrica é fechada pelo CLP. 


61 RESULTADOS PRELIMINARES & DISCUSSÃO 


Os testes experimentais preliminares foram executados no mês de outubro de 
2020, onde se encontrava na maioria das vezes dias ensolarados com incidência solar 
acentuada. As condições ambientais eram de pouco vento, pela manhã entre 10h e 12h. A 
coleta de informações ocorreu de forma automática, conforme previsto e foi disponibilizada 
em planilhas, sendo compilados os resultados preliminares nas Tabelas 6 e 7. A abertura 
da válvula elétrica foi de 65% do curso permitido. A rotação síncrona foi estabelecida sem 
problemas, apenas fazendo-se o setup no CLP. O resultado experimental para a eficiência 
global do sistema 


Mem | Energia [kWh] Eficiência Energética 
Motor Elétrico 0,1174 


Gerador Elétrico 0,0272 


Tabela 6. Eficiência energética experimental por fase do processo A-CAES 
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Potência | Energia 


Média Painéis | Corrente Potência Energia Tempo Energia Energia 


Compressor | Compressor | Geração | Gerada | Potencial 
o bm) ” [kW] [KW-h] [s] | [kWh] | [kWh] 


DT Esiênciaódiatom 


Eficiência Média Toa 





Tabela 7. Resultados preliminares experimentais de produção energética do sistema A-CAES 


A-CAES foi de 18% enquanto que o valor teórico da eficiência foi calculada previamente 
em 16% (Elmegaard & Markussen, 2011). Outra discussão sobre esses resultados 
preliminares pode ser feita em relação à eficiência de compressão de apenas 32%. Isso 
se deve principalmente na escolha comercial de um compressor recíproco mais simples e 
na forma de armazenamento (reservatório acoplado ao compressor de apenas 200 litros). 
Ressalta-se ainda, que os equipamentos usados na Alemanha e nos EUA são de portes 
diferentes e construções diferentes, tipo turbo gerador, onde é obtida uma eficiência maior e 
pressões cerca de 6 vezes maiores que a utilizada nessa bancada experimental (Budt et al., 
2016). Além disso, os sistemas são integrados com ciclos térmicos com uso de combustível 
fóssil. Em comparação com os volumes armazenados de 310.000mº às pressões da ordem 
72bar, muitas variáveis presentes nas irreversibilidades se encontram com maior nível de 
importância no aparato experimental. Mesmo assim, o sistema experimental apresentado se 


mostrou controlável com um sistema simples de CLP e poucos sensores. 


7 | CONCLUSÕES 


No presente trabalho foi apresentado um aparato experimental em escala laboratorial 
concebido e construído na Universidade Tecnológica Federal do Paraná (UTFPR), Câmpus 
Ponta Grossa, para estudo experimental da eficiência do Armazenamento de Energia em Ar 
Comprimido na Forma Adiabática — A-CAES (do inglês, Adiabatic Compressor Air Energy 
Storage) como forma de armazenamento de energia renovável proveniente de fontes solares 
para disponibilizá-la nos horários noturnos ou de baixa captação de radiação solar, de forma 
a estabelecer um regulador. A validação prática experimental do mini sistema A-CAES 
aconteceu de forma comprobatória neste experimento, usando-se de energia renovável, 
sem uso de combustíveis fósseis, sem geração de resíduos sólidos, cumprindo com seus 
objetivos principais de estabilização da geração da energia elétrica e de armazenamento 
ecologicamente orientado de energia potencial. 
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LA IMPORTANCIA DE LA FACTORIZACIÓN EN 
LA INGENIERIA MECANICA PARA DETERMINAR 
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RESUMEN: El presente trabajo consiste en 
el análisis cuantitativo a partir del teorema 
de Bernoulli involucrando a las variables que 
intervienen en los impulsores centrifugos que son 
utilizados en las turbos bombas para determinar 
en base al modelo matemático obtenido, el 
comportamiento mecânico energético de los 
impulsores pudiendo así trazar la curva de 
operación Ilamada también curva característica, 
en donde el punto máximo será la altura máxima 
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TURBOBOMBAS 


de presión (H), la cual se podrá comprobar al 
factorizar el modelo matemático de la bomba 
cuando el primer coeficiente (A) sea negativo 
se tendrá la máxima altura de presión igual a la 
obtenida cuando nuestro caudal sea igual con 
cero 

PALABRAS CLAVE: Factorización, ingenieria 
mecánica, modelos matemáticos 


11 INTRODUCCIÓN 


Es evidente que el desarrollo actual 
de modelos matemáticos, representa una 
herramienta útil, rápida y de bajo costo para 
el análisis de problemas ingenieriles reales. 
La confiabilidad y precisión de tales modelos 
es a la fecha un tema de interés científico. Lo 
anterior, debido a que se pretende que estos 
ofrezcan un resultado de tal manera que ya no 
sea necesario realizar experimentación para 
comparar la información obtenida numérica con 
respecto a la experimental. En este trabajo se 
tomaron los conceptos matemáticos del algebra 
siendo el tema de factorización que al ser 
aplicado en la ecuación de segundo grado que 
se obtuvo al aplicar los conceptos de Mecánica 
de los Fluidos, siendo la evaluación el balance 
de energia y la cinemática en los impulsores. 
Es importante mencionar que en las industrias 
de procesos químicos, petroquímicos y afines 
se utilizan compresores centrifugos para 
aire y gases de cuyo disefo y análisis están 
fundamentados en conceptos de Ingeniería 


Capítulo 5 


MECANICA. En base a la ecuación fundamental de las turbo maquinas que determinó 
Euler y Bernoulli a partir de la cinemática de los Impulsores, la cual determina la energia de 
presión , la energia cinética y la energia potencial, variables que intervienen en el modelo 


matemático de las turbobombas. 


2 | DESCRIPCIÓN DEL MÉTODO 


Ecuación de Bernoulli para una bomba centrifuga: 





Pg 29 Ec. 1 


Figura 1. Esquema de bomba centrifuga 


Donde los términos son: 


H e PP — Pa 
dic 29 Presíon de descarga de la bomba/Energia de presiíon Ec.2 
Vi-VÊ Vá | | ps ad 
di Velocidad de salida/Energia Cinética Ec3 
Z2 — L4 = Zg= Altura geodésica/Energia Potencial Ec.4 


Sustituyendo las ecuaciones (2, 3, 4) en la ecuación (1) obtenemos la ecuación de 
Bernoulli para una bomba centrifuga: 


Pa, Vi 


+ d 
pg 29 


Ec.5 


A continuación los términos de la ecuación de Bernoulli los pondremos en función del 
caudal, a excepción de la 
Energia potencial () siendo este el termino independiente “C”. 


El término que representa la Energia cinética se trabajarã con el principio de 
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continuidad: 


Q=SV Ec.6 
Donde: 
Q = Caudal 
S = Área 
V = Velocidad 


Despejando la velocidad de la ecuación (6): 
Q 
V=— 
S Ec./ 


Y sustituyendo en el término de la energia cinética: 


? 
O e 





= -——+ 


29 29 qgS? 


El coeficiente: en virtud de tener los valores de gravedad y área constantes, su 


coeficiente será “A”, por lo que se tendrá el siguiente término: 


AQ? = Energia Cinética expresada en función del caudal. 
El término que representa la energia de presión se trabajara bajo la 2a ley de Newton 
(Impulso) 
F=mV Ec.8 


F 7 Q Ho. 
po a as O md) 





P9 DP9 PI Pg pg pgs 
El coeficiente E y tener los valores de masa, área y peso específico como valores 
constantes, su coeficiente sera “B”, por lo que se tendrá el siguiente término. 
BQ = energia de presión expresada en función del caudal. 
El término que representa la energia potencial sera el término independiente C. 
Sustituyendo los términos anteriores en la ecuación de Bernoulli obtendremos: 


H, = 40º + BQ+C 


Ec.9 


Siendo esta ecuación el modelo matemático y cuya gráfica es una parabola conocida 
como “CURVA CARACTERÍSTICA DE OPERACIÓN DE LAS BOMBAS CENTRÍFUGAS” 
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Ha = AQ? + BQ+C 














0 Q, Q» Q, Q 
Cuando Q,= O Siendo esta la Altura máxima de la bomba 


DESARROLLO DEL MÉTODO 


Variación del signo del trinomio con la aplicación del valor máximo o mínimo 
del trinomio 


Vy= ax? + bx + c Ec.9.1 


1- Cada término del trinomio se multiplicará y se dividirá por el término: (4a) 


=) 4 bo) 44 
pede (ax) + a (bx) + - (c) Ec.9.2 
2 - Sumar y restar el término: b 2 en la Ec. 9.2 
3 4a? x? + 4abx + 4ac + b? — b? 
po 4a 
— 4a?x2+4abx+b?+ 4ac—b? 


ta Ec. 9.3 


3- Descomponer en la Ec.9.3 el trinomio cuadrado perfecto: (49º x2 + 4abx + b?) 


(2ax+b)2+4ac-b? 
y e — Ec. 9.4 


4a 


4.- encontrar el valor de “x” igualando a cero de la Ec. 9.4 el binomio cuadrado 
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(2ax + b2=0 


Pero como: (20x + b) (2 ax + b)=0 Ec.9.4a 
De la Ec. 9.4a despejamos uno de los dos factores, quedando la siguiente expresión: 


2 ax+ b=0 


b 
Al despejar “x” de la expresión anterior: X= — — 


2d 
5.-encontrar el valor de “y”: 
Como en la Ec. 9.4a el binomio cuadrado se igualo a cero, quedando de la siguiente 


forma: 


— 4ac — b? 
4a 


Una vez obteniendo los valores de “x” y de “y” podemos concluir que: 
a. Cuando el coeficiente(a) es positivo, tiene un valor mínimo para X = — Ea 
por lo que el trinomio tiene un valor mínima, para: 


Ahora bien si el coeficiente (a) es negativo se tendrá un valor máximo para “y 


APLICACIÓN DEL MÉTODO EN EL TRINOMIO CUADRADO PERFECTO: 
MODELO MATEMÁTICO DE LAS BOMBAS CENTRÍFUGAS 


Definiremos las variables hidráulicas de la siguiente manera: la variable “y” será la 


altura de presión “H" y la variable “x” será el caudal “Q” 


H=AQ*+BQ+C (1) 
— 24 22 RR caia Gado 
F= 40) BQO+=0) (2) 


Sumando y restando en (2) 
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pp es 44 QU HAABQRAAC HD" —B 


(3) 
4A 


Siendo 4 A?Q? + 4 ABQ + 4AC + Bº de (3) el trinomio cuadrado perfecto, al factorizarlo (2 
AQ = B)º que al sustituirla en (3), obtendremos: 


ff-= (2AQ+B)2+4AC—B? 
= 4A 


(2 AQ = B)?=0, descomponiendo este binomio cuadrado en sus factores: 


(24Q + B)(24Q + B) = 0, al despejar cualquiera de los dos factores, obtendremos la 
siguiente expresión 


2A0 + B= 0X (5) 
Despejando Q de (5): 


Interpretación de resultados: 


B 


a. H'tendrá un valor mínimo cuando el coeficiente (4) es positivo y Q = — — 


b. Htendrá un valor máximo cundo el coeficiente (A) es negativo cuyo valor má- 
4AC—Bº* 
ximo vale 
4A 
Ejemplo de aplicación 
Una bomba centrifuga que maneja los siguientes caudales (Q) con sus respectivas 


alturas de operación (H) son los siguientes: 

= m lts = 

Q, =0.045- (40)........ H, = 83.26m 
Es DE era DS: aci 

Q, = 010 (100)......Ho = 63.58m 

mé Lts 
Qs = DADO a) PRA «Ha = 11.0/m 
Cuyo modelo matemático será: 


H = -2345 + 0.375 Q+ 87 
Cuando Q=0 


H = -2345 + 0.375(0) + 87 
Por lo tanto: 
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H = 87m 
Aplicando el método de factorización descrito: 
H; tendrá su valor máximo porque el coeficiente (A) es negativo 


Por lo consiguiente: 


B .. 0.375 
ZA  2(-2345) 


= =0.000079957- 
seg 


Q= 8x1075 2 (caudal demasiado pequefio) 


F 
seg 
— 4AC-B2  4(-2345)(87)-(0.375)?  -816060-0.140625  81606060.1406 


44 4(—2345) -— 9380 -—9380 


H = 67.0001499m 


H aproximadamente igual a la que se obtiene al tener un Q=0 


H 


31 CONCLUSIONES 


En base al Modelo Matemático obtenido el cual es un trinomio cuadrado perfecto 
con el cual, podemos analizar y determinar la curva de operación llamada también CURVA 
CARACTERÍSTICA (parábola) así como también sus variables hidráulicas como las alturas o 
presiones así como también sus caudales, en esta ecuación de energia o modelo matemático 
los términos energéticos que intervienen son la energia de PRESIÓN, la energia cinética 
y la energia potencial la cual se incluirá en la ecuación de BERNOULLI. En particular las 
velocidades son las que refleja el comportamiento cinemático de los impulsores y estas a su 
vez ocasionan la fuerza de impulso debido a la segunda Ley de Newton la cual es utilizada en 
las turbo máquinas centrífugas, estas velocidades, son parte importante de las variables que 
intervienen en dicha ecuación para determinar la energia de presión expresada en alturas. De 
antemano este modelo matemático utilizado para obtener la energia es de manera analítica. 
Otra de las virtudes es el modelado y analizar el comportamiento energético de las bombas 
centrifugas, poniendo de manifiesto la aplicación de esta herramienta matemática como es 
en este caso que es el álgebra y cuyo concepto es la factorizaciones de un trinomio cuadrado 
perfecto que desde la secundaria se ensefia la cual es una herramienta fundamental para ser 
aplicado en los conceptos de Mecánica de los Fluidos, materia de especialidad que se da en 


Ingeniería mecánica y carreras afines. 
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RESUMEN: El presente trabajo consiste en 
el análisis cuantitativo a partir de la cinemática 
involucrando a las variables que intervienen en 
los impulsores tipo cerrado que son utilizados en 
las turbos bombas para determinar en base al 
modelo matemático obtenido, el comportamiento 
mecânico energético de los impulsores pudiendo 
asi obtener la energia intercambiada en el 
impulsor de esta máquinas. 

PALABRAS CLAVE: Termodinámica, cinética, 
sistemas termodinámicos, hidrodinámicos 
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11 INTRODUCCIÓN 


Es evidente que el desarrollo actual de 
modelos, representa una herramienta útil, rápida 
y de bajo costo para el análisis de problemas 
ingenieriles reales. La confiabilidad y precisión 
de tales modelos es a la fecha un tema de interés 
científico. Lo anterior, debido a que se pretende 
que estos ofrezcan un resultado de tal manera 
que ya no sea necesario realizar experimentación 
para comparar la información obtenida numérica 
con respecto a la experimental. En este trabajo 
se tomaron los conceptos de Termodinâmica y 
de Mecánica de los Fluidos, siendo la evaluación 
el balance de energia y la cinemática en los 
impulsores. Es importante mencionar que en las 
industrias de procesos químicos, petroquímicos 
y afines se utilizan compresores centrífugos para 
aire y gases de cuyo disefio y análisis están 
fundamentados en conceptos de Ingeniería 
MECÁNICA. En base a la ecuación fundamental 
de las turbo máquinas que determinó Euler a 
partir de la cinemática de los Impulsores, la cual 
determina la energia de presión (H) debido a la 
velocidad tangencial y a la componente axial de 
la velocidad absoluta. Puesto que la energia de 
presión es directamente proporcional al producto 
de las velocidades tanto axial y la componente 


axial. 
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2 | DESCRIPCIÓN DEL MÉTODO 


Arista de salida 
* de los álabes 






= brista de entrada 
de los álabes 





Corte transversal iai 
(b) 


Corte meridional 
(a) 


Figura 1. Rodete de una bomba centrifuga: (a) corte. 


DESARROLLO DEL MÉTODO 


Los sistemas termodinâmicos e hidrodinâmicos son descritos por las leyes de la 
Termodinâmica y de la mecánica de los fluidos, siendo la ecuación de balance de energia la 
misma para estos sistemas como lo indica la primera ley de la termodinámica y la ecuación de 
Bernoulli para la hidrodinâmica, la ecuación de balance de energia para estos sistemas está 
dada por el siguiente modelo matemático (Ec. 1) en base al modelo mecánico del impulsor 
(Fig. 1b) 


Ecit Ec— q = Eç (1) 


Donde E es la Energia Cinética del Impulsor, Ec es la Energia Adquirida en su 
recorrido, q es el Calor debido a la Fricción gaseosa, remolinos y Transferencia de Calor al 
medio Circundante, £.., es la Energia Cinética a la Descarga del Impulsor meridional, (b) 
corte transversal. En este último se han dibujado los triângulos de velocidad a la entrada y 


a la salida. 
Yty —q=Y (1.3) 
Donde: 
Y,=Ec ,= C* +CpT1 (J/Kg) (2) 
Yu = UiCiu - UzCou (JIKg) (3) Ecuación de Euler 
Y, =E,, = CÊ +CpT, (J/Kg) (4) 
2 
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Como no hay intercambio de energia las Ec (2) y (4) son igual a cero, ya que en las 
bombas centrífugas la temperatura a la succión y descarga es constante, por lo que las 
entalpias son iguales: 

HH, 


De la cinemática del impulsor (ver Fig. 1b). Euler obtuvo la ecuación de la energia 
específica Ec. (3) intercambiada entre el rodete y el fluido, siendo esta: 


Ph = E c=UCy— UC (5) 


a. Deducción de la Ecuación de Euler 


Esta deducción se hará con relación a la misma Figura 1.Que representa, como ya 
hemos dicho, el rodete de una bomba centrifuga; pero todo el razonamiento y por tanto la 
fórmula de Euler deducida mediante él, será válido para todas las turbo-máquinas. 

Las tres velocidades c,, u & w están relacionadas segun la mecânica del movimiento 


relativo, por la ecuación vectorial: 


Wi =C —U (6) 


Suponemos que el álabe (o su tangente) tiene la dirección del vector W,, con lo que 
la partícula entra sin choque en el álabe. La partícula guiada por el álabe sale del rodete con 
una velocidad relativa a la salida Ws, que será tangente al álabe en el punto 2. En dicho 
punto, el álabe tiene la velocidad periférica Uz. La misma composición de velocidades de 
la Ecuación (6) nos proporciona la velocidad absoluta a la salida Co, 


C» — W> -|- Us, (7) 


La partícula de fluido ha sufrido, pues, en su paso por el rodete un cambio de 
velocidad de ca c;. Del teorema de la cantidad de movimiento se deduce el teorema 
del movimiento cinético o del impulse o cantidad de movimiento. En efecto, la ecuación 


aplicada al hilo de corriente a que pertenece la partícula de fluido considerada, 


F = oQAv (8) 


Será: 


dF = dQp(C, — E) (9) 


Tomando momentos en la Ec. (9) con relación al eje de la máquina tendremos: 


aiM = dQp(lzc, — lyc4) (10) 
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Que es el teorema del movimiento cinético. 

Donde: 

dM Es el momento resultante con relación al eje de la máquina de todas las 
fuerzas que el rodete ha ejercido sobre las partículas que integran el filamento de corriente 
considerado para hacerle variar su movimiento cinético. 

dQ Es el caudal del filamento. 

2,1 Son los brazos de momento de los vectores cy c, respectivamente. 

Suponemos ahora que todas las partículas de fluido entran en el rodete en un 
diâmetro D con la misma velocidad c,, y salen a un diâmetro D, con la misma velocidad c,. 
Esto equivale a suponer que todos los filamentos de corriente sufren la misma desviación, 
lo cual a su vez implica que el número de álabes es infinito para que el rodete guie al fluido 
perfectamente. Aplicando esta hipótesis Ilamada teoría unidimensional, o teoria del número 
infinito de álabes, al hacer la integral de la Ec. (10) el parêntesis del segundo miembro será 


constante, obteniéêndose finalmente: 


M = Qp(lac, — bic1) (11) 
Donde: 
M es el momento total comunicado al fluido o momento hidráulico. 
Q es el caudal total de la bomba. 


Pero de la Figura 1 b, se deduce fácilmente que 
L=ncosa Y ls =1, Cosa; 


Luego 


M = Qp (rc) cosa, — 1r,C, Cosa) (12) 


Este momento multiplicado por será igual a la potencia que el rodete comunica al 
fluido. Por tanto 


Pw = MN = QoN(r cosa, — rm, cosx,) W,SI (13) 


Donde: 


Coleção desafios das engenharias: Engenharia mecânica 2 Capítulo 6 2 o. 


o Zn 


rad 
= —— — velocidad angular del rodete, — 
60 seg 


Por otra parte, si llamamos Y a la energia especifica intercambiada entre el rodete 
y el fluido, en nuestro caso la energia específica que el rodete de la bomba comunica al 
fluido, y G al caudal másico que atraviesa el rodete, se tendrá en el Sl: 
pie Di sola ta TEM 
RW) = 6 ( 5 ) Tu (=) RA ! 5 Jp (as) 9 (5) Hu (Mm) (14) 
Donde: 


H, Es la altura equivalente a la energia intercambiada en el fluido: 
2 


fa (75) = Ma (==) = Humg (=) 


Igualando las dos expresiones de la potencia de las Ecuaciones (13) y (14) se tiene 


QpY, = Qpw(rc; cosa; — r/C; COS Q4) 


(15) 
Pero 
nw=uU&rao =u, 
C4 COS — Cy g 02 COS, = Coy 
Donde 
o Son proyecciones de c y c sobre u,u, 0 componentes periféricas de las 
14 2uU 1 2 12 


velocidades absolutas a la entrada y a la salida de los álabes. 
Sustituyendo estos valores en la Ecuación (15) y simplificando, se obtiene la 
ecuación de Euler: 


h, = UzCay — Ui Cry (16) 


Se comprueba la Ecuación (3) 
(Ecuación de Euler: Bombas, ventiladores y turbocompresores) 


De la variable Y se pasa a la variable H por la ecuación: 


fu &) =9 (5) H (m) (17) 
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Por tanto, dividiendo los dos términos de la Ec. (16) por se tendrá: 


* | Primera forma de la ecuación de Euler (Expresión en alturas). 


Us CyymlUoCo, 
H,, — + 1 E” Eau (17.a) 
(Ecuación de Euler, primera forma: bombas, ventiladores, turbocompresores, 
turbinas hidráulicas, turbinas de vapor y turbinas de gas: signo [positivo +) maquinas 
motoras y signo [negativo — ] máquinas generadoras; unidades m, SI). 


b. Modelo Matemático del impulsor 


La altura teórica que da una bomba es la expresada por la ecuación (17.a) donde los 
puntos 1 y 2 (ver fig.2) se refieren a la entrada y salida del impulsor. 
H es la altura que el impulsor imparte al fluido. 






= 
e al 
e a SP 


Figura 2 


De la Ecuación. (16) tenemos que: por ser entrada axial se tendrá la expresión 
energética siendo esta la Energia Cinética 


Ka — UsCouy (18) 


La ecuación (18) expresada en alturas: 


H a UsCoy (1 9) 
Uu g 


Como el término 
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Rs | 
(uzCau) =( =) = V* y este al ser dividido entre (9) se obtiene: 


(19.a) 


Cuya expresión es la energia cinética. 
Ahora bien, la expresión anterior se pondrá en función del Caudal para obtener el 
Modelo Matematico de la Energia Cinética: 
Q= AV (20) 


Donde 

Q es caudal que maneja el impulsor 

A es el área del impulsor 

V es la velocidad del fluido que es impulsado por el impulsor 


Despejando: 


V=Q/A (20.9) 


Sustituyendo la Ec. (20.a) en la Ec. (19.a) 





Finalmente obtenemos el MODELO MATEMÁTICO 
De la ENERGÍA CINÉTICA 


H = BQ? 
En donde: 


B es una constante 


Q es el caudal 


Gráfica 
La representación gráfica del modelo matemático resulta ser una parábola, debido a 


que es una ecuación cuadrática. 
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Figura 3. 


31 CONCLUSIONES 


En base al Modelo Matemático obtenido con ello, poder determinar la ecuación de 
energia o término matemático conocida como la energia de Cinetica, la cual se incluirá en la 
ecuación de BERNOULLI. En particular las velocidades son las que refleja elcomportamiento 
cinemático de los impulsores que son utilizados en las turbomáquinas centrífugas de las 
turbomáquinas, estas velocidades, son parte importante de las variables que intervienen en 
dicha ecuación para determinar la energia de presión expresada en alturas. De antemano 
este modelo matemático utilizado para obtener la energia está de manera analítica. Otra de 
las virtudes es el modelado y analizar el comportamiento energético de las turbomáquinas, 
poniendo de manifiesto la aplicación de la herramienta matemática como es en este caso 
la trigonometria que desde la secundaria se enseha para las soluciones de problemas de 
trigonometria para ser aplicada en mecânica vectorial. 

Por otra parte, podemos observar la ley de la conservación de la energia ya que 
la energia cinética de un fluido, se convierte en energia de presión la cual puede ser 


cuantificada o medible por el instrumento de medición Ilamado manómetro de presión. 
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RESUMO: A aerodinâmica experimental é 
essencial para o desenvolvimento de estruturas 
aeroespaciais. Ela complementa a modelagem 
analítica e as simulações numéricas. A principal 
ferramenta para estudos experimentais 
aerodinâmicos é o túnel de vento. Esse trabalho 
busca projetar e construir um túnel de vento de 
bancada, preferencialmente com um escoamento 
laminar na seção de testes. A melhor opção em 
relação a custo e tamanho da estrutura, é um 
túnel de vento de circuito aberto com a seção 
de testes fechada. As estruturas presentes em 
um túnel de circuito aberto e seção de testes 
fechada são: a câmara de estabilização, o bocal, 
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VELOCIDADE 


a seção de testes, o difusor e o ventilador. O ar 
entra na câmara de estabilização, onde a tela de 
honeycomb torna o escoamento mais uniforme. 
O bocal acelera o escoamento e sua geometria 
influência na qualidade do fluxo na seção de 
testes. Nessa estrutura, duas funções cúbicas 
descrevem a superfície. Um ventilador com pás 
retas adiciona energia ao sistema. A simulação 
numérica baseada nas equações de Navier- 
Stokes transientes e incompressíveis estuda 
turbulência e efeitos viscosos no túnel. A malha 
tridimensional estruturada usada na simulação 
leva em consideração os gradientes de fluxo. O 
modelo de turbulência k-w é usado para avaliar a 
camada limite ao longo do túnel. Como resultado 
dos trabalhos em curso, foi concluído o projeto do 
túnel de vento. A construção do sistema permitirá 
a validação de estudos teóricos e numéricos. 
Para a fabricação do difusor e do bocal, moldes 
foram criados para serem cortados em router 
CNC. O túnel de vento permitirá estudar a 
distribuição de pressão na superfície dos corpos 
aerodinâmicos e a visualização do fluxo em baixa 
velocidade usando uma máquina de fumaça. 
PALAVRAS-CHAVE: Túnel de vento, 
aerodinâmica, fluido incompressível, 
fluidodinâmica computacional, router cnc. 


LOW-SPEED WIND TUNNEL PROJECT 


ABSTRACT: Experimental aerodynamics is 
essential for the development of aerospace 
structures. It complements analytical modeling 
and numerical simulations. A wind tunnel is a 
tool used for the first experimental analyses. 
This work seeks to design and build a bench 
wind tunnel, preferably with a laminar flow in 
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the test section. The best option in terms of cost and structure size is an open circuit wind 
tunnel with a closed test section. The structures in an open circuit tunnel and closed test 
section are the stabilization chamber, the nozzle, the test section, the diffuser, and the fan. 
Air enters the stabilization chamber, where the honeycomb-type flow straightener improves 
its uniformity. The nozzle accelerates the flow, and its geometry influences the quality of the 
flow in the test section. In this structure, two cubic functions describe the surface. A fan with 
straight blades adds power to the system. The numerical simulation based on transient and 
incompressible Navier-Stokes equations coupled with k-w turbulence model. The structured 
three-dimensional mesh for the simulation takes into account the flow gradients. The k-w 
turbulence model validated for low-velocity applications describes the propagation of the 
boundary layer. As a result of the current work, the authors completedthe wind tunnel project 
according to specified requirements. The construction of the system will allow the validation of 
theoretical and numerical studies. For the manufacturing of the diffuser and the nozzle molds 
are made with the CNC router. The wind tunnel will allow studying the pressure distribution 
on the surface of aerodynamic bodies and visualizing the flow at low speed using a smoke 
machine. 

KEYWORDS: Wind tunnel, aerodynamics, incompressible fluid, computational fluid dynamics, 
router cnc. 


11 INTRODUÇÃO 


Pesquisas em aerodinâmica são aplicadas em diversas áreas, como: construções (p. 
ex. edifícios e pontes), veículos (p. ex. automóveis e aviões), equipamentos especiais (p. ex. 
aero geradores e paraquedas), dentre outras. Desde meados do século XX, as pesquisas 
na área de aerodinâmica foram motivadas, principalmente, pela necessidade de melhorar 
o desempenho e o consumo de combustíveis de carros e aviões, principalmente, quando 
os fabricantes constataram que os carros perdiam potência pela resistência do ar (Mehta 
and Bradshaw, 1979). As pesquisas na área aerodinâmica evoluíram essencialmente com 
o desenvolvimento dos túneis de vento. Túneis de vento são equipamentos que permitem 
aos engenheiros estudar o fluxo de ar sobre objetos de interesse, forças que atuam sobre 
eles e sua interação com o fluxo. Os túneis de vento são utilizados para investigar teorias 
aerodinâmicas, e facilitar, principalmente, o projeto de aeronaves. Atualmente, a pesquisa 
aerodinâmica tem expandido para outros campos, como indústria automotiva, arquitetura, 
meio ambiente e educação, fazendo com que os testes realizados em túneis de vento de 
baixa velocidade sejam ainda mais importantes (Carminatti and Konrath, 2019). A prática da 
aerodinâmica experimental de baixa velocidade foi fundamental no desenvolvimento de uma 
série de veículos, e outros dispositivos que devem desempenhar suas funções na presença 
de forças geradas por elevados fluxos de ar ou água. 

Os túneis de vento podem ser de circuito aberto ou fechado. Os sistemas de circuito 
aberto não recirculam o ar e são mais viáveis economicamente, pois eles são compostos 


por três estruturas principais: um bocal, uma seção de teste e um difusor. Já em um túnel 
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de circuito fechado o ar recircula dentro do sistema, portanto ele tem mais componentes e 
é uma sistema mais complexo e custoso. Portanto, o projeto é um túnel de vento de circuito 
aberto, onde as simulações de CFD foram feitas e analisadas para compreender o problema, 
com isso, a geometria do túnel foi definida para facilitar a construção das geometrias mais 
complexas. 

Outra tecnologia importante para a engenharia é a fabricação de peças usando a 
router CNC que proporciona precisão e velocidade na fabricação de peças. Usando essa 
máquina os moldes do bocal de difusor podem ser fabricados de forma precisa e eficiente. 


21 METODOLOGIA 


A metodologia do projeto proposto segue as etapas descritas a seguir. 


2.1 Projeto e Simulação das Partes Constituintes do Túnel de Vento 


O objetivo desta etapa será projetar as partes que constituem um túnel de vento 
subsônico na escala 1:10, contendo: bocal, seção de teste, difusor, telas, colmeias e câmara 
de estabilização. Será projetado um túnel de circuito aberto em que o fluido não circula, pois 
não possuirá passagem de retorno. Os túneis de circuito abertos são amplamente utilizados 
fins de instrução e para calibração de dispositivos de fluxo. Para todas as partes constituintes 
do túnel de vento proposto será calculada sua geometria e perdas de pressão. A retidão 
nestes cálculos será crucial para a construção de um túnel de vento preciso e eficaz, que 
atenda o escopo deste projeto. A seção de teste, local em que os corpos a serem testados 
são posicionados para os experimentos, será retangular confeccionada em acrílico com 
iluminação adequada para aquisição de imagens. A seção de teste será projetada de maneira 
que as interferências no corpo de prova sejam mínimas. O difusor terá como finalidade reduzir 
a velocidade com pequenas perdas, tanto quanto possível, convertendo energia cinética em 
energia de pressão. 

A partir das literaturas as seções do túnel e a potência do motor foram dimensionadas. 
O bocal e o difusor serão feitos de fibra de vidro e a seção de teste de acrílico. Simulações de 
CFD foram realizadas com o ANSYS Fluent para uma melhor compreensão do comportamento 
do escoamento. 

Uma das primeiras coisas a ser feita é definir as dimensões da seção de teste e a 
velocidade do escoamento no seu interior, onde seu diâmetro hidráulico será a referência 
para a para projetar o bocal e o difusor. 

Para a fabricação de componentes complexos como o bocal e o difusor, é necessário 
a confecção de moldes feitos em router CNC. Esses moldes serão desenhados usando o 
software SolidWorks. 
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2.2 Câmara de Estabilização 


O processo de uniformização do fluxo começa na câmara de estabilização, reduzindo 
as vorticidades. Uma tela de honeycomb é instalada na entradade de ar do túnel, onde 
possui uma seção de área constante. Esta placa deve seguir três parâmetros para seu 
projeto: o tamanho da célula, a espessura da parede da célula e a espessura da placa do 
honeycomb. 

Inicialmente, a pesquisa de mercado nos permitiu descobrir a disponibilidade do 
honeycomb. Saber o tamanho das células e a espessura da parede permite calcular a 
àrea de superfície do honeycomb. A espessura da estrutura que compõe as células devem 
ser pequena, para que a área aberta estejaentre 95% e 99%, o que resulta em pequena 
queda de pressão, redução da turbulência e do ruído. Outro fator crítico é a relação entre 
o tamanho da célula e a espessura do núcleo. O tamanho da célula deve ser sete vezes 
maior que a espessura do núcleo. Portanto, a célula não pode ser muito grande porque 
a espessura da placa do honeycomb seria muito grande. Assim, uma placa com células 
muito grandes torna a peça muito custosa, pois requer mais material para a sua fabricação 
(Bitzer, 1997). 


2.3 Seção de Testes 


A seção de testes do túnel de vento é um dos parâmetros iniciais do projeto, onde 
a velocidade e qualidade do escoamento são requisitos importantes ao dimensionar a 
estrutura. É nessa parte onde os experimentos são feitos e observados. O formato da 
seção de teste que melhor se adequa para túneis de vento de pequeno porte é a que possui 
uma seção transversal quadrada. Para que o escoamento seja adequado, essa estrutura 
deve ter um comprimento duas vezes o seu diâmetro hidráulico (Barlow et al., 1999). 


2.4 Ventilador Controle de Velocidade 


Inicialmente é necessário ter uma noção da magnitude da potência do sistema. 
Usando a Tabela 1, a potência do ventilador é determinada para a velocidade desejada no 
túnel de vento. Essa potência é relacionada com a área da seção transversal da seção de 
testes. Além da potência, outra característica importante do ventilador axial é o formato da 
pás que devem ser retas (Barlow et al., 1999). 





Tabela 1 — Potência do ventilador por unidade de área para velocidades variadas. 


Fonte: (Barlow et al., 1999). 
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O controle de velocidade do ventilador é feito usando um modulador de largura de 
pulso, conhecido como PWM. Para escolher o melhor controlador é necessário calcular a 


corrente máxima do ventilador usando a Equação 1 apresentada a seguir 


(1) 


em que P é a potência do ventilador, U a tensão de operação e i a corrente. Com o valor da 


corrente máxima do sistema o PWM é escolhido. 


2.5 Difusor 


O difusor deve ter um ângulo de contração entre 2º e 3,5º. Esse ângulo é calculado 
usando a Equação 2 e o comprimento do difusor é calculado usando a Equação 3, respeitando 
o intervalo de O (Barlow et al., 1999). Outro fator importante para o tamanho do difusor, são 
os ventiladores presentes no mercado, que atendam a dimensão estipulada para o projeto. 


Um exemplo de difusor está apresentado na Figura 1. 





0 = arctan (E) (2) 
(58 i 





Figura 1 — Geometria do Difusor. 


2.6 Bocal 


O bocal é dimensionado seguindo a mesma ideia do difusor, porém o ângulo de 
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contração O deve ser de 12º. Além desse ângulo, a razão das áreas da entrada e saída do 
bocal é utilizada para obter o comprimento da estrutura. Para estudos aeronáuticos a razão 
de áreas deve estar na faixa de 8 a 9. Essa razão, N, é calculada usando a Equação 4 onde 


A, é a área de entrada e A, a área de saída. 


- Ae (4) 
= 


Outra consideração importante é o formato do bocal, que deve ser formado por duas 
equações cúbicas (Ahmed, 2013). Geometrias retas do bocal apresentam um escoamento 


menos uniforme na seção de testes (Driss, 2019). 


2.7 Simulação CFD 


A simulação CFD foi realizada com o software ANSYS Fluent, usada para compreender 
melhor o problema e avaliar qualitativamente o escoamento na seção de testes. A simetria 
do meio permite que apenas W do fluxo seja analisado. Dessa forma, a discretização foi 
feita e pode ser observado na Figura 2. Nessa figura é possível observar um refino na malha 


próximo as paredes e quinas onde é formada a camada limite. 





Figura 2 — Malha computacional. 


O modelo usado para resolver as equações de Navier-Stokes é o k-w SST que tem 
um bom compromisso entre o esforço computacional e precisão dos resultados. A diferença 
de pressão de 30 Pa entre a entrada e a saída, foi a condição de contorno utilizada. 
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31 RESULTADOS E DISCUSSÃO 


O resultado final é um túnel de vento de bancada com 2,3 m de comprimento, 
apresentado na Figura 8. 


Câmara de 
Estabilização 





| Ventilador | 


Figura 3 — Túnel de vento. 


3.1 Câmara de Estabilização 


A câmara de estabilização é composta por uma placa de honeycomb feita de alumínio 
3003H18. Essa placa pssui um tamanho de 0,58 m X 0,58 m e uma espessura de 0,04 m. A 
tela de honeycomb possui 14355 células, e cada célula possui uma espessura de parede de 
0,05 10% m, e o tamanho de cada célula é de 5,2 10º m 


3.2 Seção de Testes 


A seção de testes quadrada, na Figura 4, possui um diâmetro hidráulico de 0,2 m e 
seu comprimento é de 0,4 m. Essa peça possui uma abertura lateral de 0,3 m X 0,15 m para 
manuseio de experimentos. 
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Figura 4 — Seção de Testes. 


3.3 Ventilador 


A potência do ventilador calculada é de 24 Watts, para uma velocidade de 8 m/s 
na seção de testes, com uma área transversal de 0,04 m?. Ao pesquisar o mercado de 
ventiladores axiais foram encontrados dois modelos aptos ao problema, um com 0,4 m de 
diâmetro e outro com 0,3 m de diâmetro. Porém, o ventilador com o maior diâmetro tornaria 
o túnel muito grande. Dessa forma, um ventilador com 0,3 m de diâmetro, bifásico, com uma 
tensão de 12 Volts e 80 Watts de potência foi escolhido para o túnel. 

Com a Equação 1 e as informações do ventilador, a corrente máxima é 6,6 Ampere. 
Assim, uma fonte de 12 Volts e 5 Ampere de corrente foi escolhida, junto com um controlador 
PWM de 12 Volts e corrente máxima de 5 Ampere. Esta fonte de alimentação é capaz de 


fornecer 60 W de potência. 


3.4 Difusor 


A partir do diâmetro do ventilador e usando a Equação 3, o comprimento do difusor 
é de 0,954 m. A região de entrada do ar é a parte quadrada na Figura 5, que possui 0,2 m X 


0,2 m de dimensão, e a porção onde o ventilador é fixado têm um diâmetro de 0,3 metros. 
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Figura 5 - Difusor. 


3.5 Bocal 


A razão de contração do bocal escolhida é de 8,41, ou seja, a área de entrada do bocal 
é de 0,3364 m* e a área de saída é de 0,04m?. Portanto, as dimensões da entrada e da saída 
são 0,58m X0,58m e 0,2m X 0,2 m, respectivamente. Usando a Equação 3, o comprimento 
do bocal é de 0,874 metros. A estrutura final do bocal está apresentada na Figura 6. 





Figura 6 - Bocal. 


3.6 Simulação CFD 


O resultado da simulação em regime permanente é apresentado a seguir na Figura 7, 
em que é possível observar o perfil de velocidade dentro da seção de testes. Na simulação, o 
ar entra com uma velocidade aproximada de 1 m/s no bocal e na seção de testes chega com 
uma velocidade média de 8.71 m/s. 
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Velocidade 
8.71 


7.84 
6.97 
6.10 
5.23 
4.36 
3.48 
2.61 
1.74 
0.87 


0.00 
[m's] 





Figura 7 — Velocidade do escoamento. 


Para melhorar o comportamento do escoamento na seção de testes é aconselhável 
colocar chanfros de 45º nas quinas da seção de teste (Barlow, 1999). Porém, a dificuldade 
construtiva do túnel ficaria maior, dessa forma ao avaliar a Figura 8 é possível observar uma 
maior turbulência nas quinas da seção, contudo o interior do escoamento apresenta uma 
velocidade uniforme. Sendo assim, optou-se por construir uma seção de testes sem chanfros 


em sua quinas. 


Velocidade 
8.71 


7.84 
6.97 
6.10 
5.23 
4.36 
348 
2.61 
1.74 
0.87 


0.00 
[m/s] 





Figura 8 — Velocidade do escoamento na seção de testes. 
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41 CONCLUSÃO 


As dimensões e geometrias do tunel de vento foram determinadas, com essa estrutura 
os estudos aerodinâmicos poderão ser realizados para diferentes velocidades, usando o 
controlador PWM. A partir da simulação, a geometria quadrada foi implementada na seção 
de teste. Além disso, é possível observar na simulação a região da seção de teste com as 
velocidades mais uniformes, com isso, a melhor posição para acoplar os experimentos dentro 
do túnel. 
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RESUMEN: Se presenta el modelo y simulación 
de un Mecanismo de Yugo Escocés, este será 
disehado de manera analítica y se llegará al 
Disenio del Modelo en Solidworks. Proporcionando 
el diseho del mismo se estudiarán y comprobarán 
con el software Matlab a través de una interfaz 
SolidWork-Simmechanic para obtener un 
diagrama de control del mismo con la finalidad de 
manipular las variables físicas, de movimiento y 
demostrar la trasformación del movimiento. 
PALABRAS-CLAVE: Disefio, Yugo escocés, 
Simulación, SOLIDWORKS, Ecuaciones de 
movimiento. 
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SOLIDWORKS-MATLAB 


ANALYSIS AND SIMULATION OF A 
SCOTCH YOKE MECHANISM THROUGH 
THE SOLIDWORKS-MATLAB INTERFACE 


ABSTRACT: In this article is presented the model 
and simulation of a Scotch Yoke Mechanism, this 
willbe analytically designed and the Model Design 
will be generated in the Software Solidworks. 
Providing the design of it, they will be studied 
and verified with the Matlab software through 
a SolidWorks-Simmechanic interface to obtain 
a control diagram of it in order to manipulate 
the physical and movement variables and 
demonstrate the transformation of the movement. 
KEYWORDS: Design, Scotch yoke, Simulation, 
SOLIDWORKS, Equations of motion. 


11 INTRODUCCIÓN 


Esta investigación presenta el disefio y 
simulación de un mecanismo de yugo escocês 
realizado en Simulink mediante la interfaz 
SolidWork-Simmechanics con el fin de analizar 
su funcionamiento en las condiciones más 
reales posibles, con el objetivo de que los 
resultados obtenidos sean los más cercanos 
al funcionamiento real del mecanismo y así 
implementar una plataforma de simulación 


confiable para el análisis posterior del 
funcionamiento del prototipo. [1] 

En este trabajo se profundiza sobre la 
realización de dicha plataforma, dando especial 
atención a las ecuaciones de movimiento, modelo 


dinâmico de mecanismo, y la interacción entre 
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ellos formando así un solo modelo completo del mecanismo de Yugo Escocés, visualizando 
como afecta los distintos parâmetros de cada uno de los modelos al comportamiento global 
del sistema y a los parâmetros de interés en cuanto al control. 


21 MECANISMOS 


Un mecanismo es aquella parte de una máquina que contiene dos o más piezas 
dispuestas de tal manera que el movimiento de una de ellas obliga a moverse a las otras, 
según una ley definida que depende de la naturaleza de la combinación [12]. 

Los mecanismos pueden clasificarse de la siguiente manera: 


*- | Mecanismo de Transformación de Movimiento: Es aquel en el que el elemento 
motriz y el elemento conducido tienen distinto tipo de movimiento, por ejemplo, 
Mecanismo de Manivela-Biela-Corredera y Mecanismo de Yugo Escocés. 


*- | Mecanismo de Transmisión de Movimiento: Es aquel en el que el elemento 
motriz y el elemento conducido tienen el mismo movimiento, por ejemplo, Me- 
canismo de Engrane Recto, Mecanismo de Engrane Helicoidal, Mecanismo de 
Polea con Banda en V y Mecanismo de Catarina con Cadena de Rodillos. 


MECANISMO DE YUGO ESCOCÊS 


El presente artículo se base en el mecanismo de Yugo Escocés que es un mecanismo 
de transformación de movimiento, ya que, convierte un movimiento rotatorio en un movimiento 
lineal. La manera en que este mecanismo realiza su movimiento es el siguiente, cuando el 
disco gira mediante un movimiento rotatorio hace que el yugo tenga un movimiento lineal. En 


la figura 1, se muestra una imagen de este mecanismo. 


LE e 15 





Figura 1 Mecanismo de Yugo Escocés de una Máquina de Vibraciones (ITSTR). 


DESCRIPCIÓN DEL MÉTODO 


Para determinar las características cinemáticas, de fuerzas y de esfuerzos, se 
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determinaron ciertas características geométricas del mecanismo de Yugo Escocés, las 


cuales se muestran en la Tabla 1: 


Factor de Disefio 
Posición Angular de la Manivela () 
Velocidad Rotacional de la Manivela (n) 


Tabla 1. Características geométricas del mecanismo de Yugo Escocés. 





Los diversos análisis que se realizaron en el mecanismo de Yugo Escocés, fueron 
los siguientes: 


* Análisis Geométrico. En esta parte se dan las ecuaciones para determinar la 
masa del Disco. 


* Análisis Cinemático. En esta parte se dan las ecuaciones para determinar la 
posición, velocidad y aceleración del Yugo, así como la velocidad del Disco. 


* Análisis de Fuerzas. En esta parte se da la ecuación para determinar la fuerza 
que soporta el perno. 


* Análisis de Esfuerzo. En esta parte se da la ecuación del esfuerzo que soporta 
el perno, además de la teoría de disefio para determinar el material del perno. 


* || Disefio y Modelado del Mecanismo en Solidwork. En esta apartado se conjuntan 
los resultados obtenidos de los análisis y la modelación del mecanismo en 3D. 


* | Simulación. En este apartado se obtiene el diagrama de control del mecanismo 
a través de la interfaz con Solidwork-Simmechanic para la manipulación de las 
variables físicas y verificar la trasformación del movimiento. 


Análisis Geométrico 
Para determinar el volumen del disco V se emplea la Ecuación (1), donde r es radio 
del disco y t es el espesor del disco, de esta forma: 
F= rt 
= m(5)*(0.5) 


= gê 
V = 39.27 plg Ecuación (1) 
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Para determinar la masa del disco m se emplea la Ecuación (2), donde p es la densidad 
del material y V es el volumen del disco, de esta forma: 
m = pV 
= (.283(39.27) 


n = 32.17 slugs lÓ 
; slugs Ecuación (2) 


Análisis Cinemático 
Para el análisis cinemático, se utilizaron las Ecuaciones (3), (4) y (5) donde X es la 
posición, v es la velocidad lineal, A es la aceleración lineal, r es el radio, w es la velocidad 


angular. Por lo tanto: 


X=r(1-—-cos0)+ (r?/2L)(sen20) Ecuación (3) 
 =WT Ecuación (4) 
de=: - bd Ecuación (5) 


En la Tabla 2, se muestran los resultados de la posición, velocidad y aceleración del 
yugo escocés para las dos formas que se deben enlazar en el disco y el yugo, a decir, a un 
radio de 3 plg y a uno de 4 plg: 


Radio del Perno (r) 3 plg 4 plg 


Posición del Yugo (xY) o Brg 8 plg 


Velocidad Lineal del Yugo (vY) oO 


Tabla 2 Resultados del análisis cinemático del mecanismo de Yugo Escocés. 





Análisis de Fuerzas 
A partir del análisis cinemático el siguiente paso es el análisis de fuerzas que actuará 
en el perno y se aplica la segunda ley de Newton que dice que la magnitud de la aceleración 
es proporcional a la magnitud de la fuerza e inversamente proporcional a la masa de la 
partícula [24]. Si una fuerza F se aplica a una partícula de masa m, la segunda ley de Newton 
expresa mediante la Ecuación (6): 
F=ma Ecuación (6) 


En la Tabla 3, se muestran los resultados del análisis de fuerzas en el Yugo Escocês 
para las dos formas que se deben enlazar en el disco y el yugo, a decir, a un radio de 3 plg 


y a uno de 4 plg: 
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3 plg 4 plg 
32.17 slugs 32.17 slugs 


82.25 plg/seg2 109.66 plg/seg2 
28.5 Lbs 38 Lbs 


Tabla 3. Resultados del análisis de fuerzas del mecanismo de Yugo Escocés. 





Análisis de Esfuerzos 

Al estar en funcionamiento el mecanismo de Yugo Escocés, es necesario un análisis 
de esfuerzos ya que el perno que sujeta el Disco y el Yugo está sometido a Esfuerzos de 
Corte Normal. Por lo tanto, es necesario determinarlo junto con el Factor de Disehio que es 
una parte fundamental para decir que la pieza es segura. De esta forma, el Esfuerzo de Corte 
Normal se obtiene a partir de la Ecuación (7) donde ?P es la fuerza que soporta el perno des 
el diâmetro del perno: 


4P 

O =— e 
md? Ecuación (7) 
Para determinar que el perno sea seguro se aplica la teoria de la Energia de 
Distorsión mediante la Ecuación (8), donde 5, ES la resistencia de cedencia del material 


y nes el factor de diseho: 





n=—? Ecuación (8) 


En la Tabla 4, se muestran los resultados del análisis de esfuerzos en el Yugo Escocés 
para las dos formas que se deben enlazar en el disco y el yugo, a decir, a un radio de 3 plg 
y a uno de 4 plo: 


Radio del Perno (r) 3 plg 4 plg 
Fuerza en el Perno 28.5 Lbs 38 Lbs 


145.15 Psi 193.53 Psi 
Esfuerzo Permisible del Perno 290.3 Psi 580.59 Psi 


Tabla 4. Resultados del análisis de esfuerzos del mecanismo de Yugo Escocés. 





De acuerdo con los resultados obtenidos de la Tabla 4, el material que se debe 
utilizar para el perno debe de ser un material que soporte un esfuerzo de cedencia de 
580.59 Psi, este material es un ACERO AISI 1040 Rolado en Caliente ya que su esfuerzo 
de cedencia es de 42000 Psi. 
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Simulación Mecanismo de Yugo Escocés 
La primera fase de la simulación es el modelado del dispositivo al que va a ir 


acoplado el mecanismo, este se realiza mediante un software de diseho. En la Figura 1, se 


muestra el Modelo 3D Generado a través de Solidworks. 





Figura 1. Mecanismo de Yugo Escocés modelado en Solidworks. 


Una vez que se culminó el modelo 3D se realiza la interfaz a Matlab con la ayuda 
de la herramienta Simmechanic, para poder obtener un diagrama de control con el cual se 
pueden interactuar con las diferentes variables físicas del mecanismo. En la Figura 2, se 


muestra el diagrama de control del ensamble del mecanismo: 
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Figura 2. Diagrama de control mecanismo de Yugo Escocés simulado a través de Matlab. 


En la Figura 3, se muesitra la validación del Modelo Simulado con Matlab: 
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BE] EnsamblePrototipo0x2820x29 


View Simulation Model Help 
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Figura 3. Mecanismo de Yugo Escocés simulado a través de Matlab. 


-En las Figuras 4, 5 y 6 se pueden observar el comportamiento del mecanismo 
y la transformación del movimiento de rotatorio del disco a un movimiento lineal con la 


implementación de un Yugo Escocês: 


E] Scope3 
aolavialiEk Dá 





Figura 4. Gráfica que muestra el movimiento del eje al disco del Yugo Escocés. 
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Figura 5. Gráfica que muestra la transición del movimiento del disco al perno el cual le trasmite el 
movimiento al Yugo a través del contacto en sus paredes. 





Figura 6. Gráfica que muestra la transición del movimiento de la corredera a la superficie de 
desplazamiento se observa un deslazamiento lineal en la gráfica. 


31 COMENTARIOS FINALES 


En este artículo se presentó un modelo en Simulink de un mecanismo Yugo Escocés 
con el objetivo de observar el rendimiento de este a través de diferentes parâmetros que 
se muestran con la ayuda de la simulación en el entorno Simulink. Se mostró un diagrama 


de bloques junto con la parte matemática en la cual se basó el diseho para la elección de 
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materiales y simulación para el comportamiento con el fin de abarcar de mejor manera 
el funcionamiento del sistema. Con respecto a los objetivos propuestos, estos se han 
alcanzado en su totalidad teniendo en cuenta que la simulación responde de manera 
satisfactoria a las condiciones de simulación pre-establecidas. 

Considerando los resultados obtenidos se disehió un mecanismo Yugo-Escocés y se 
validó con la interfaz SolidWork-Matlab para manipular las condiciones de cada elemento, 
y determinar si todos los aspectos necesarios para su función son correctos y se mostraron 
los resultado de la simulaciones y su diagramas de control obteniendo así para un óptimo 


funcionamiento como pauta para mejor en el âmbito del diseno y el control de mecanismos. 


41 RECOMENDACIONES 


Se recomienda darle seguimiento a este proyecto con el disefo de otros mecanismos 
que se puedan montar en el prototipo, que pueda ayudar a alumnos a tener un mejor 
entendimiento de los Mecanismo de Transmisión de Potencia ya que su finalidad desde 
siempre ha sido ayudar a alumnos de Ingeniería Mecânica y áreas a fines en la Materia de 
Análisis y Síntesis de Mecanismo, ya que con la interfaz que se muestra ahora no solo se 
puede operar si no también llegar a un control y que este sea aprovechado en todos los 


aspectos de la Ingenieria. 
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RESUMO: Muitos estudos que realizam 
simulações estruturais de ossos utilizando o 
método dos elementos finitos consideram o 
tecido ósseo como sendo um material isotrópico, 
porém este é um material ortotrópico de acordo 
com a literatura. O presente trabalho faz um 
comparativo entre esses dois comportamentos 
mecânicos, com o objetivo de analisar se é viável 
utilizar este tipo de aproximação, levando em 
consideração os parâmetros de deslocamentos 
dos sistemas e o tempo de duração da 
simulação. Para isso, é utilizado uma fratura 
coronal do côndilo femoral que foi fixada com 
uso de 2 parafusos esponjosos anteriores e uma 
placa reta do tipo DCP com uso de 4 parafusos 
na face lateral do fêmur distal. Observou-se que 
os resultados de deslocamentos não tiveram 
grandes variações entre o osso isotrópico e o 
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ortotrópico, com exceção dos deslocamentos 
dos parafusos esponjosos, onde houve uma 
diferença expressiva de 106,7% e 135,1%. Então 
é indicado que sejam utilizadas propriedades 
mecânicas isotrópicas apenas para resultados 
iniciais. Para as conclusões finais é mais 
apropriado considerar o osso como material 
ortotrópico, mesmo que prolongue o tempo 
de simulação, fazendo com que não haja uma 
diferença muito alta nos resultados, pois deste 
modo acabaria fugindo da realidade. 
PALAVRAS-CHAVE: Método dos Elementos 
Finitos. Ortotrópico. Fratura Coronal do Côndilo 
Femural. Biomecânica. 


COMPARISON BETWEEN A BONE 
WITH ISOTROPIC AND ORTHOTROPIC 
MECHANICAL PROPERTIES BY THE 
METHOD OF THE FINITE ELEMENTS IN A 
FRACTURE 


ABSTRACT: Many studies that perform 
structural simulations of bones using the finite 
element method consider bone tissue to be an 
isotropic material, however this is an orthotropic 
material according to the literature. The present 
work makes a comparison between these two 
mechanical behaviors, in order to analyze 
whether it is feasible to use this type of approach, 
taking into account the displacement parameters 
of the systems and the duration of the simulation. 
For this, a coronal fracture of the femoral condyle 
is used, which was fixed using 2 cancellous 
screws and a straight DCP type plate using 4 
screws on the lateral face of the distal femur. It 
was observed that the results of displacements 
did not have great variations between isotropic 
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and orthotropic bone, with the exception of displacements of cancellous screws, where there 
was a significant difference of 106.7% and 135.1%. So it is indicated that isotropic mechanical 
properties are used only for initial results. For the final conclusions, it is more appropriate to 
consider bone as an orthotropic material, even if it prolongs the simulation time, making sure 
that there is not a very high difference in the results, as this would end up escaping reality. 
KEYWORDS: Finite Element Method. Orthotropic. Hoffa's fracture. Biomechanics. 


11 INTRODUÇÃO 


Os estudos para o melhoramento e desenvolvimentos de novos implantes vem 
crescendo progressivamente com a utilização da tecnologia de simulação computacional 
pelo método de análise dos elementos finitos (ANSYS, 2019). Muitos estudos utilizam uma 
aproximação das propriedades mecânicas do osso para isotrópica, exemplos podem ser 
vistos com Roseiro et al. (2014) e Duarte e Ramos (2014). Porém, o osso é composto por sua 
maioria de hidroxiapatita, colágeno, proteínas e água, isso faz com ele seja um compósito, 
possuindo propriedades mecânicas ortotrópicas (CURREY, 2016). 

Um elemento de um corpo com propriedades mecânicas isotrópicas apresenta as 
mesmas propriedades em todas as direções observadas, isso ocorre devido a distribuição 
aleatória de seus constituintes. Por outro lado, um elemento ortotrópico possui as mesmas 
propriedades mecânicas em cada plano apenas, considerando 3 planos ortogonais no 
espaço, esse comportamento pode ser observado em madeira, chapas de aço laminadas 
em apenas uma direção e em materiais compósitos por exemplo. Um comportamento 
anisotrópico é caracterizado por possuir propriedades dissemelhantes em todas as direções, 
não se limitando apenas nas direções dos planos ortogonais como ocorrem com os materiais 
ortotrópicos (BORESI; SCHMIDT, 2003). 

No presente trabalho é feito um comparativo dos resultados obtidos de uma 
simulação estrutural de uma fratura utilizando as propriedades mecânicas do osso com os 
comportamentos isotrópico e ortotrópico, para melhor compreensão de suas particularidades. 
Com isso, é verificado se as divergências nos resultados de seus deslocamentos são 
relevantes e se o tempo de duração da simulação justifica essa aproximação. 

Para a comparação dos resultados do osso isotrópico e ortotrópico é utilizado um 
estudo de caso de um homem de 44 anos, que sofreu um acidente de moto, dando origem 
a fratura coronal do côndilo femoral conhecida como fratura de Hoffa, conforme Fig 1. 
Essa fratura fica localizada na parte distal do fêmur, na região do joelho onde o fêmur tem 
uma forma trapezoidal, região cnamada de côndilo, aonde o fêmur se articula com a tíbia 
(TORTORA, 2013). 
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Figura 1. Fratura de Hoffa. 


Esse tipo de fratura é comumente causado devido a traumas de alta energia, sendo 
uma fratura difícil de ser identificada se não for realizada uma Tomografia Computadorizada 
(TC), devendo ser realizado um tratamento cirúrgico criterioso para se evitar complicações 
como lesão articular, principalmente na inserção de parafusos, ou perda de redução por 
fixação insuficiente (RUEDI; MURPHY, 2001) 


21 METODOLOGIA 


Na modelagem da placa e a realização das simulações, foi utilizado um computador 
com processador Intel? Core”” ;5-8250U CPU Q 1.60GHz 1.80 GHz, memória RAM 8,00 
GB, sistema operacional de 64 bits, processador com base em x64, Windows 10 Home 
Single Language, NVIDIA GeForce MX150. Para a modelagem da óriese e dos parafusos foi 
utilizado o programa SolidworksO), versão 19.0, da empresa Dassault Systêémes S.A e para as 
simulações estruturais pelo método dos elementos finitos foi utilizado o programa ANSY SO. 

Foi realizada uma Tomografia Computadorizada (TC) para a criação do modelo 
geométrico do osso fraturado, utilizando o programa Invesalius? (v3.1.1, Centro de 
Tecnologia da Informação Renato Acher — CTI, Campinas, Brasil), no qual foi obtido o objeto 
3D por Mendonça (2018) em sua dissertação. Após isso, o modelo foi tratado no programa 
MeshmixerQ), versão 3.5 para Windows, da empresa Autodesk, reduzindo sua a qualidade, 
inspecionando e corrigindo seus defeitos, e por fim transformando-o em solido, com a 
finalidade de otimizar o modelo para a montagem e simulação. 

Após isso, foi realizado o redesenho da placa reta do tipo DCP e dos parafusos, sendo 
quatro parafusos bloqueados de comprimento 16 mm, 26 mm, 26 mm e 28 mm, todos com 
diâmetro de 3,5 mm, e dois parafusos esponjosos de comprimento 50 mm e 60 mm, ambos 


com diâmetro de 4 mm, respectivamente conforme a Fig 2. 
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Figura 2. Placa reta e parafusos 


Em seguida, é feita a modelagem da placa no osso e, após isso, são inseridos 
os parafusos no sistema, sendo os bloqueados fixados na placa e no osso enquanto os 
esponjosos apenas no osso. Na Figura 3 é possível observar a montagem da placa e dos 
parafusos no osso fraturado. 





Figura 3. Montagem da placa e dos parafusos no osso fraturado 


Para a realização da simulação estrutural, os valores dos módulos de elasticidade 
e coeficientes de Poisson utilizados por Currey J. (2016) em seu capítulo de livro, no qual 
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realiza um estudo bibliográfico de testes feitos em osso cortical. Para o presente estudo, foi 
utilizado os resultados obtidos de ensaio mecânico de compressão, visto que é a situação 
de trabalho de um fêmur na parte distal. Para os dados de entrada dos dados de engenharia 
do AnsysêQ, é considerado as direções X e Y como radias e Z como axial. Os dados foram 
preenchidos conforme consta na Fig. 4 para a elasticidade isotrópica e na Fig. 5 para a 
elasticidade ortotrópica. Na Figura 6 estão as propriedades mecânicas da placa e dos 
parafusos, no qual é utilizada a própria biblioteca do Ansys, considerando-os como stainless 
steel (aço inoxidável). 
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Figura 5. Propriedades mecânicas do osso ortotrópico 
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Figura 6. Propriedades mecânicas do aço inoxidável 


Para as condições de contato entre os corpos, é considerado bonded (sem graus de 
liberdade) em todas as faces que estão em contato com os parafusos bloqueados, sendo 
a placa, o fragmento e o fêmur. Também é considerada essa condição para as faces dos 
parafusos esponjosos que tocam o fragmento. A condição frictionless (sem atrito) é usada 
na parte do corpo do parafuso esponjoso que toca o fêmur. As faces da cabeça do parafuso 
esponjoso e da placa que estão em contato com o osso, é considerado frictional (com atrito) 
com coeficiente de atrito igual a 0,2. 

Na etapa de discretização do modelo, é utilizada a função sizing para definir os 
elementos dos parafusos bloqueados e da placa com tamanho de 2 mm e dos parafusos 
esponjosos com tamanho de 3 mm. Depois disso, é utilizada função face meshing nos 
parafusos esponjosos para melhor distribuição, aplicando nas faces da rosca e nas faces da 


cabeça paralela a ela. A malha resultante pode ser observada na Fig 7. 
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Figura 7. Malha resultante da placa e dos parafusos 


Para o fragmento e o fêmur, não foram realizados refinamentos, apenas é utilizada 


a malha automática que o próprio software gera. A malha do sistema inteiro com todos os 
seus componentes está na Fig. 8. 
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Figura 8. Malha do modelo com todos seus componentes 


É considerado um uma carga de 4000 N de compressão na extremidade do osso, 
conforme Taylor, Heller, Bergmann e Duda (2004), a carga atuante no joelho é igual a 5,1 
vezes o peso do corpo, levando em consideração uma pessoa subindo uma escada, então 


considerando uma pessoa de 80 Kg é encontrada a carga apresentada anteriormente. O local 
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e a direção em que é aplicada essa carga foram escolhidas devido ser o caso mais crítico, no 
qual ocorrem maiores deslocamentos devido a força estar sendo aplicada paralelamente ao 
plano da fratura, conforme verificado por Loureiro e Ramin (2019). 

Além da carga de compressão que atua no osso, a extremidade proximal do modelo 
foi fixada e aplicou-se uma pré-carga nos parafusos esponjosos de 2500 N, utilizando a 
média dos valores de trabalho 2000 N e 3000 N, conforme Perren (2002). Todas as condições 
de contorno no modelo podem ser observadas na Fig. 9. 
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Figura 9. Condições de contorno do modelo 


Para uma melhor comparação dos resultados, foram escolhidos quatro pontos 
nos dois modelos: dois nas extremidades inferior e superior da placa e dois localizados 
no plano da fratura, sendo um localizado no fragmento e o outro no fêmur. Os pontos 
informados podem ser observados nas Fig. 10. Além desses pontos, foram selecionadas as 
faces superiores dos parafusos esponjosos para realizar a comparação dos deslocamentos 
máximos destas áreas também. 
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Figura 10. Ponto 2 localizado na extremidade superior da placa de osteossíntese 


31 RESULTADOS 


Ambos os resultados possuem valores próximos para os deslocamentos máximos do 
sistema, com variação de 5,2% em relação ao osso ortotrópico. Além disso, o deslocamento 
ocorre exatamente no mesmo local, na extremidade inferior da placa de osteossiíntese, 
conforme Fig. 11 e 12. 
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Figura 11. Deslocamento máximo para osso isotrópico 
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Figura 12. Deslocamento máximo para osso ortotrópico 


É notável que o gradiente de deslocamento entre os dois são muito semelhantes, porém 
há uma divergência nos resultados dos parafusos esponjosos. Então foram selecionadas as 
faces superiores desses dois parafusos para uma melhor análise dos resultados. Na Tabela 
1 consta o comparativo dos valores dos deslocamentos máximos nas faces selecionadas do 
osso isotrópico e ortotrópico. Observa-se que há uma variação grande dos resultados obtidos 
dos deslocamentos desses parafusos entre os dois materiais, sendo 106,7% para o parafuso 


superior e 135,1% para o parafuso inferior. 


Parafuso Deslocamento osso Deslocamento osso Diferença | Diferença em relação ao 
esponjoso isotrópico (mm) ortotrópico (mm) (mm) osso ortotrópico (%) 


Superior 0,884 0,428 0,456 106,/% 
Inferior 1,397 0,594 0,803 135,1% 


Tabela 1. Deslocamentos das faces superiores dos parafusos esponjosos 





Na Tabela 2 é verificado os valores dos deslocamentos referentes aos quatro pontos 
selecionados citados anteriormente e a comparação entre eles. Nota-se que as diferenças 
dos valores dos deslocamentos são muito inferiores aos deslocamentos que ocorrem nos 


parafusos esponjosos. 


Deslocamento osso | Deslocamento osso | Diferença | Diferença em relação ao 
isotrópico (mm) ortotrópico (mm) (mm) osso ortotrópico (%) 





Tabela 2. Deslocamentos dos 4 pontos escolhidos 
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O tempo de duração da simulação da análise do osso isotrópico foi de 10 minutos 
e 54 segundos enquanto do osso ortotrópico foi de 26 minutos e 29 segundos, sendo uma 
diferença de tempo de 15 minutos e 35 segundos, que representa uma porcentagem de 
58,8% do resultado do material ortotrópico e 143% do resultado do material isotrópico. 


4 | CONCLUSÕES 


Os resultados da distribuição de deslocamentos de ambos os casos possuem 
resultados semelhantes, porém acaba divergindo para os resultados dos parafusos 
esponjosos. Isso se dá pelo fato das propriedades mecânicas nas direções radiais serem 
diferentes, afetando o resultado final. 

Levando em consideração o estudo realizado por Loureiro e Ramin (2019), que 
utilizaram o mesmo modelo e as análises foram feitas em pontos localizados no plano da 
fratura do osso, os valores dos resultados não tiveram uma variação significativa, pois o 
principal objetivo era analisar a estabilidade mecânica da osteossíntese, ou seja, analisar O 
deslocamento dos fragmentos da fratura. Caso fosse feita a análise nos parafusos esponjosos, 
os resultados teriam divergidos da realidade, visto que uma simulação é uma aproximação. 

Então, para uma interpretação inicial dos resultados e otimização do tempo, é 
interessante utilizar uma aproximação para o osso isotrópico, visto que usualmente são 
realizados ajustes antes de chegar no resultado final, pois é necessário fazer alterações e 
considerações durante o processo. Após os ajustes necessários, é indicado que altere as 
propriedades do osso para ortotrópica para obter um resultado mais próximo do real. 
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ABSTRACT: Computer  simulations are 
performed in order to analyze systems dynamic 
behavior. However, some subcomponents may be 
too complex to be modeled, which occurs due to 
nonlinearities and coupling of subsystems. In this 
scenario, a hardware-in-the-loop (HiL) test can be 
done, which consists of separating the system into 
two subsystems: a physical one, containing the 
complex component, and a numerical one, to be 
simulated. The numerical subsystem calculates 
the variable to be imposed on the hardware, 
such as a displacement. The hardware responds 
to the input and its response is measured and 
sent back to the numerical subsystem. Thus, 
the HiL test is a real-time closed-loop system, 
which has been used in the automotive area in 
applications such as anti-lock brake systems and 
semi-active suspensions. In order to close the 
loop, a transfer system composed of actuators 
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and sensors is required. Since it is a dynamic 
system, the transfer system cannot respond 
instantaneously to a command, causing a time 
lag that deteriorates the closed-loop response 
and may lead to an unstable condition. In this 
paper, a HiL test is performed on a mass-spring- 
damper system using the spring as the hardware. 
Methods for delay compensation are presented, 
simulated and compared experimentally. 
KEYWORDS: hardware-in-the-loop, 
dynamics, delay 


control, 


11 INTRODUCTION 


In the development stage of new 


technologies, system simulation before its 
conception is indispensable, leading to lower 
cost with prototypes and safer product validation. 
However, a mathematical model of the system of 
interest may not be easily obtained, since non- 
linearities and complexities may be present. In 
this context, a Hardware-in-the-Loop (HiL) test 
may be employed, which consists of dividing 
the overall systems into two subsystems, one 
that is a numerical model (the software) and the 
other that is the physical component of interest 
(the hardware). With actuators and sensors, the 
two subsystems are coupled and a real-time 
simulation can be employed, thus analyzing 
the overall system dynamics. Figure 1 shows a 
generic Hil simulation block diagram. Signal is 
the external situation imposed in the system, 
like a force or displacement, and signal is the 
physical subsystem's measured response. Signal 
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is the numerical subsystem input, which is calculated based on signals and . The software 
output is converted to a continuous signal through a D/A converter which is then sent to 
an actuator, who will impose it to the physical component through signal . The component's 
response is measured by a sensor, yielding , which is converted through an A/D converter to 
. From this analysis it is clear that the HiL test is a closed-loop real-time system. 
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Figure 1: Typical HiL test block diagram. 


HiL testing begun in aerospace industry, through flight simulators and missile location 
systems (BACIC, 2005). In the past 20 years, the technique has been used mainly in automotive 
applications. Svenson et al. (2009) developed a HiL simulator in order to test heavy trucks 
stability control systems, where the braking system and its controllers were the hardware and, 
as software, the truck's model was employed. One of the main applications of HiL is within 
suspension systems testing. Hong et al. (2002) used HiL test to evaluate a Skyhook control 
method, developed by the authors, for McPherson semi-active suspensions. The hardware 
in this situation was the variable damper, which is a very non-linear component, while the 
controller module and the other suspension components were the modeled software. The 
main problem that arise from HiL simulations is a time delay inherent to the interface between 
software and hardware. Since the actuator cannot respond instantaneously to a command, 
there is an inevitable time delay which is inside the closed-loop system. Hence, the response 
presents an error and, depending on the amount of delay and in dynamic characteristics of 
the subsystems, it may also lead to instability. Delay compensation methods were proposed 
and many works on this area were performed in order to improve HiL results, but no method 
is perfect. The Forward Prediction Method (WALLACE, 2005) adjusts a polynomial to the past 
calculated output values of the software and extrapolates this curve, forward predicting the 
signal's value in the future, in a time equal to the time delay. Thus, the closed-loop system is 
compensated. 

Another method was proposed by Smith (1959), but not for the HiL context. Known as 
the Smith Predictor, the technique uses a model of the system under test in order to predict 
its future response, thus compensating the time delay. However, a model of the plant is not 
available when HiL test is employed. By employing an approximation of the plant, errors arise, 
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thus this is not the most indicated method. A similar idea was proposed by Gawthrop et al. 
(2008) for real-time simulations, where a model is needed but, again, this is not the best 
solution for this application. A different approach where an analysis of the discrete operations 
performed by the acquisition boards was employed in order to characterize the time delay 
(CARRILLO, 2012). By choosing the simulation's sample time as larger than the time delay, 
the author noticed that the equivalent delay becomes equal to one time step, which allows for 
a different kind of correction based on the board's tasking order alteration. However, a high 
sampling period may not be suitable for some applications, leading even to stability issues 
(CARRILLO, 2014). 

In this paper a Hardware-in-the-Loop test is applied to a mass-spring-damper system. 
In order to investigate how the time delay affects the system's dynamics, the hardware is a 
spring, because it can be easily identified and modeled. Delay compensation methods are 
presented, applied and, because the hardware is parametrized, it is possible to compare the 
results to a reference response, thus validating the experiment. A new compensation method 
is then proposed and compared to the previous ones. Numerical simulations are run by using 
a virtual spring and, then, experiments are conducted with the real spring by using MATLAB/ 
Simulink software, dSPACE acquisition board and SCHENCK hydraulic actuator. 


2 | EXPERIMENTAL SETUP 
The following equipment was used: 
. MATLAB/SIMULINK software, to model the numerical subsystem; 


* | dSPACE DS-1102 board, which contains A/D and D/A converters and the Con- 
trolDesk software to monitor the signals; 


* — Instron/SCHENCK actuator system, used for displacement control, equipped 
with a LVDT sensor for displacement sensing and load cell for force sensing, 
shown in Figure 2. 





Figure 2: Instron/SCHENCK actuator. 
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The dSPACE board works as a zero-order holder system, which means that the 


digital input signal in the D/A board is held constant during one sampling interval. 


31 SYSTEM MODELING 

A mass-spring-damper system is analyzed, with the spring being the physical 
component. 
3.1 Full mass-spring-damper model 

The mass-spring-damper system is indicated in Figure 3a, where is the weight force 
in N, which refers to mass in kg, is the stifífness of the spring in N/m and is the viscous 
damping parameter in Ns/m. Figure 3b shows a free-body diagram for this problem, where 


and are the forces due to the spring and damper, respectively, in N. 














(a) (b) 


Figure 3: Model of the system under analysis. (a) Mass-spring-damper, (b) free-body diagram. 


The equation of motion is then given by: 


my+cy+ky=F-—b. (1) 


In order to represent the system in state-space form, these equations need to be 
re-written in terms of states as indicated by Equations 2 and 3, where is the state vector, 


contains the system's outputs and represents the inpuis. 


x = Ax+ Bu, (2) 
y = Cx + Du. (3) 


The following states are adopted: x, = y and x, = y. The system's input is U = = Iy 
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while the output is displacement, Y = X1. This leads to: 


0 1 O 1 
E —hk im sétml RT il U, (4) 
y=[1 0]x+[0u. (5) 


3.2 HiL model with physical spring 


In this case, the term F,, corresponding to the spring restoring force, is modeled as an 
input, which implies that u = É — É, — kw. Then, Equation 1 can be re-written in state-space 
form, adopting the same state vector, as: 


+=[9 ccjm)**[1/m) (6) 


y=[1 o0lx+[0u. (7) 


3.3 Hardware identification 


Knowledge about the spring, sensor and delay characteristics is needed in order to 
validate the HiL simulation results. A displacement sine wave was imposed to the spring and 
the force/displacement behavior is shown in Figure 4(a). Since non-linearity was noted, a third 
order polynomial was adjusted to the curve, which is described by Equation 8, where the unit 
of is kN and the unit of is mm. Note that for the spring is under compression. 


Fn = —1.28x107'yº — 3.441x107*y* + 0.01238y — 0.01173. (8) 


The system delay is defined as the time needed by the actuator to execute a 
commanded displacement and it was modeled as a transport delay, as is commonly done 
in the literature, although other methods such as Padé approximation could be used. By 
inputting a sine wave and measuring the response, it was possible to quantify the time delay 
as e 6s. Figure 4(b) shows the actuator's behavior due to the input sine wave with a zoom 
applied in order to visualize it. 
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Figure 4: System characteristics. (a) Force and displacement behavior of the spring under test, (b) 
Delay between commanded and executed displacement. 


Finally, in order to determine the gains involved in the voltage to displacement 
conversion and vice-versa, a voltage sine wave was generated in SIMULINK, sent to the 
D/A board with the help of ControlDesk software and then sent to the actuator, which led to 
a 40 mm/V gain. Then, by measuring force rather than displacement, a 10 kN/V gain was 
determined for the force to voltage conversion. 


41 DELAY COMPENSATION METHODS AND SIMULATIONS 


Three different methods are presented and simulated in this work. The following 
subsections describe their functionalities. 


4.1 Tasking order sequence alteration method 


The response delay is exactly T, when T,. T, where T, is the sampling time. The 
discrete state equations are given below, where the subscript in the matrices represent their 
discrete equivalent (CARRILLO, 2014). 


x[k E 1] — AgX|k | + Bqulk|, (9) 
v|k] = Caxlk] + Dqulk). (10) 
By analyzing the output equation in time step X, it can be seen that the calculation of y 


[X] depends on x[k] and u[k]. For this implementation, the tasking order states that vector 
vIk] = f(x[k] u[k]) is calculated before the state update equation x[k + 1] = f(x[k], u[k]). 
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However, as there is one time step delay when measuring the physical subsystem's reaction 
force, the state equations input is ulk| = Flk| — Fp — Falk — 1], which leads to an incorrect 
value for u [k]. The correction lies on the tasking order alteration where the state vector is 
updated before the output equation, which means that x[k| = f(x[k — 1 ulk — 1]) is 
calculated and, then, Y[k| is computed. Note that, in order to achieve ulk — 1], the input force 
must also be delayed. This method is valid as long as the output y[k|] does not depend on 
ulk]|, so matrix D, must be zero, and it is given by: 


x|k| = Agxlk — 1) + Bqulk — 1), (11) 


Ylk] = Caxlk]. (12) 


4.2 Polynomial extrapolation method 


This method consists of using the previous values of a signal x(t) to predict its future 
value (WALLACE, 2005). In order to do so, an N-th order polynomial is calculated each 
time step by using and its previous values vlk — 1], ylk — 2], ..., y|k — (n — 1), where is the 
number of data points and it is related to the polynomial order according to n = 'N'+ 1. The 
value of y predicted time steps ahead is denoted y' = y[k + P] and is given by: 


v[k] 


y'[k) = [00 + Qu] (13) 








y[k — N] 


where each a term is assembled into vector a, which is given by: 


l Ú q () 

=— " nm A | 

g=[L PR em) Es E ( Ts) | (14) 
| +A(n-Dre (DT 


4.3 Taylor series method 
A novel method is proposed here based on the Taylor series, which approximates 
a function f(x) near a point à based on values of its derivatives evaluated at this point, as 


indicated below, where is the -th derivative of f(x). 


n 


(n) 
f(x) = >. fi = ar (15) 


n! 
i=0 





When applied to a real-time simulation, the objective is to predict the numerical 


Coleção desafios das engenharias: Engenharia mecânica 2 Capítulo 10 106 


subsystem output at a future time y (t + T), where t is the current time. By adopting a = t and 
x=t + T on Equation 15, the prediction is given by: 


eia (E 


| (16) 
IL: 


vC+T) = y(0 + y8P(OT + + 


By modeling the numerical subsystem on state-space form, one has access to u(b), 
yo (t) and vê (t) through the states and inputs of the system. It is known that y(t) = x,(t) 
and y? (t) = x, (t). The second derivative can be obtained from the input u and the states. 
Through Newton's Second Law, it is known that y? = x, = u/m — cx,/m. By substituting 
the previous expressions into Equation 16, the predictor's expression is given by: 


YC+D = (0 + xo + (SED) eo, (17) 


Im Mm 


where E(t) represents the truncation error. The discrete version of the predictor is given below: 








m m 2 + elkl. (18) 


vlk]= xlk|+ x5|kI(TGP) + (e pi 5 


Adopting C'=[|1 T;P-=c(T;P)Z/(2m)] and D' = (TsP)*/(2m), Equation 18 can be 


written in matrix form as: 


y'[k] = C'x[k] + D'ulk] + elk]. (18) 


4.4 Numerical simulations 


In order to compare the methods listed before, numerical Hil simulations were 
conducted by using a model of the spring as the virtual hardware. The effect of the sampling 
time was analyzed through two simulation setups: Setup 1 with T,> T, and Setup 2 with T, 
nd Ts Block diagrams were created in SIMULINK in order to simulate the HiL test with delay. 
The one used for T,> T, is shown in Figure 5a, while the one used for T,< T, is shown 
in Figure 5b. When adopting T,> T.., the effective delay becomes one time step, which is 
modeled by the 1/z block. For T,< T,, the delay was modeled as being proportional to 7. 
by a factor of d. The K term refers to the value of F (x = 0) from Equation 8. An uniformly 
distributed random wave signal with amplitude 4, is the external force acting on the mass. 
Physical parameters chosen for the simulations are m = 50kg, F,= 490.5 N, 4, = 100 N, w,. 
=rad/s. For the damping characteristic, it was used c = 600 Ns/m for Setup 1 and c = 200 
Ns/m for Setup 2. 
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Figure 5: Block diagrams used for numerical simulations for (a) Setup 1 and (b) Setup 2. 


The time delay used for the simulations was equal to the measured one. T. = 0.055 
was chosen for Setup 1 and T. = 2 ms for Setup 2. Figures 6a and 6b show the results for 


Setups 1 and 2, respectively, compared to the non-delayed system. 
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Figure 6: Simulation results for (a) Setup 1 and (b) Setup 2. 


The delay compensation methods were implemented on SIMULINK. A one time step 
delay was used because the method was derived considering To a In order to simulate 
Setup 2, the block “Unit Delay” is replaced by a general 'Delay' block. Figures 7a, 7b and 
7c show the simulation results for Setup 1 using the Tasking Order Sequence Alteration, 
Polynomial Extrapolation and Taylor Expansion Methods, respectively, and Figure 7d shows 
the difference between all the previous compensated signals compared and the correct one. 
lt can be seen that the Tasking Order Sequence Alteration Method provides, in simulation 
level, a perfect delay compensation. But it must be noted that it is valid when T.> T, and D 
= O. Thus, its use may not be suitable due to sampling period limitation, since it cannot be 
used for HiL testing of rapidly varying systems. For the Polynomial Extrapolation and Taylor 
Expansion Methods, it is possible to notice that the second one achieved a good result, 
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while the first did not. This fact indicates a limitation of the polynomial method related to 
sampling time, while the Taylor method provides better results. Finally, the results for Setup 
2 are given in terms of the difference between correct and compensated signals, as shown 
in Figure 8. It can be seen that the Tasking Order method presented a bad performance 
compared to the other methods, which are more suitable for situations where Ts a 


Compensated System 
— (COrTect signal 


Compensated System 
— (COrTect signal 


y [mm] 
y [mm| 


Compensated System 
— Cormect signal 


— Taylor Expansion 
- — = Tasking Order 
Polynomial Extrapolation 


MA 


Im 


0 2 4 6 8 10 
t [s] 
(d) 
Figure 7: Simulation results for Setup 1. (a) Tasking Order Sequence Alteration Method, (b) Polynomial 


Extrapolation Method, (c) Taylor Expansion Method, (d) Difference between correct and compensated 
signal. 
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Figure 8: Difference between correct and compensated signal for all methods for Setup 2. 


5 | EXPERIMENTAL RESULTS 


The methods presented in the previous section were implemented on the real 
hardware. By using the equipment listed in Experimental Setup section, the HiL test was 
conducted. A damping of c = 300 Ns/m and time steps of 0.1 s and 2 ms were chosen. 
For = 0.1 s, Figure 9a shows a comparison between compensated, uncompensated 
and correct signals. Due to noise in the signals, the compensation method is not able to 
provide a perfect test, but it can be seen that the response was improved by employing 
the method. Figure 9b shows the results for the Polynomial Extrapolation Method with Um 
= 2 ms. The measured signal tracked the correct one almost perfectly when compared to 
the uncompensated response. For this case, a second-order polynomial was used. Finally, 
Figure 9c shows the test result when using the Taylor Expansion Method, which clearly 


differs from the correct result. 
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Figure 9: Experimental results with delay compensation methods. (a) Tasking Order Sequence 
Alteration Method, (b) Polynomial Extrapolation Method, (c) Taylor Expansion Method. 
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The compensated systems did not follow the simulation because, in a real scenario, 
there is measurement noise and, also, time delay may not hold constant. The Tasking Order 
Sequence Alteration Method presented good result. However, it can only be used when large 
time steps are chosen, which may not be possible in certain applications since it influences 
the stability of discrete systems. In addition, large time steps may lead to abrupt changes, as 
can be seen in the results for t < 1 s from Figure 9a. The Polynomial Extrapolation Method 
was the most efficient, due to the proximity between correct and compensated signals. A small 
steady state error can still be perceived, which may have occurred due noise and variable 
time delay. Although it was not investigated in this paper, a limit for the TT, relationship 
may exist in order to maintain stability and accuracy, and it should be investigated in future 
works. Finally, the Taylor Expansion Method presented a very low performance compared to 
its excellent simulation behavior. This may have happened mainly due to the error estimation 
€ [X]), which is affected by noise in the force measurement, and also due to variable time 
delay. This method should be further improved in future works. Although not analyzed in this 
work, it is believed that the methods can be used for more complex systems including non- 
linearities and more degrees of freedom, since no approximation was done to derive such 
compensation techniques. Different test setups can also be used, such as different actuator 
and sensors, where the only difference would be regarding the amount of time delay and the 
dynamic characteristics of the sensors, which would have to be identified. 


61 CONCLUSIONS 


In this paper, the Hardware-in-the-Loop simulation technique was presented, studied 
and implemented in a mass-spring-damper system, where the hardware chosen was an 
already parametrized spring, in order to validate the results. A time delay was identified and 
its effects on the closed-loop behavior were analyzed. It was shown the sampling rate of the 
discrete HiL system influences stability and the relation between delay and sampling time must 
be analyzed prior to the test. Two delay compensation methods, Polynomial Extrapolation 
Method and Tasking Order Sequence Alteration Method, were studied and implemented. A 
new technique, based on Taylor series, was presented and compared. In terms of numerical 
simulations of the HiL test, it was seen that the proposed method overcame the other two. In 
terms of experimental tests, the Polynomial Extrapolation Method showed better performance. 
New ways for estimating the truncation error and also studies of the measurement noise 
influence on the methods should be considered for future works. Since no approximation 
regarding the system's dynamics was made, the techniques can be employed for more 
complex systems and this should be tested in future works. 
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RESUMO: Como resultado da disputa cada dia 
mais complicada no mercado e na procura de 
maiores lucros, o foco de qualquer companhia 
tem sido na redução de retrabalhos e perdas, 
sendo assim, com a evolução constante dos 
processos e sistemas se tornou uma prioridade 
dos gestores modernos atuarem conforme a 
necessidade. Com tudo, identificar, analisar e 
sanar este problema vem a ser um processo 
complexo e que envolve várias equipes de 
trabalho. A empresa de fabricação de rolhas 
metálicas localizada no polo industrial de 
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INDUSTRIAL DE MANAUS 


Manaus, propõe-se melhorias nos indicadores de 
manutenção e para isso orientou a aplicação da 
metodologia dos seis sigmas como ferramenta 
para o gerenciamento da manutenção. O método 
DMAIC (Definir, Medir, Analisar, Melhorar e 
Controlar), que é parte integrante da metodologia 
Seis Sigma, fornece um meio ordenado a ser 
seguido por meio de ferramentas precisas de 
análises e soluções de problemas, a utilização 
do DMAIC traz uma confiabilidade e precisão na 
execução dos projetos realizados nas empresas. 
Os resultados mostrados com a implementação 
na empresa, permitem visualizar os principais 
problemas que limitam a produção, analisar 
estes e desenvolver soluções para a melhoria 
do processo. A metodologia mostrou ser válida 
no trabalho ao se obter avanços em termos de 
eficiência do processo e diminuição das perdas 
por paradas não planejadas em um 8 %. 
PALAVRAS-CHAVE: Metodologia seis sigma, 
perdas metálicas, paradas não planejadas, 
confiabilidade. 


SIX SIGMA METHODOLOGY: 
APPLICATION IN À COMPANY OF METAL 
STOPPERS OF THE INDUSTRIAL POLE 
OF MANAUS 


ABSTRACT: As a result of the increasingly 
complicated dispute inthe market andinthe search 
for greater profits, the focus of any company has 
been on reducing rework and losses, so, with the 
constant evolution of processes and systems, it 
has become a priority for modern managers to 
act. as required. However, identifying, analyzing, 
and solving this problem is a complex process that 
involves several work teams. The metal stopper 
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manufacturing company located in the industrial hub of Manaus, proposes improvements in 
maintenance indicators and, for this, guided the application of the six sigma methodology as 
a tool for maintenance management. The DMAIC (Define, Measure, Analyze, Improve and 
Control) method, which is an integral part of the Six Sigma methodology, provides an orderly 
way to be followed through accurate analysis and problem solving tools, the use of DMAIC 
brings reliability and precision in the execution of projects carried out in the companies. The 
results shown with the implementation in the company, allow to visualize the main problems 
that limit production, analyze these, and develop solutions to improve the process. 

The methodology proved to be valid at work by obtaining advances in terms of process 
efficiency and reducing losses from unplanned downtime by 8%. 

KEYWORDS: Six sigma methodology, metallic losses, unplanned downtime, reliability. 


11 INTRODUÇÃO 


Planejar o impacto de paradas não planejadas em um processo é de grande importância 
para o setor de Planejamento de manutenção. Os envolvidos com esse processo avaliam as 
variáveis para obter os resultados desejados pela companhia. Atualmente, é importante para 
qualquer empresa, simplificar os seus processos tendo uma queda significativa nas perdas 
e desperdícios. Com isso, autuar com processos de produção mais eficientes tem sido o 
grande objetivo das pequenas e grandes empresas. 

As paradas não planejadas em quase todos os casos têm perdas bem grandes. 
Em desafios como este, é normal existir problemas entre (quantidade de problemas x 
quantidade de pessoas para resolver problemas), sendo assim, se ter um processo de 
soluções ágeis, se tem por análises rápidas que procuram atacar o efeito que o problema faz 
e não a causa do mesmo, em grande parte acaba gerando custos de retrabalhos e com isso 
se ter vários reparos nos processos. Para Hardy, DL, Kundu, S. e Latif, M. (2021), encontrar 
uma maneira de zerar perda sempre será melhor que demitir colaboradores, além do mais, 
ter o desafio de baixar o custo de uma companhia é uma tarefa bem satisfatória 

Por vez a resolução de problemas a curto prazo é necessária diante da necessidade. 
Neste caso, muitas estratégias são criadas pelas companhias em busca da melhor aceitação 
de um sistema de produção com a sua realidade, como por exemplo a metodologia SIX 
SIGMA apresentada por Dias (2011). 

Para Dias o método SIX SIGMA, tem grande foco em melhoria continua tendo uma 
evolução da gestão da qualidade, levando-a por várias iniciativas a fim de melhorar os níveis 
de eficiência e eficácia das atividades. Para Schroeder et al. (2007) o método das seis 
sigmas é definido como “uma meso-estrutura paralela, organizada para reduzir a variação 
de processos utilizando-se de especialistas em melhoria, um método estruturado e métricas 
de desempenho com a meta de atingir objetivos estratégicos”. Terner (2008) diz, ressalta que 
o mercado já não tem mais espaço para ineficiência, para ele a concorrência já ultrapassa 
o critério de preço, o foco agora está voltado para o valor ligado a padrões de qualidade, 


eficácia dos serviços e preço justo. 
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O método DMAIC (Definir, Medir, Analisar, Melhorar e Controlar), que é parte 
integrante da metodologia Seis Sigma, fornece um meio ordenado a ser seguido por meio de 
ferramentas precisas de análises e soluções de problemas, a utilização do DMAIC traz uma 
confiabilidade e precisão na execução dos projetos realizados nas empresas. 

Desta forma, é muito importante as empresas começarem a desenvolver projetos que 
utilizam metodologias e que tenham ferramentas para análise de problemas, para que o 
envolvido obtenha dados e possam estruturar o processo de tomada de decisões para poder 


ter uma opção que seja assertiva para a companhia. 


21 REFRENCIAL TEÓRICO 


2.1 Método das seis sigma 


O método das seis sigma (SIX SIGMA) se originou na Motorola no ano de 1987 com foco 
na redução de perdas de qualidade dos produtos e tendo um maior retorno na lucratividade. 
Segundo Werkema (2004), conceituo o Six Sigma como: “medir o desempenho da linha atual 
calculando quantos sigmas existem até que ocorra a insatisfação de um cliente, ou seja, que 
um defeito é identificado pelo cliente”. Um processo será considerado uma metodologia Six 
Sigma quando não gerar mais de 3,4 dpmo (defeitos por milhão de oportunidades). 

Corrêa e Corrêa (2009) citam que essa variabilidade está relacionada aos resultados 
gerados pelos processos, e que alto índice dessa variabilidade está ligado a má qualidade, 
altos custos e insatisfação do cliente. 


2.2 Metodologia DMAIC 


Atualmente o método mais utilizado para assessorar esse desenvolvimento do Six 
Sigma é o DMAIC, selecionada por Cleto (2011) para conduzir a gestão de projetos em 
uma empresa automobilística, uma vez que para ele esse método possibilita a organização 
apropriada de implementação, desenvolvimento e conclusão da maior parte dos projetos Six 
Sigma. 

A metodologia DMAIC é mundialmente conhecida como evolução do programa Six 
Sigma pois sua abordagem principal é focada em melhorias de processos. Para Schroeder 
(2007), as etapas da metodologia DMAIC classificam os desvios e apontam oportunidades 
para melhorias através dos dados coletados e analisados. 

De acordo com Lucas (2002), a base do DMAIC é o ciclo de melhorias PDCA (Plan, 
Do, Check e Act), sendo que ambos possuem características semelhantes, tais como escolha 
dos pontos a serem melhorados, análise, planos de ação, verificação/controle. 

A utilização do método de análise e solução de problemas DMAIC incorporado no 
gerenciamento de projetos traz uma sistemática de condução das etapas e ainda proporciona 
a utilização de ferramentas estatísticas para identificar um problema ou variação de um 


processo. para melhorar cada vez mais os produtos e com isso atender as exigências dos 
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clientes internos e externos. 

Lucas (2002) coloca o DMAIC como uma estratégia de ruptura (breakthrough 
strategy) que fornece ao processo de melhoria uma forma organizada de ação, com técnicas 
específicas para cada uma de suas etapas. 

O ciclo de melhoria do DMAIC é utilizado também como um método para gerenciar os 
projetos uma vez que cada uma das fases apresenta entregas bem claras, tanto com relação 
ao formato da análise quanto em relação aos métodos a serem utilizados nessa análise. Esse 
entendimento de entregas e de gerenciamento de fases é demonstrado por George (2002) e 


Pande; Neuman e Cavanagh (2001). 


31 METODOLOGIA 


O presente trabalho utilizará a metodologia do Estudo de Caso, pois, a obtenção de 
dados se fará com o uso de documentação indireta e direta in loco, por meio de pesquisas 
qualitativas no qual faz aplicação de modelos teóricos à situação prática da vivência 
organizacional mediante a observação sistemática e levantamentos de dados que irão 
futuramente servir para medir quantitativamente o caso estudado. (LAKATOS e MARCONI, 
20083). 

As etapas do DMAIC, Figura 1, são definidas da seguinte forma: 


* | Definir (Define): se refere a finalidade do projeto e o as necessidades para sua 
resolução. Tem como objetivo estabelecer um desenho detalhado do processo 
com seus pontos críticos e relevantes para o problema. As ferramentas geral- 
mente utilizadas nesta etapa são brainstorming, histograma e fluxograma para 
auxiliar na identificação dos riscos, comportamento e percepção operacional do 
problema. 


*- Medir (Mensure): tem como finalidade determinar, de forma quantitativa, as va- 
riáveis associadas ao problema. Nesta fase é feito uma coleta mais detalhada 
dos dados do problema e estratificá-los é preciso para entender como ocorrem. 
Para realizar as atividades, podem ser utilizadas ferramentas como Folha de 
Verificação para fazer a coleta de dados e gráfico de Pareto para estratificação. 
Caso haja problemas mais complexos, pode ser utilizado o software Minitab. 
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Figura 1. Sequência do método DEMAIC 


Fonte: Autoria própria. 


* Analisar (Analyse): abrange o reconhecimento das causas potenciais dos pro- 
blemas que estão impactando no processo. Nesta fase, auxiliam na análise dos 
dados, ferramentas como Diagrama de Ishikawa e análise de 5 PQ's. 


* | Implementar (Improve): consiste na implementação das melhorias identifica- 
das para solucionar a causa dos problemas. Nem todas as melhorias e ideias 
propostas podem ou devem ser implementadas. Para priorizar as ações a ma- 
triz Impacto x Esforço é comumente utilizada. 


* | Controlar (Control): Monitorar o desenvolvimento da estratégia de execução 
para que o grupo não se perca e que venha estabelecer critérios de controle 
(Check-Lists, Estatísticas, além disso, podendo criar padrões, instruções e dia- 
gramas de trabalho para que independentemente haja mudanças na mão de 
obra. 


41 RESULTADOS E DISCUSSÃO 


4.1 Generalidades 


O processo de fabricação de rolhas, apesar de “simples”, possuem certos detalhes, 
que se não observados com cautela podem acarretar problemas ainda maiores como por 
exemplo problema de vedação nas garrafas, ferrugem na rolha podendo misturar com o 
líquido. As variáveis que impactam em um processo de fabricação de rolhas vão desde a 
qualidade da folha do aço (AÇO T4) até qualificação técnica-operacional sobre como operar 
os equipamentos. 

As rolhas metálicas desse estudo são compostas por 2 partes, sendo elas, a folha com 
formato de rolha (tampa) e o vedante, que tem aproximadamente 0,240g + 0,0159, porém 
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o corpo da rolha recebe uma camada de verniz para proteção contra oxidação. A partir as 
rolhas metálicas possuem duas variações que são as seguintes: TWIST-OFF, PRY-OFF que 
são para determinados tipos de cervejas. As linhas de produção de embalagens de rolhas 
metálicas operam por 24h/dia por cerca de 350 dias ao ano. O tempo de setup (mudança de 
produto) e tem em média 1 hora pois precisa ser feito a limpeza de toda a linha e verificação 
se não há nenhum resíduo. 

O planejamento de produção é calculado e planejado conforme a necessidade 
e disponibilidade da linha, portanto, qualquer imprevisto seja ele relacionado a máquina, 
material, método e mão de obra, comprometendo as demandas necessárias. Fazendo uma 
análise inicial do problema, em determinado mês, foi notado que mais de 25% do tempo 
disponível para produção estava ocupado por paradas imprevistas conforme é mostrado na 
Figura 2. 


Paradas Planejadas 
14% 


Paradas não 
Planejadas 
26% 





Figura 2 - Disposição de paradas não planejadas em linhas de produção de rolhas metálicas 


Fonte: CENG 2019- AROSUCOS 


Neste caso, foi escolhido a realização de um projeto pelo método DMAIC- White Belt 
para se ter uma análise mais concreta do problema devido ao fato de ser um grande ponto 
negativo que afeta a produção de modo que gera retrabalho, custo e gasto de mão de obra, 
com isso, devido a facilidade de e usar as ferramentas e os resultados terem valores bastante 
significativos o método DMAIC foi assegurado como forma de sanar esse problema. 


4.2 Definir 


As paradas não planejadas têm total influencia na eficiência do processo e no valor 
gerado pelas perdas. As perdas podem ser facilmente calculadas pelo tempo de produção 
perdido e também pela matéria prima (aço), perdida nas paradas de linhas inesperadas, 
com isso chama-se na produção em estudo por “perda de GLY”. Pode ser citado também as 
perdas com os gastos de mão de obra que pode sem calculas considerando a hora/lhomem 
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sem produzir um produto de boa qualidade. 

Tendo base na Figura 3, as paradas imprevistas representam 25,6 % da produção. 
Em um dia de uma grande demanda de produção, esse valor pode representar cerca de 6 
horas de produção perdida, considerando apenas a perda que acontece no equipamento, 
em que ele não está em funcionamento por possíveis ajustes ou então por uma parada 
não planejada, durante um mês a empresa tem um custo negativo de aproximadamente R$ 
700.000 mil reais. A equipe necessária para a solução deste problema foi composta por um 
Técnico Mecânico do processo e dois supervisores sendo eles: Supervisor de Manutenção 
e Supervisor de Perda de GLY (perda de produção), com a orientação deles para ter uma 
continua melhoria do processo. 

O propósito do projeto limitou-se a estudar quatro grupos de paradas não planejadas: 
1. Pequenas Paradas 2. Espera no processo (SET-UP) 3. Quebras 4. Falha de processo. 

Como as pequenas paradas teve uma maior tendência representativa, o proposito 
estreitou apenas para o estudo deste grupo. Para representar as paradas não planejada, ouvi 
a divisão sem subcategorias, como mostra a Figura 3. 

Para ter noção do problema foram coletados dados do Sistema SAP ERP sobre os 
gatilhos de paradas, no período de janeiro a maio de 2021, com isso os dados selecionados 
foram expostos em um histograma, para se ter os valores referentes a % das paradas não 
planejadas sendo calculado através da equação 1: 

% PNP = ((Total de horas de Paradas não planejadas) / Tempo total disponível- 
Paradas Planejadas)> *100 





Paradas nao 
planejadas 


[eee peer) 





Pequenas Esperas de None Falhas do 
paradas processo processo 


Enroscamentos 


Regulagens 


E rompimento do 
fio de cobre 


Figura 3. Paradas não planejadas. Estratificação dos problemas (Pequenas Paradas) 





Fonte: Autoria própria. 
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A partir disto, pode se notar na Figura 4, o histórico do problema que vem ser 
destacado a média do percentual (%) de paradas não planejadas no período analisado foi de 
aproximadamente 7,6% e o melhor valor encontrado foi de aproximadamente 6,4%. 


8,08% 











le Indicador 1 Media lu Benchmark 


Figura 4- Histórico do problema 


Fonte: CENG 2019- AROSUCOS. 


O objetivo deste projeto, alinhado com a liderança da empresa, foi reduzir mais de 46% 
do problema que vem ser a diferença entre a média e o menor valor, obtendo um percentual 
mensal de paradas de 7%. As fases do DMAIC, foram sintetizadas em um cronograma para 
se ter o planejamento das atividades e entrega no prazo proposto, conforme a Figura 5. 
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Figura 5- Cronograma das Etapas do DMAIC 
Fonte: CENG 2019- AROSUCOS. 


A efetivação para se iniciar um projeto DMAIC é feita através de um “contrato”, com 


as informações do projeto como: equipe, dados, meta, prazos e assinado por todas as 


partes que nele estarão envolvidas. 
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4.3 Medir 


Como os dados obtidos não demonstram confiabilidade, devido a distribuição dos 
dados serem feitos pelos próprios operadores das linhas onde o problema ocorre, foi verificado 
também alta dispersão dos dados, com isso, foi escolhido pelo uso da folha de aceite. Quando 
se tem dados não confiáveis isso se torna um problema para toda tomada de decisão, pois se 
torna suscetível a erros e possíveis retrabalhos. Durante um determinado período proposto, 
foi feita a coleta de dados manual das paras ocorridas nas linhas de produção, o Quadro 1 


mostra um modelo da folha de verificação utilizada pelo método. 


Quando medir? | Como coletar? | Porque coletar? | Responsável 
pela coleta? 


Quadro 1- Modelo de folha de Verificação 


Fonte: CENG 2019- AROSUCOS. 





Com a obtenção de dados confiáveis, o próximo passo seria estratificá-los. Pois a 
estratificação é o agrupamento de todos os dados em diferentes fatores. Como, máquina, 
lote, turno, tipo de parada e material. Com isso, o principal objetivo da estratificação e indicar 
os potenciais problemas que são prioridade ou não. Basicamente, serão realizadas análises 
dos dados das máquinas e a partir disso, foi feita a estratificação por tipo de parada ou 
ocorrência. 

As estratificações foram feitas com a utilização do Gráfico de Pareto e o princípio 
80/20, em que você tem 80% dos problemas são gerados por 20% das causas. Sendo 
assim, tanto para estratificação por problemas, quando na do tipo de parada, percebeu-se 
que foi 80% dos problemas potenciais, com isso a Figura 6 mostra o Gráfico de Pareto 
da estratificação do problema geral. Na Figura 7 foi possível perceber que os problemas 
que mais impactam no período, seguindo o conceito citado acima 80/20, foram: Amassando 
chapa, Descarte de rolhas pelo CVS, Mesa de esquadro desalinhada e Má formação. Com 
todas os problemas identificados, basicamente se resume a duas máquinas que vem a ser a 
Prensa (PTC) e o LINER. Com isso, a perda estimada por hora na máquina estar perto das 


198 horas durante o período de análise. 
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Mapeamento Geral dos Problemas 
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Figura 6- Gráfico de Pareto da estratificação do problema geral 


Fonte: Autoria própria 


Se estratificar-se o problema da prensa (PTC) em prioridades obtiverão se os 
seguintes valores conforme a Figura 7. 

Na máquina LINER, o período analisado e as paradas não planejadas somam 
aproximadamente 140 horas, sendo os potenciais problemas, descarte de rolhas pelo CVS, 
Má formação no vedante, dosagem do vedante na rolha, 1. Estrela de transferência, faca 
desgastada, Defeito na rádio frequência e punção formador. Abaixo na Figura 8 é ilustrado 
a Estratificação do problema na Máquina LINER. 


Principais problemas na Prensa (PTC) 
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Figura 7. Principais problemas na prensa PTC. 


Fonte: Autoria própria. 


A etapa medir é concluída com o levantamento das possíveis causas. Os problemas 
que podem ocorrer com o passar do tempo, em todos e qualquer projeto, deve ser feito com 
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Oo maior número de pessoas que estejam envolvidos com os problemas citados acima. 


Principais problemas no LINER 
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Figura 8- Principais Problemas no LINER. 


Fonte: Autoria própria. 


4.4 Analisar 


Nesta etapa, são estudadas e analisadas cada uma das causas potenciais levantadas 
pela etapa anterior. No brainstorming realizado no final da etapa anterior, foram levantadas 
ações imediatas para serem implantadas, mais conhecido como o plano de ação VER E 
AGIR. Para analisar, os demais problemas, utilizou-se a ferramenta Diagrama de Ishikawa, 
com isso distribuímos as possíveis causas nos 6Ms, sendo: Máquina, Material, Mao de 
Obra, Método, Meio Ambiente e Medida. Após ser feito isso, todas os possíveis problemas 
foram analisados atrás da ferramenta 5 Por Quês, para se ter a identificação das causas 


raízes. 


4.5 Implementar 


As soluções propostas devem ser implementadas para garantir três entregas 
fundamentais, sendo elas, padronização, resultado e conhecimento. O resultado está 
atrelado ao alcance da meta, a padronização deve garantir que não haja erros nas execuções 
das atividades, sejam elas procedimentos técnicos ou operacionais. Para certificar-se que os 
padrões estão sendo cumpridos à risca, pode-se fazer uso das ferramentas a prova de erro 
ou Poka-Yoke. 

A entrega do conhecimento está relacionada ao entendimento dos operadores e 
técnicos sobre como executar as determinadas tarefas em suas funções. Quando houver 
problemas que dependem da habilidade técnica ou operacional, a entrega de conhecimento 
se torna a mais importante para se chegar a grandes resultados. Para se transmitir 


conhecimento, se tem diversas formas de fazer o mesmo, com treinamentos em salas de aula, 
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treinamentos técnicos com fornecedores, reuni-os entre turnos para tratar dos problemas da 
empresa. Algumas vezes, como foi o caso desse estudo, devido ao número e a complexidade 
das ações foi preciso criar uma matriz de prioridade, como mostra a Figura 9. 







Prioritários (Plano de Ação 5W2H) 
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Figura 9. Matriz de prioridade. 
Fonte: CENG 2019- AROSUCOS. 


Na pesquisa foram mapeadas 57 ações e todas foram priorizadas através da matriz 
de prioridade. Com auxílio do plano de ação 5W2H, do total de ações, 51 foram levadas para 
o plano, 3 descartadas e 3 foram classificadas como complexas para serem tratadas em um 
futuro. 


4.6 Controlar 


Nesta etapa, a parte de controlar teve uma duração de quase 3 meses e seu principal 
foco foi acompanhar se os resultados obtidos foram e estão sendo contínuos. Para tomar 
nota, os resultados planejados foram obtidos, foi utilizado novamente o gráfico sequencial 
para visualizar o percentual de paradas não planejadas obtidas ao longo do tempo, conforme 
está na Figura 10. 
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Figura 10- Gráfico Antes contra Depois da Ações implementadas. 


Fonte: CENG 2019- AROSUCOS. 


Nesta etapa, além de se ter o controle propriamente dito dos resultados obtidos, foram 
atualizados os procedimentos operacionais e teve um aumento da eliminação de ajustes 
desnecessários no processo. 

Por fim, com os resultados obtidos teve uma redução bem significativa de 
aproximadamente 22 horas/mês de paradas imprevistas. Para se ter um processo contínuo 
desses resultados outras ações foram implementadas, como uma maior capacitação dos 
colaboradores sobre os procedimentos e também se tem um sistema para se ter acesso aos 


pontos de ajustes e regulagens necessários das máquinas para possíveis erros. 


5 | CONCLUSÃO 


Com base nos resultados desta pesquisa destacam-se as seguintes conclusões: 

A redução média e de 20% ao mês do problema, com isso teve um impacto de 8% 
de queda em paradas não planejadas das linhas em geral saindo de aproximadamente 25% 
para 1/%. 

As conclusões desse artigo e os resultados obtidos dentro da empresa afirmam que as 
citações da literatura de falam sobre a metodologia DMAIC se obteve um resultado positivo 


através do mesmo. 
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ABSTRACT: In this research work a proposal 
is made for a device to measure vital signs for 
elderly people, mainly in relation to their heart 
rate and respiration. The device also includes a 
GPS (Global Positioning System) system and a 
help button, which can be activated if the user 
falls or needs special help, remotely monitoring 
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user activity. The methodology used to carry 
out the analysis and selection of the model is 
based on the Innovation Cells program of the 
Confederation of Industrial Chambers of the 
United Mexican States (CONCAMIN). As a final 
product, a concept with local market innovation 
characteristics is presented. 

KEYWORDS: Electronic wristband, vital sign 
measurement, remote system, innovation cells. 


PROPOSTA DE UM DISPOSITIVO DE 
PULSEIRA ELETRONICA PARA MEDIÇÃO 
DE SINAIS VITAIS 


RESUMO: Neste trabalho de investigação é 
feita uma proposta de um dispositivo para medir 
sinais vitais de pessoas idosas, principalmente 
em relação à sua frequência cardíaca e 
respiração. O aparelho ainda tem um sistema 
GPS (sistema de posicionamiento global) e 
um botão de ajuda, que pode ser acionado 
caso que o usuário caia ou necessite de ajuda 
especial, monitorando remotamente a atividade 
do usuário. A metodologia utilizada para realizar 
a análise e seleção do modelo é baseada no 
programa Células de Inovação da Confederação 
das Câmaras Industriais dos Estados Unidos 
Mexicanos (CONCAMIN). Como produto final, é 
apresentado um conceito com características de 
inovação do mercado local. 
PALAVRAS-CHAVE: Pulseira eletrônica, 
medição de sinais vitais, sistema remoto, células 
de inovação. 
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11 INTRODUCTION 


The total population in Mexico amountis to 123.6 million people, of which 17,246 million 
are between 50 and 64 years old and 9,735 million are 64 years old or older (NATIONAL 
INSTITUTE OF STATISTICS AND GEOGRAPHY, 2018). In 2050 they will represent about 
30 percent of the Mexican population, for that year our life expectancy will be 86 years 
(HERNÁNDEZ, 2018). Globally, by 2050, the number of people aged 60 and over is expected 
to increase from 600 million to almost 2 billion, and the percentage of people aged 60 and 
over is projected to double from 10% to 21% (ALVARADO GARCÍA; SALAZAR MAYA, 2014). 
This increase will be greater and faster in developing countries, where the elderly population 
is expected to multiply by four in the next 50 years (ONU, 2002). The figures are considerably 
alarming, and they must be present in public policies, socioeconomic development and the 
potential market in order to minimize the needs of the elderly, who must adapt to them, to 
maximize the health and functional capacity of older people, as well as their social participation 
and safety (OMS, 2019). 

There are multiple definitions of aging, according to several authors, they agree 
on treating it as a dynamic, multifactorial process inherent to all human beings (GÔMEZ; 
CURSIO, 2002). Aging can be defined as the sum of all the changes that occur in humans 
over time and that lead to functional deterioration and death, and includes those changes that 
occur in any cell or organic system as a function of the time, regardless of any external or 
pathological influence such as disease (BAZO, 1998). 

Among the signs that determine the aging of people, we have (INM, 2021): 

Physical age: physical and biological changes that occur at different rates, which 
depend on sex, place of residence, economy, culture, food, type of activities carried out and 
emotions. 

Psychological age: changes in emotions, feelings, thoughts and the meaning that old 
age has for each person. Additionally, there are changes in psychological processes, such as 
memory or learning. 

Social age: related to the meanings of old age, different for each human group, 
according to their history, their culture and their social organization. 

Regarding the state of the art of electronic devices for monitoring vital signs, there 
are multiple devices developed or in design stages. Gutiérrez (2016), develops a device that 
consists of a continuous monitoring system of vital signs and environmental variables, with 
non-invasive sensors and wireless data transmission, to determine the current health status 
of the user and the conditions of their environment, in order to provide medical assistance if 
required. (REYES; ASTO, 2018), designed a monitoring system with an interface for doctor- 
patient interaction, which allows monitoring vital signs of patients in emergencies treated 
in the working conditions of the Cia. Salvadora Trujillo Fire Department No. 26. A portable 


device was built with sensors that measure vital signs (heart rate, oxygen saturation, blood 
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pressure and body temperature) of patients. Velartus, (2018), make use of wearable devices 
for prevention, well-being and the current state of available medical clouds. 

Biomedical technology aims to improve the efficiency and quality of the provision 
of health services, through a series of resources, especially of the technological type, that 
provide greater economic benefits and high effectiveness of services (RUIZ, 2005). On the 
other hand, biometric technology is responsible for taking standard measurements of living 
beings or biological processes, and as a result, data can be obtained through the use of 
invasive systems such as blood samples, or non-invasive, such as sensing systems. 

In this research work, a proposal is made for a device to measure vital signs for elderly 
people, mainly in relation to their heart rate and respiration. The device also includes a GPS 
(Global Positioning System) system and a help button, which can be activated if the user 
falls or needs special help, remotely monitoring user activity. The methodology used to carry 
out the analysis and selection of the model is based on the Innovation Cells program of the 
Confederation of Industrial Chambers of the United Mexican States (CONCAMIN). As a final 


product, a concept with local market innovation characteristics is presented. 


21 METHOD DESCRIPTION 


For the proposal of the non-invasive glucose meter, the Innovation Cells methodology 
was applied, coordinated by CONCAMIN, through a structured process that allows putting into 
practice skills, methodologies and innovation tools in solving real problems of companies or 
entrepreneurs. 

The Innovation Cells program is shown in Figure 1 and consists of ten methodological 
stages (HERRERA, 2018). 

Stage 1. Define: the objective of this stage is to attract innovation projecis from 
companies and support them in clearly defining an innovation challenge and the profiles 
required for its development. 

Stage 2. Recruit: attract university students with desirable profiles based on the 
definition of the challenge. 

Stage 3. Kick off: program start; coexistence between student-company where work 
roles are discussed and known. 

Stage 4. Explore: Teams discover, design, and implement exploration tools and 
techniques to generate valuable information. 

Stage 5. Devise: teams discover, design and put into practice creativity tools and 
techniques to generate a large number of solution proposals to their challenge. 

Stage 6. Feedback: event in which teams present their solution proposals to a large 
number of people to provide feedback and select the best solutions. 

Stage 7. Prototyping: teams discover, design and put into practice prototyping tools 


and techniques to materialize and be able to show their solutions. 
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Stage 8. Validate: teams discover, design and put into practice validation and 
experimentation tools and techniques to assess the relevance of their solutions and improve 
them. 

Stage 9. Final presentation: event in which the teams present the final results of the 
process and show their prototypes or concepis. 

Stage 10. Evaluate: evaluation of the performance at the team level and at the individual 
level of all the participants. Identification of the next steps of the project. 


4. Define E» E> 3. Kick-OF 
6. Feed back € 4 4. Explore 


v 


7. Prototype B 8. Validate o E> 10. Evaluate 


Figure 1. Stages of the methodology used in this research. 





31 RESULTS 


Based on the proposed methodology, the results were as follows. 

Stage 1. The company Aarson Productos Hospitalarios de Hermosillo Sonora, 
launches the challenge of proposing a biomedical device with significant impact on the 
national and international market. 

Stage 2. An interdisciplinary group of university students was formed, made up of two 
biomedical engineering students, one mechatronics engineering student and one business 
management engineering student, accompanied by a professor from the metal-mechanics 
area with a profile in mechatronics. 

Stage 3. A visit was made to the company Aarson Hospital Products by the students 
and advisor, where the general aspects of the program, the scope and limitations of the 
proposed challenge and a global analysis of the problem were discussed. 

Stage 4. The basic information collected in this stage was as follows. 

4.1. Innovation opportunities.For the development of this section, the research 
focused on facts, through information obtained through interviews, observation or any 
research method, among which the following stand out: 1) the device will have to be practical 


and lightweight, 2) the first source of help is the youngest child or the partner and 3) when 
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being in a risky situation, the user chooses to lie down in the closest place. The preliminary 
findings from the facts are shown in Table 1. 


Topics to continue 


Main conclusions Surprises Opportunities 





exploring 


Most users have more 


than two children. Most The user usually The device can be 
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Table 1. Most important findings found in the investigation. 


4.2. Description of the trend. Table 1 shows the most relevant aspects of the trend, 
figure 2 represents a schematization. Figure 3 shows the trend and impact of the technology, 
according to the perception of the local market. 
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Illustrations: What will the product look like”? 


Figure 2. Electronic bracelet to measure vital signs. 


Description of the trend Impact on the market 
What exactly is it or does it do? Electronic | What is it used for? Currently there are 
bracelet that measures the pulse and heart rate, | already several ” wellness “ applications 
also provides the location via GPS of the wearer | and devices that allow monitoring vital 
If the gyroscopes detect an abnormal movement | signs through a device linked to a certain 
to people designated as sources of help that; If | app on a smartphone, as well as Apple 
the person is transferred in an emergency, the | Watch, which is more directly, without 
device will have the database of the patient, his external devices. 
medical history and family. 
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Potential applications or possible innovations. How can it work for the project? 

* Strategic arrangement of the sensors to achieve a practical, simple and functional design. 
* Specify the location of the sensors and reader in the device for a correct measurement. 

* Geolocation of the person in case of loss. 

* Faster response speed in an emergency. 

* Inclusion of medical and contact information in the event of an accident. 

* Provide a grip that is not prone to slack or loosening. 


Evaluation and consequences Where do you position this trend? 


High 


Technological Complexity 





Figure 3. Product trend. 


Table 1. Description of the product trend. 


4.3. Description and specification of the concept. Table 2 shows the most relevant 
aspecis regarding the description of the concept. Important points are established such as the 
sector of focus of the research, the details of the concept, USP (Unique Selling Proposition) 
and benefits for the user, USP and benefits for the client, promotion, distribution and allies, 
the impact for the business and the value statement. Table 3 shows the specifications of 
the concept. 









Target segment 1070) afe7=]o) Mo [FT] E 


Who? Elderly people with conditions that promote What does it contain or what is it about ”?: 
immobility of the user, memory problems and in Bracelet with pressure monitoring, location 
need of constant care and user database functions. Distribution in 
exclusive pharmacies (such as YZA, Benavides, 
etc.) and online. It will be packed in a box with 
user manual and security seals 








USP and user benefits USP and customer benefits 
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Why does the user need and will buy the 
product” It satisfies the need for constant 
monitoring and speeds up response time in the 
event of an unforeseen event. It will be focused 
on the elderly and people with Alzheimer's 


or a trusted person (Source of help) 








Promotion, distribution and allies 


Mobile advertising, infomercials will be given 

in the news and on social networks. Alliances 
with pharmaceutical companies. Distribution to 
pharmacological chains and home deliveries by 
shipping companies. 


disease and its monitoring will be done by a tutor 


Why will the suppliers or allies support or sell 
this product? The potential market is very 
extensive and there is no direct competition 
The product will provide comfort to the user's 
relatives 











Business impact 


Creation of adjacent products for medical 
institutions. A very wide market segment 
with very litile competition that is aimed at a 
completely different type of clients. 








Statement or value proposition: SECURITY AND CONTROL AT YOUR FINGERTIPS! 





Table 2. Description of the concept. 


Concept criteria 


User or client requirements 






































Timely, immediate 
4 Product A.- Timely synchronization of the 10 8 12 
reliability synchronization gyroscope is needed to alert in 
the event of an emergency. 
B.- Accurate The margin of error of the 
e 10 9 11 
readings readings must be minimal 
Provide the user's location to 
2 Functionality C.- Geolocation the person assigned as the 10 7 13 
Of the device first source of help. 
Include the user's medical 
a ando history, as wellasthe contact 9 5 13 
details of the first aid source 
The product must be optimal 
for daily use and of a 
3 a a considerable size (no larger 8 8 8 
y p y than the palm of the hand) and 
weighing less than 150gr 
EL piscraion The product must use neutral 4 7 9 
colors and be ergonomic. 

















Table 3. Concept specifications. 


4.4. Competition map. Figure 4 shows the main competitors in the product market, as 


well as the qualitative and quantitative characteristics of each of them. With this image, we 


can verify which are the strengths, opportunities, weaknesses and possible threats of the 


product with respect to the competition. 
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Ml — Wellness 
2 
—— Apple watch 
, 


Conventional 





Complexity Discretion Accessibility Practicality Precision Durability 


Figure 4. Competition map. 


4.5. Map of innovation opportunities. In figure 5, some technological and innovation 


characteristics of the system are shown, all of them represented in an opportunity map. 


Innovation opportunities: 
r oB A.- Timely synchronization. 
B.- Accurate readings. 
C.- Geolocation 
D.- Medical and contact information 
' EÉ.- Size and portability 


F.- Discretion 


Satisfaction 





Importance 


Figure 5. Map of innovation opportunities. 


41 CONCLUSIONS 


Multidisciplinary teams of diverse engineering, in cooperation with school academy and 
industry, can solve problems of a great diversity of branches, for example, those belonging to 
the health branch, based on well-structured, proven and well-founded methodologies. 

In the biomedical sector, the areas of opportunity for innovation and patenting of new 
products open up to a wide panorama of possibilities, with the characteristic that it is possible to 
work in synergy with engineers from various disciplines, who can be biomedical, mechatronic 
and engineering in business management (to name a few), to result in value propositions for 
the industry, which can later be converted into products or patenis. It is proposed, in future 
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research, to complement stages 7 and 8, to validate the proposed model. 

Innovations in projects of this nature can provide solutions to problems of great impact 
in older adults, for example, in the prevention of cardiac arrest, which according to the OMS 
(2020), is one of the two causes of mortality worldwide in the last 15 years, a problem that has 
increased within the health sector, since there are no systems or devices that could predict 


cardiac arrests in time. 
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